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In the present paper we shall be concerned with eigenfunction expansions 
associated with the Schroedinger operator — A+ q(x), where A denotes the 
3-dimensional Laplacian and q(x) is a real-valued potential function defined on 
the whole 3-dimensional Euclidean space E = £?, which tends to 0 at infinity. 
Solving the expansion problem, we go forward to apply the results obtained to 
clarify some properties of the spectrum of the Schroedinger operator and to 
show the unitary character of the S-matrix that plays an important part in 
describing the scattering process in quantum mechanics. 
Arch. Rational Mech. Anal., Vol. 5 
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to 


So far as ordinary differential equations are concerned, the eigenfunction 
expansion theory that originated from a study of WEYL [36] has been developed 
to a satisfactory extent by the efforts of Stone [31], TITCHMARSH [34], KODAIRA 
[19], Yosipa [38], Krein, Levitan and others 2 

As for partial differential equations, however, it seems that no complete 
theory, comparable with the one for ordinary differential equations, has been 
presented. One of the main difficulties appears to consist in the presence of 
infinite multiplicities of the spectrum, which is not the case with ordinary differ- 
ential equations. 

In 1934 CarLEMAN [6] studied the Schroedinger operator in E under the 
assumption that g(x) is locally square integrable, and obtained the following 
result: If the unique self-adjoint extension of the operator —A-+ q(x) exists, 
then there exist functions #(x, y; w), called spectral functions, such that for any 
/ (x) €L,(E) the expansion formula 


() f(x) = fd, J 0%, 954) £0) dy 


holds, where the formal derivatives with respect to uw of #(x, y; uw) satisfy the 
Schroedinger equation 


(2) —Agpt+ q(x) ~="y 


as functions of x and also of y. The same problem was investigated also by 
TITCHMARSH, mainly in the 2-dimensional case?. POVZNER [28] extended CARLE- 
MAN’S result to every self-adjoint extension of the operator —A-+q(x) and, 
moreover, using the Radon-Nikodym theorem, proved the existence of the 
spectral density o(u), by which #(x, y; jw) can be represented as 


d,O(x, 95 u) = p(x, v5 Mu) do(u). 
Thus (1) reduces to 


6) (0) = Sado Sve viw) fo) ay, 


where p(x, y; 4) serve as eigenfunctions, 7.e. they satisfy (2) as functions of x 
and also of y. But these eigenfunctions are not separated in the form of linear 
combinations of the products of the eigenfunctions of one space-variable, whereas 
this is the case with the ordinary differential equations. Such a separation of 
y (x, 5 4) was effected by MAUTNER [27] and GARDING [9] in the form 


co Nu 
(4) Ma) =S 2 v(x, mw) do(u) Se(y.u)f() dy, 


where y, (x, 4) are eigenfunctions (N,, Sov). 


Other approaches to eigenfunction expansions connected with more general 
partial differential operators were made by MAUTNER [21], GARDING [9], BRow- 
DER (2, 3, 4], GEL’FAND and KostyucHENKO [10], BEREZANSKI [1], FLEKSER [8] 


and Ito [73]. Their results can be summed up i 
S in the f 
the Schroedinger operator in E. ‘ iter cc 


Woe, 
3 a Snes (34, 35], NAIMARK [27] and the references cited in these books 
ee [35]. It is also abundant in references up to about 1956. 
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We should remark here that these authors have not specified the eigen- 
functions needed for expansion. In other words, a prescription is wanting for 
the conditions under which we should solve (2) in order to obtain the eigenfunctions 
y, (x, #) in question. Another point to be noted is that in (4) a countable number 
of eigenfunctions associated with a given spectral point are used, while a con- 
tinuous family of such functions is also known to exist in simple examples. 


We shall consider the case where q(x) tends to 0 at infinity. More precisely, 
we assume that g(x) is Hoelder continuous except for a finite number of singu- 
larities, is square integrable and behaves like O(|x|~?~") (h>0) as |x| 00. 
Many an important problem of potential scattering is included in this case. 
Now we can find bounded solutions p(x, k) of the Schroedinger equation (2) 
for 4>0, where k denotes a 3-dimensional wave vector such that |k|?=, as 
unique solutions of the integral equation 


7 |k| |~—y| 
(5) p(x k= eths— 1 f g(x) oy, b) dy, 


4m (a — | 
E 


k-x denoting the scalar product of k and x. g(x, k) represents the distorted 
plane wave, 1.e. the plane wave plus the outgoing scattered wave. On the other 
hand, in case of negative yw, (2) is solvable only for particular values of wv, called 
eigenvalues, which we denote by mw, (n=1,2,...). The associated solutions 
~, (x) CL, (E) are called ezgenvectors, to be distinguished from the eigenfunctions 
p(x, k)’. Here we agree to count mw, as many times as its multiplicity if mw, is 
degenerate. The @,,(x) can be regarded as forming an orthonormal system. In 
terms of the eigenfunctions p(x, k) and the eigenvectors q,,(x) our expansion 
formula? reads: 


(6) f(x) = (2)~! f v(x, R) f(y dk+ di, Pn (%), 


M n= 
where 


/(k) = (22) 4 f g(x, #) fa) ax, ( =f pn (2) f(x) dx 


and M is the 3-dimensional space formed by vectors k, which is not essentially 
different from E. Formula (6) is a natural generalization of the ordinary Fourier 
expansion, whereas (4) is not. (6) shows that the system of functions 9 (x, k) 
and @,,(x) is complete, but it will also be shown that these functions form an 
orthonormal set (in a sense to be specified later). 


Now we shall turn to the problem concerning the S-matvix that is known 
as a physical quantity which governs a scattering process®. The existence of 
the S-matrix, however, is not self-evident from the mathematical point of view. 
In order to study the S-matrix S, it is convienent to introduce the isometric opera- 


tors W,, called wave operators, by W..= lim e''” e~'!*>, where H and Hy are, in 


t— +0co 


our case, —A-+g(x) and —A respectively. S can be defined in terms of W, 
as S=W-* W_ and is generally believed to be unitary. 


8 We agree to call (x, k) an eigenfunction belonging to the eigennumber |k|. 
Cf. also SHNOL’ [30]. 

* For more precise expression see Theorem 5 (§ 8). 

5 See e.g. MOLLER [26]. 


1* 
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A mathematical formulation of the S-matrix theory has been proposed by 
Jaucu [14] and Kuropa [20]. According to their results, the existence and the 
unitary character of S will be cleared up if we prove the existence of W, and 
that W,. have one and the same range. The existence problem of the wave opera- 
tors has been investigated by Cook [7], JAucH & ZinNEs [15] and Hack [//]. 
Cook assumed that g(x) €L,(E). JaucH & ZinnEs treated a spherically symmetric 
potential: g(x) =|x|~° (1<c<#). Hacx considered the case where q(x) is locally 
square integrable and g(x)=O(|x|~°) (c>1) at infinity. No results, however, 
have been reported by them on the unitary property of S. Recently Kuropa [20] 
has solved both the problems under the assumption that g(x) €L,(E)°L,(£). 
We shall give a proof that S is unitary under the conditions stated above; these 
conditions are weaker than Kuropa’s assumption in one respect. Our proof 
depends partly on the so-called time-dependent theory and partly on the eigen- 
function expansion. 

We shall outline here the contents of the present paper. In I (§§1, 2, 3,4 
and 5) we give some relations between the resolvent of H and its kernel function, 
whose conjugate Fourier transform is deeply connected with the eigenfunction 
p(x, k), and introduce the kernel equations and investigate their properties. 
II (§§ 6—9) deals with the eigenfunction expansion. §6 is of a preparatory 
character, where the eigenfunctions g(x, k), tools for expansion, will be intro- 
duced. In §7 we shall comment on some properties of the spectrum of the 
Schroedinger operator and prove that the negative part of the spectrum consists 
only of discrete eigenvalues. In §§ 8 and 9 we shall state the expansion theorem 
with a partial proof thereof and show the absolute continuity of the positive 
part of the spectrum. In III (§§10 and 11) we shall show that S-matrix is 
unitary and at the same time complete the proof of the expansion theorem 
given in §8. IV is reserved for remarks on the 2-dimensional and _higher- 
dimensional cases. 


I. Resolvent kernel 

§ 1. Assumptions. The resolvent kernel. We shall consider the Schroedinger 
operator —A-+q(x) with the potential function g(x) defined on E=E%, where 
x denotes a position vector in E, its length being |x|. Throughout the present 
paper q(x) is assumed to satisfy the following conditions: 

(A) q(x) ts a@ real-valued function which is locally Hoelder continuous except 
at a finite number of singularities. Furthermore, q(x) is square integrable (9(x)EL,(B)) 
and behaves like O(\x|~*~") (h>0) at infinity, i.e. there exist positive numbers h, 


Cy and Ry such that 
: a a [ata lee Col ei toy | ade, 


Let us first define the operator A by A f(x) = —Af(x)-+q(x)f(x) for f€ Cote 
where Co°(E) consists of all functions which are infinitely differentiable and 
have compact carriers. Cj°(E) is contained and dense in the Hilbert space L,(E) 
in the sense of the L,-norm || || (|| || will be used exclusively for the L,-norm). 
Then it is known® that, under the conditions (A) imposed on q(x), A is lower 
semi-bounded and essentially self-adjoint in L,(E). Moreover, if we denote by 


g See Kato [16], StuMMEL [32] and WrenHoLtTs [37]. Our conditions (A) are 
more stringent than required for essential self-adjointness. 
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H the unique self-adjoint extension of A and by H, the corresponding operator 
for the case g(x) =0, then it is known? that D(H) =D/(H,). 

We shall study the resolvent R,=(H—A)+ of H and its kernel G(x, y; 2) 
for non-real 4; later this kernel will be extended to real 2 outside the point 
spectrum of H. The main results are summarized in 


Theorem 18, Let ImA=-0%. i) R, ts an integral operator of the Carleman 
type’ and its kernel G(x, y; 2), called the resolvent kernel, satisfies the integral 
equation 


oiVAlx—yl 


-_ 1 eil4le— ie 4 
Pe ree ge aie ee Ne 
E 


4n|x— 


as a function of x ae. in E for ae. fixed yCE. ii) If F(x, y; 4) is a solution 
of (1.1) such that F(-, y; A) €L,(E) for each fixed y, then F(x, y; A) is the resolvent 
kernel of R,: F(x, y; A)=G(x, y; A) forae. x, yin EXE. 

We shall call (1.1) the kernel equation. 


Now let B be the Banach space of all continuous functions u(x) defined 
on E, tending uniformly to 0 as |%|—> 0, with the norm 


ell» = max |u|" 
We also put for Imx>0 dsm 
ei %l*— iy 


ike= 1) 5 dz 
AM (x, y; 2) Malia z| 2 (z, isao) ) 
2) 
(P= 40254) 
Then we have the following 
Theorem 2. A(x, y; x) are continuous in x for x+y and A“ (-, y;x)EB 
(Imx=0). Further, let Imx>0 and Imx?+0. Then there exists a umque solution 
)(x, y; x) in B of the modified kernel sai 
(4) re ee AC 5 
eee ate G ial AO yeu abs : 


ei ¥le— i 


H(z, y;x) dz 


xa] 4 


7 See Karo [16]. We denote by D(T) the domain of the operator T. 
8 For the proof of Theorem 1, i) it suffices to assume only that qg(x)€L,(F). 
® Im A means the imaginary part of /. 
10 An integral operator T as well as its kernel T(x, y) is said to be of the Carleman 
type if T(¥, y) satisfies the conditions {| T(x, y)|?dy < oo fora.e. * and T(x, y)|?dx< 00 
E 


for a.e. y. See also STONE [31]. PovzNErR [28] showed that if g(¥) is continuous 
on E, then the resolvent R, of every self-adjoint extension of A is of the Carleman 
type. This result was extended by Burnat [4] so as to include the case where g (*) 
has a countable number of singularities x,,, in the neighbourhoods of which 


q(¥) =O (|x — %,)—) (Cn <2). 

1 Here and in the sequel by |/A is meant the branch of the square root of A with 
Im i> 0. 

12 “qe,” means ‘“‘almost every”’ or ‘almost everywhere”. 

13 Our Banach space is somewhat different from the one ee by PovzNER 
[28]. He used the Banach space B, (v being fixed) of all functions u (7) = =al4—y|4- 
u(#), where a is a constant and v(* ) is continuous and tends a rota e 0 as |x| + oo, 
with the norm ||u||z,=|a| + max |v(x)|. Cf. also Burnart [5]. 
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for each fixed yCE. If we put 


3 
(1.4) H(x, y; x) = A(x, ys x) + HO (x, y; x), 


1=0 


then H(x, y;*)=G(x, y; x) for ae. x, y mn EXE and, moreover, 
(1.5) Ped Cree eae all Sta 


where C is a constant dependent on x but not on yCE. 

As a consequence of Theorem 2, we can regard the resolvent kernel G(x, y; A) = 
Coy: \/4) as a function of x defined everywhere in E except at x=y for each 
fixed y. 

§ 2. Proof of Theorem 1. The kernel equation. Let us first show that R, 
is an integral operator of the Carleman type. Denoting by R,, the resolvent 
of H, and by V the operator of multiplication by g(x), we have the following 


operator equations (the so-called second resolvent equations 4) in virtue of the 
relation D(H)=D(H,): 


Ry, is known to be an integral operator of the Carleman type with the kernel 
(42|x%—y|) 1 eV4lx—yl, Let us put /A=a+ib (65>0). Then we have 


eilaix—y| 2 
[f 2 Sateor 
EE 


here and henceforth we denote by C any constant, not necessarily the same 
(C in the second member of the above equality is a pure constant and C in the 
third depends on bd). This implies that V R,, is a completely continuous operator 
of the Hilbert-Schmidt type. R, being a bounded operator, R,V Ry, is also of 
the Hilbert-Schmidt type and should be an integral operator with a Hilbert- 
Schmidt kernel}, say K(x, v; A). Thus we have 


—2b|x—y| 


dxdy=C f \g@)pax f © dy=C|lq|P; 
E E 


Ga y|? 


R,VRoal(x) =f K(x, yi Aty)dy ae., 
E 
where 


Sf |K(x, 9; Pdxdy <o. 
EE 


Since Ry, and R,VR,, a fortiori by the above argument, are of the Carleman 
type, we see from (2.1) that R, is of the same type, too, and is representable as 


(2.2) Rit(s) =f G(x, 9; Afy)dy ae. 
E 
by means of the resolvent kernel G(x, y; A), where 


S|G(%, 9; a)Pdy<oo for ae. x, Jf |G(%, v3 a)Padx<oco for ae. y. 
E E 

# See HItLe & PuILiips [12], p. 126. 

1° See e.g. Riesz & Nacy [29], Chapter VI. 
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Let us proceed to the deduction of (1.1). For any f(x)€L,(E) we obtain 
from (2.1) and (2.2) the equation 


[ews aoa f Toray 


A 4 = 
(2.3) ; 7 |~—y| 


; { ; eilAle—yl 
— Gq | Olea Malas f iy ay ake: 
E E 


By repeated use of SCHWARz’s inequality we have 


Nee 2a a 


ei VAlz—yI 


Tea 10] 49 S CMG «Alltel l/l < ee 


fOrea.es 4. We can, ea interchange the order of integration in the last 
term of (2.3) by Fusini’s theorem. In view of the arbitrariness of {/CL,(E), 
it follows that 


eiVAlx—y ei l4le— y| 
GAM A) Piel eA pe Oa I ae | — (the 


for a.e. y with a.e. fixed x. By use of Lemma1.1 just below we can replace 
G(x, y; Aa) by G(y, x; Aa) and G(x, z; 4) by G(z, x; 4) and, on interchanging x 
and y, we obtain (1.1). Thus we have proved i) of Theorem 1. 


Lemma 2.1. G(x, y; A) ts symmetric in x and y: 
(2.4) Gen AG ve Ale OF Ase. iy 07 ot, 


ae Let us introduce an indefinite “inner product” <,> in L,(E) b 
a a x) g(x) =(f, g). Then what we have to show is that A, is symmetric 


with Fe ack to this inner product, 7.e. <R,f, g> =<f, R,g>. For this purpose it 
suffices to show that <H/, g> =</, Hg>. But this is obvious, since we have 
Hf, E> —=<Aof, EMV 1) = (ot 8) + Fex= Ui, Hos) + Vero =< Hog> 
+<f,Vg> =<}, Hg), noting that Hy is a real operator and that V is an operator 
of multiplying by ¢ (x). 

Let us now prove ii) of Theorem1. To this end consider the difference 
f(x) =F (x, v; A) —G(x, y; 4) which is by assumption in L,(E) (for a.e. fixed y) 
and satisfies the equation 


(2.5) fips ate” [ eae) fe) de a.e. 


The kernel (42|x — |) e*/4I*—4l q(z) is of the Hilbert-Schmidt type (see above) 
and hence defines a completely continuous operator 7, which is an extension 
of R),V. Therefore, we can rewrite (2.5) in the form 


(2.6) Maia ores le i= (== 1, 2; aa) 


where 7” are also integral operators of the Hilbert-Schmidt type. Consequently 
we have, noting (1.2), 


f(a)=— en z; VA)q dz. 
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The function A® (x, oe A) is bounded in x and z, as will be shown in Lemma 2.2 
that follows. This fact together with ScHWARz’s inequality gives 


2) SCS ale) 1 |42<C lal I 


Thus f(x) is bounded and hence V/€L,(E), which permits us to rewrite (2.6) 
as 


(2.7) f=—RaV{ED(A). 

Application of Hy — A to both sides of (2.7) yields 

(2.8) (A, -Afs=—Ui on, His sAf. 

Since H is self-adjoint and A is non-real, we have /=0, which proves ii) of 


Theorem 1. 


Lemma 2.2. I/ Imx=0, then A®)(x, y; x) and A(x, y;%) are bounded in 
x, y and x. 


Proof. Let us show first that A(x, y;x)=O(|x—y|*). By means of 
SCHWARz’S inequality we have for |x — y|=d>0 


dz 
Al(x, yx SC [oa 2—C/d 
|A1 (x, 9; )| [x—2]? |2—9]? lI 2 
E 
in consideration of homogeneity of the above integral in x—y. Thus we have 


(2.9) JAM (a, ¥; 6) |= Cle — yl 4, 


where C depends neither on x, y nor on x. 


Next we show that A®)(x, y; #) =O (|log|x—y||). We have 


JAC yal] C/ fee 2) p AM fans be) de +f \as4s- 

K (Ry) K (Ry)’ 

Here and in the sequal K(x, R) denotes the sphere of radius R with its centre 
at «(K(R)=K(0, R)) and K(x, R)'=E—K(x, R) the exterior of Kp. (For Ry 
see (A) of §1). We also denote by K(R)K(R,) the intersection of k(R) and 
K(R,). We shall now estimate J’ as follows: By (2.9) and (A) we get 


Fac it rapa ih Ja I 


a— Z| I 
B(Ry) AK (y | | K(Ry)’ —K(y, 1) 
1) 8° 


An estimation of J; is 


Pate ° dz ee dz 
J = | Rath Fe ca Pac po, a yy — 
ih) 0 | 4 ally |? 13 


W 
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For Jz we have, introducing spherical coordinates, 


vsin ddr db 
J aoe |z me =cff pth 


Reh (\x]2 +72—2| x] 7c0s 8)! 


= chef 7! (|x| +7 —||x| —r|) dr. 


If |x| <R, and x= 0, we obtain 
YL CN eh Os 


Ry 


the case x=0 can be dealt with separately, and a similar result follows. 


|x| > Ro 


|x| 


De eter er flelr Pars OC; 


Thus J; is always bounded; hence so is ]’. 


Now we have to estimate /: 


K(Ro) 
, sane tS )=hth (6=|x— yh), 
K(Ry) 0 K(x, 6/2) K(R,)—K(x, 6/2) 
d 
fisess dh aa anes 


K(x, 6/2) 


If 6=4, we have for J, 


hee | poaty=a ayaa] = salle a Cf ip=a ‘lel dete 


If 


*, }) K(Ry) 
If 6<, we first assume that |*|=R)+1. Then 
TG pants C: 
K(Ry 
For |«|<R)+1, on translating the origin, we have 
R+1 
{RRO Apeeaes re ih r134(r + 6 — |r — 6|) dr SC log 4}. 
K(2. Ry +1)—K(0/2) oe 

Thus we get 

eGo elec ho [aay as 
(2.10) 


SC|log|x—y|| if |x—y|<2 


where C is independent of x, y and x. 
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Now we shall show that A®)(x, y;x) is bounded. We have from (2.10) 


|A® (x, ys] <C fie ar -Hogle— yl] 24 Chia Og os 


K(y,4) 


L and L’ can be estimated as follows: 


Lc] [eee fie atta, 


Z| 2 
K(y, 4) 


Bede eae ie fe ee ae 
[re ad flat ra) Bai Op 


ssf ae Pe : a (See Jz above.) 


ws 


L,35C 


These estimates show that 
(2.11) |A@ (x,y; g/ SC, 
where C is independent of x, y and x. 


From (2.11) and the estimation performed for L’ we get 


(2.12) 


A® (x53 2)| SC, 


C being independent of x, y and x. (2.11) and (2.12) prove the assertion of 
Lemma 2.2. 

The following lemma will be needed for the proof of Theorem 2 that we 
shall give in § 5. 


Lemma 2.3. If Imx>0, then ||AM(-, y;x)||SC (¢=0, 1, 2,3, 4), where C 
depends on x but not on y. 


y; x)|| SC, where C does not depend 
on y. We assume that this is the case for A“ (-, y; x) and prove that the assertion 
is true for A“+4)(., y; x). 

Using the operator T defined before (see (2.8) and (2.9)), ™ we have Alt) (., y; x) 
=T A" (-, y; ), whence follows 


(2.13) AMP, 95 |S CI A(, ys |] SC 
in view of the boundedness of T. 


ix | = | 
§ 3. Asymptotic behavior of the function ae LOY, ) dy. In 
4 ae y)ay 
Be |a—y | 


order to solve the kernel and modified kernel equations, we shall study some 
asymptotic properties of functions of the form 


{ ei * |v — | 
ere ep u(y)dy (Imx = 0). 
E 


Expansions for the Theory of Scattering M4 


Lemma 3.1. Let Imx=0, let v(x) be locally integrable and v(x) =O (\x|~?~*) 
(|x|—> co) (e€>0). Then as |x| 00 
(3.4) p(x) =O(lx|7) +0 (|x|~*). 

Proof. Let R be such that |v(x)|SC|«|~?~* for |x| =R. If |x| >R we have 


1 ei *l*—y| r 
ee ee ie 


K(R) K(RY 


Y pOue fain bdr an 
I< LOM ay<sc | y? sin 
oo |~—y| ihe ye 


p(x) = 


|a|?+72—2 |x| cos 9)! 
=C|x|4 fr-*-*([x| +r —||x| —7|)dr=Clx[2+C]x|-* 
R 


In order to give an estimate for J, it suffices to note that v(x) is integrable over 
K(R) and that |y—y|=|*|—R. Then we have J=O(|x|+). By putting J and 
I’ together, (3.1) follows. 

Lemma 3.216, Let v(x)=O(|x|~3~*) (|x| 0) (e>.0), v(x) CL,(E) and let 
x=a be real. Then we have 


ef lA 


aa 4a|x| 
E 


(3.2) p(x) = e*402Yy(y)dy +O(|x|2-*?) +O (|x|), 
where w, denotes the unit vector with the direction of x and w,- y denotes the scalar 
product of w, and y. 

Proof. Let R, be fixed so that |v(x)|<C|«|~°~* for |x| 2R, and let |x| be 
so large that |«|#=R>R,. Then we have 


ta\|x—y| . 
(33) v(x) = fal é oopayt f\ar+r, 
K(R) 


4n |a—y| 
K(R)’ 
An argument similar to the one used in proving Lemma 3.1 gives 
eS Pg Oa £ ea Mame — a On Wy a ao la 
We proceed to estimate J. Considering | y|<R we have 
poh ie falta ay) | (n= Olly |/lxl)2), 


Je yl =|2b [t+ (oe) +m (a= Olly I/l1)), 


and hence 
ta|x—y| : : : : : 7 
aged & ae ei @lx|—idog-ytials|m ae "2 ei @le—y fe ee WO, ¥ ef ale—yl, 
[v—y| 4 || |7| 


16 Cf, POVZNER [28]. He obtained the result in which the last two terms in (3.2) 


é 
are replaced by O (jaj777 Bae i Our result enables us to proceed along his line without 
assuming that the potential has the asymptotic form g(¥v) =O(|¥|~3°~*), which was 
assumed in [28] together with the continuous differentiability of q¢(-). 
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Thus we have 


ta|x| 


= t oi LA Wy+V d é fetsveto(y) dy + 
If 4 m|%| fe v(y ) er “4n|4 [ 
E K(R)’ 
ta|x ; : 
Gy bs ath Comeau es PLS 
4n|¥| . 
K(R) 
ont ee ia|x—y| pees ieee [ o-: etl*—"y(y\dy 
4 |x| [me u(y) ay, 4n|x|? J 4 y 
K(R) K(k) 
‘ala 
=— fo fetter oo)dyththt hth 
4n|x| : 


I; ((=1, 2, 3, 4) are estimated as follows: 


|Ty| Chal [u(y dy SC |x| AR Clan), 
K(k)’ 


als Chel fe belle aay =Clor*( fF eo Jae | 


K(R,) —K(R 


ONC = seme 2R2-# =. C |x|? 46 ]x|-- 


aleck fe kody=cll?( f brleoay+ f 
oe K(Rj) K(R)—K(R,) 
SC|xf P+ Clx[ PR 8 = Chale Caf, 
Hal Clot? f lleorlay=Clxl9(_f 


K(R,) K(R)—K(R, ra 


Clair +.Clalc* Rey Gla Ae 
In these estimates the constants C may or may not depend on a, € and v(-). 
(3.2) now follows from (3.3) and the estimates for J and I’. 


Lemma 3.3. Let v(x) and p(x) be as in Lemma 3.2. Then (x) satisfies the 
vadvation condition 17, i.e. 


(3.4) lim | cet aq 


|x| co 


S=(06 


For the proof see PovzNER [28], Chapter II, Lemma 2 


Lemma 3.4. Let p(x) €B be a solution of the integral equation 


, is 1 + eile —y| d 

(3.5) P= — Ze | yp 70) eWay, 
E 

where ars veal. Then 

(3.0) fem" g(x) p(x)dx=0, 


E 


where w ts an arbitrary unit vector. 


There are definitions of the radiation condition other than the one adopted 
here. See e.g. MIRANKER [22] and MUELLER [24]. Here we follow PovzNeR [28]. 
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For the proof see PovzNneR [28], Chapter II, Lemma 5. But it seems necessary 
to add some remarks. The proof relies on Lemmas 3.2 and 3.3 and on GREEN’S 
formula. We maintain that if y(x)CB satisfies (3.5), then p(x) =O (|x|~4) (|x| >0x), 
so that Lemmas 3.2 and 3.3 may be applied. In fact w(x) =O (|x|) +0 (|x|~") 
by Lemma 3.1 because g(x) =O(|x|~?~”). Hence repeated application of Lemma 
3.1 to (3.5) furnishes us with the result p(x) =O (|x|7). 

Another remark applies to the use of GREEN’s formula. An explicit use of 
the Hoelder continuity of q(x) and q(x) is made to obtain from (3.5) the dif- 
ferential equation 


(3.7) —Ag (x) +4(*) v(x) =a@ p(x), 


which holds for x different from the singularities of g(x) (Hoelder continuity 
of v(x) follows from (3.5) and the proof of Lemma 4.1 given below). We apply 
GREEN’S formula to g and @, which also satisfies (3.7), and get 


(3.8) 0= [ (4p P—9-Apax= [O22 — 9 \as 
S(R) 


K(R) 


where S(R) is the surface of K(k) and n is the outer normal. (3.8) is valid if 
K(R) contains no singularity of g(x), but even if singularities of g(x) exist in 
K(k), it is seen to be true. For we first exclude from K(k) the singularities 
which are finite in number, and then by a limiting process we obtain (3.8) valid for 
K(R), when we note that dy/én as well as 6G/0n is O (|x — x,|~4) («,: singularities 
of g(x)), as is seen from differentiation of (3.5) under the integral sign, so that 
the surface integrals taken around the singularities are arbitrarily small. 

The above remarks are enough to get to (3.6), following the proof by Povz- 
NER [28]. 


§ 4. Solution of the kernel equation. Let us define an operator 7 for func- 
tions in B by 


Eee ed fee | 


7 AS aa a 
E 


ina —y| 


7] q(y) f(y)ay (Imz 20, fEB). 


Lemma 4.1. TJ, is a bounded linear operator on B to B (Imx=0). 

Proof. Let f€ B. Then q(x) f(x) =O(|x|~?~") (|x| ce), so Lemma 3.1 shows 
that JT, f(x) +0 uniformly as |x|— o. 

T, f(x) is bounded in x. In fact we obtain from (4.1) 


IZA) SC Ile f LE ay. 


The boundedness in x of the right member follows from the estimation carried 
out for L’ in the proof of Lemma 2.2. 

Next we have to show the continuity in x of Z/f(x), and for this purpose 
we consider the difference 

in|x—y| tn\|x’ —y| 
1 e =€ 
a HS Conese : Die 
T.f(z) — Tf) =— [ ——-== a tty) ay 


5 E 
Bee ae ee 
DV 4Agdole—y) l= 7 
E 


él lg(y) {)dy=h +h 
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Considering the inequalities |e’*'*~*! — e'*!”~9!| < |x| Jw — x’| and 
sabe aes | = eer e 
[ol #91 | = Tole —9 


we can give estimates for /, and J): 


Alsi — x] f HOLL ay sclifilsx —*'|, 
E 


y 


|Jo] SC lx — x’ 


lq (y)| Tere es 

frase = se oIh Fe 

Iflle lx—y| |x —y| Neale | 
E 

as is seen from the estimations worked out for L’ and A(x, y; x) in the proof 

of Lemma 2.2. Hence it follows that 


(4.3) T, f(x) — Ef x’) =O (|x — x’ |8) for |x — x’|-0. 


(4.3) shows that Z/(x) is not merely continuous, but also Hoelder continuous 
in x. Consequently 7f turns out to be an element of B. This completes the 
proof of the lemma. 


Remark. The proof given above uses only the boundedness of f(x). Thus 
T,f belongs to B and is even Hoelder continuous if / is only bounded. 


Lemma 4.2. TZ) is completely continuous }°. 


Proof. Let }, (n=1,2,...) be any sequence in B such that ||f,||,; 1. We 
have to show that we can take out of the sequence u,— Tf, a subsequence 
convergent in the norm of B. 


Re-examining the proof of Lemma 4.1, we find out firstly that ||w,||p are 
uniformly bounded with respect to m, because ||f,,||z, <1, and secondly that on 
the same account the continuity of w,,(«) is uniform with respect to 7. Moreover, 
we see that w,(x)—>0 uniformly in as |x|—> oo. Thus u, (x) are equi-bounded 
and uniformly equi-continuous with respect to » in any compact domain of E. 
On applying the Ascoli-Arzela theorem to wu, (x), we can choose a subsequence 
u,(x) converging uniformly to a continuous function w(x) in any compact 
domain. Since u,,(*%)—>0 uniformly in n’ as |x| 00, u(x) also tends to 0 at 
infinity, and hence we have ||w,,—||,—>0, which was to be proved. 


We are now in a position to make use of the Riesz-Schauder theory of com- 
pletely continuous operators in a Banach space ¥*. If T is a completely continuous 
operator in B, the equation f=g+ Tf is solvable for any given gCB, if and 
only if = Tf implies that f/=0. Thus we have the following 


Lemma 4.3. Let g€B. Then the integral equation 


ei%|*—9| 


(4.4) i) =e) — fF | a tT o)ay me 20 


8 Cf. also PovzNER [28], Chapter II, Lemma 3. 
19 See e.g. Riesz & Nacy [29], Chapter VI or Yosipa (39) S44: 
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has a umque solution in B if and only if the homogeneous equation 


tn|x—y| 
(4.5) f(x) =— ff —— a0) i 0) ay 
has the unique solution { (x) =0. 


Lemma 4.4. If x? is real and positive, then (4.5) has the unique solution f (x) =0. 


Proof. If f(x) satisfies (4.5), (3.7) follows from the assumed Hoelder continutiy 
of g(x) and that of / (x), remarked in the proof of Lemma 4.1 (see also the remarks 
to Lemma 3.4). On the other hand, Lemmas 3.2 and 3.4 give the asymptotic 
order 


(4.7) f(x) =o(|x|7) (|x| > 0). 


According to Kato [18], however, any solution of (3.7) subject to the condition 
(4.7) vanishes indentically at least outside a sufficiently large sphere containing 
all the singularities of g(x). Hence the unique continuation theorem for solution 
of an elliptic differential equation?® is applied to obtain the result that /(«) = 
everywhere. 


Lemma 4.5. Let Imx>0. Then (4.5) has non-trivial solutions in B if and 
only if x? 1s an eigenvalue of H. In particular, if Imx? +0 in addition, (4.5) has 
only the trivial solution f(x) =0. 

Proof. Let {(x)€B satisfy, (4.5). Since g(x) €L,(E), g(x) f(x) EL,(£). Then 
from (4.5) we have /=—Ry,:Vf, which is nothing but (2. a with A=x?, and 
hence we get /=0 if x? is non-real, in exactly the same way as in § 2. This proves 
the second assertion. 

Let /(x) be a non-trivial solution of (4.5) (Imx>0). Then, as we have seen 
before (the proof of Theorem 1, ii) in § 2), it follows that /, V/€L,(E) and Hf=#*/. 
This proves the necessity of the first assertion. 

Conversely let x? be an eigenvalue of H with an eigenvector /€ D(H) =D/(H,). 
Then V{€L,(E£) and (4.5) readily follows. Accordingly, by use of SCHWARz’S 
inequality, we have 


POSS (ler ee a9 f Ia y) /oikay) = ()=Im»), 


from which and (4.5) we see that /€B, as in the proof of Lemma 4.1. This 
proves the sufficiency. 


Lemma 4.6. J(Imx=0) depends continuously on x; i.e. given any e>0, 
there exists a d=6(e) such that ||T,, —T.,||p <7) if |%,—| <0. 

For the proof see PovzNER [28], Chapter II, Lemma 8. We have adopted a 
norm different from the one used by PovzNER (see footnote 13), but this makes 
the proof of the present lemma rather simpler, for we have avoided inclusion 
of singular functions in B. 

20 See e.g. MUELLER [25]. 


21 The norm of an operator T in B is defined by ||7||z = sup ||7/llz- 


Ilflla=1 
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All the bounded operators on B form a Banach algebra B, in which ||Z, — T ||, > 0 
(n> cc) and the existence of Z7+ and T+¢€B imply ||JZ,*— T7||,+0*. Thus 
from Lemma 4.6 follows 

Lemma 4.7. I/ x2 is not a non-positive eigenvalue of H, (I—T)1€B exists 
and depends continuously on x in the sense of the norm || ||p- 

Lemma 4.8. Let D be a compact domain of the upper x-plane (ImxZ0) such 
that D does not contain the square roots of the non-positive eigenvalues. Let fc B 
be a unique solution of 


(4.9) fe= bx t Ly, 


where g,€B is strongly continuous in xCD. Then f€B is strongly continuous 
innxeD. 

Proof. Since f,=(I— T,)4g, and (I— JZ) + is uniformly continuous in x € D 
(Lemma 4.7) and also g, is by assumption strongly continuous in x, the assertion 
follows immediately. 

§ 5. Proof of Theorem 2. We shall show first that A“ (x, y; x) are continuous 
in x unless x=y. By definition we have for 7=0, 14, 2, 3 


AWD (x, y; 2) — AGED (x’, y; 2) 


ve 7 eit lx—2| eit i% =2| 
ie fg 
13) 
1 1 1 ix|a’—z (Gi 
qa f ap eae Oh 4°] ale 
E 


=f,+ Je. 


Proceeding as in the proof of Lemma 4.1, we arrive at the following inequalities 
for, [nang je: 


\A|<C|x — x iy |q(z) ee rN) Gk 


als Clea] f pee ee dx. 
E 


Since A(z, y; x) are bounded in z outside some sphere with its centre at y, 
as is seen in the proof of Lemma 2.2, we get 


[Als C|x—x’| Rhee: ge) d td ge ae * Jie )AM(z, y; x)| dz, 


K(x, 0) 
ze Giooes 
[Fal = lx 1 f eae =4 Be “| ear dz+ 


+= # [ \2@ JAG yale, 


K (x, 6)’ \ K(x’,n)’ 


*2 See any book on Banach algebras, e.g. HILLE & PHILtips [12] and Yosipa [39]. 
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where we have chosen 6 and 7 so that K(x, 6) and K(x’, 7) do not contain the 
point y. Then from the estimations in the proofs of Lemmas 2.3 and 4.1 we 
obtain by ScHwarkz’s inequality 


AlSClemx'|4e|e— a 
|Jo]S C|x—x' | + Clx—x' 


lal, 950 
b4Clx—a'|[fallAC, 95, 


which prove the asserted continuity of A" (x, v; x) ((=1, 2,3, 4), while A (x, y; x) 
is obviously continuous in x for x+y. 

That A“(-, y; x) €B follows from (2.11) and the remark after Lemma 4.1. 

The above facts together with Lemmas 4.3 and 4.5 show that a unique solution 
H (x, y; x) of the modified kernel equation exists for fixed y. If we prove that 
|| (., vy; x)|| <C, where C depends only on x, we can complete the proof of the 
theorem in view of Theorem 1, ii). 

From (1.4) we see that H)(x, y;) is the sum of A“)(x, y; x) and f(x)= 
—Ry 2~VH™ (x,y; x) for fixed y. It is easily seen that ||f|| < || Ro,2||||¢||||H@( ¥:%)|lp 
is independent of y, for we have 


|2C, 95 lle SE — Fy AC»: ale. 


Now in virtue of (2.12) we have proved the desired result. 


II. Eigenfunction expansions 
§ 6. Existence of the eigenfunctions. In quantum-mechanical problems of 
potential scattering, the distorted plane wave plays an important role. It is 
generally assumed to be a bounded solution of the Schroedinger equation 


(6.1) —A p(x) +9 (x) p(x) =|RP g(x) 


with the asymptotic conditions v (x) = g (x) — e***=0 (|x|) and |x| (du (x)/4|x| — 
i|k|v(x)) + 0 (radiation condition) for |«|—> oo. Here k is a 3-dimensional vector 
called the wave vector and |k| denotes its length. We shall denote by M the 
totality of the wave vectors k, which is, of course, isomorphic to E. g(x) has 
the form (plane wave) + (outgoing wave). It is also possible to consider the 
distorted plane wave which has the form (plane wave) + (incoming wave) in 
quite a similar way. As is described in most of the physical literature ?8, the 
distorted plane wave g(x) is obtained as a solution of the integral equation 


i{al|x—y| 
— pik-% 1 ea : is d 
(6.2) p(x, k) =e real ene Pay 
E 
which is more convenient for mathematically rigorous treatment than (6.1), for 


(6.1) is not valid at singularities of g(x). Let us introduce the function 


+ || |x—y| 
(6:3) p(k) =— ZL f Salyer. 
E 


|~—y| 
23 See e.g. Motr & Massey [23]. 
Arch, Rational Mech, Anal., Vol. 5 


bo 
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Then, if we put v(x, k) =@(x, k) —e'*'*, (6.2) becomes 


» t|kl|x—yl 
(0.4) v(x, R) = p(x, Rk) = a | oe =i g(y)v(y,k) dy. 
E 
Since p(-, k) € B as is easily seen from the remark after Lemma 4.1, (6.4) has 
a unique solution v(-,k)@B for each fixed k (\|k|>0) in virtue of Lemmas 4.3 
and 4.4. Moreover, we see from Lemma 4.1 and the remark thereafter that 
v(x, k) as well as p(x, k) is Hoelder continuous in x. Since g(x) is Hoelder con- 
tinuous, g(x, k) is then seen to satisfy (6.1). Thus we arrive at the following 


Theorem 3. There exists a unique solution v(x, k)€ B of (6.4) for |k| >0 and 
p(x, k)=e'® *+ v(x, k) is a solution of (6.1) and (6.2). Moreover, p(x, k) as well 
as v(x, k) is bounded and uniformly continuous in x and k for xCE and kED, 
where D is any compact domain of M not containing the origin. 

p(x, k) is called an eigenfunction associated with the ezgennumber |k|*, and 
we can speak of an co®-continuous family of eigenfunctions in conformity with 
various magnitudes and directions of the wave vectors k. However, we should 
remark here that under our assumptions (A) on g(x), v(x, k) does not always 
satisfy the asymptotic conditions mentioned above 4. 

Proof of Theorem 3. We have only to prove the last statement. From the 
remark after Lemma 4.1 it follows that |p(x, k)|<C, where C is independent 
of x and k. Since (6.4) can be written as v(-, k)=(1 —T},) 1 p(-, &), we see from 
Lemma 4.7 that v(x, k) is bounded in x¢€ FE and k€D. Similarly it follows from 
Lemma 4.8 that v(x, k) is uniformly continuous in x and k in the same domain, 
if we show that #(-, k) is strongly continuous in k CD. 


Given any 6>0, we have for a sufficiently large R=R(6) and for |x|=R 
(6.5) jot asc f 2" ay <a 
; j 


by Lemma 3.1. For |x| <R we have 


4 [Ril |e—y9|__ pt [Rel lay 


la(y)|a@y+ 


yl 


leek) pte asc flé 


5 BIR, ae ey 
ref 20h errata 


SC f[ |Iml—lAallla)lay +¢ ee Pore 
K(R,) Acie 
+C flab, y| [Ol ay +C feet ay 
K(R,) : K(R,)' 3, 
=e lee 


Here we can choose an R,=R,(d) sufficiently large so that each of J, and fy 
may be bounded by 6/4 (see the estimation of J; in the proof of Lemma 2.2). 
On the other hand, j, and /,; can be made arbitrarily small J1< 6/4, J5< 6/4) 


** These conditions are satisfied if we assume that q(x) SO(W| 4 *) (ale oo) 
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by taking jk; —k,| =||%,|—| || small. Together with (6.5) we see that 
|b (x, hy) — ple, be) |<d (BEE, hy, BM) 


if |k; —k,| is sufficiently small. This proves the desired continuity of p(-, k) 
and completes the proof of Theorem 3. 


§7. The spectrum of H. It is generally believed that the positive real line 
is occupied by what are called the continuous eigenvalues of H which are, in 
our terminology, ezgennumbers, while in the negative half there exist possibly 
discrete eigenvalues but no continuous spectrum. This is certainly the case 
under our assumptions (A) on g(x). But from the mathematical standpoint the 
non-existence of positive eigenvalues is not self-evident. This problem has been 
raised by PovZzNER [28], though he only commented on its plausibility and 
proved it for the case where q(x) is a function of | x| at least for large | «|. Recently, 
however, Kato [18] has proved a theorem®® which implies the non-existence of 
positive eigenvalues in our case. We shall give a proof of this fact later (Theorem 6). 
Concerning the features of the negative part of the spectrum of H, we have 


Theorem 4. On the negative real line the continuous spectrum of H is absent®®. 
The negative eigenvalues, if they exist, are of finite multiplicity and are discrete 
in the sense that they form an isolated set having no limit point other than the origin 0. 

Proof. Let a be sufficiently large so that we may have H+a=H,+V+a>a' 
>0. (Note that H is bounded from below; see § 1.) Since the spectrum of H, 
occupies the full interval [0, oo), that of (Hy+-a)7™ is the interval [0, 1/a]. Now 
the second resolvent equation gives (H+ a)1=(H,+a)1—(H+a)1V(H,+a)4, 
where (H+ a)1V(H,+ a) is a completely continuous operator, for —a lies 
in the resolvent sets of both H and H, and, as we have seen in § 2, V(H)+ <a) 
is of the Hilbert-Schmidt type. According to a well-known theorem of WEyYL27 
that asserts the invariance of the set of limit points?8 of the spectrum of a self- 
adjoint operator by the addition of a completely continuous operator, the spectrum 
of (H+ a) consists of the whole interval [0, 1/a] and, possibly, of isolated points 
outside [0, 1/a] (which may have limit point 1/a). This means that, on the 
negative real axis, H may have only discrete eigenvalues of finite multiplicity 
having, possibly, limit point at 0. This completes the proof of Theorem 6. 

Hereafter we shall denote by w,, the discrete negative eigenvalues counted 
according to their multiplicity and by ¢, the orthonormalized eigenvectors (see 
footnote 3) associated with u,,. 


§ 8. Expansion theorem. As we have pointed out in Introduction, we can 
use the distorted plane waves p(x, k) constructed in §6 as the eigenfunctions 
for expanding an arbitrary function. 


Theorem 5. Let f(x) be an arbitrary L,-function. i) Then the generalized 
Fourier transform 


(8.1) j(k) = (2A Lim. f p(B) f(a) dx 


25 Cf. Lemma 4.4. 

26 Cf. also PovzNER [28], Chapter II, Theorem 6. 

27 See Riesz & Nacy [29], p. 367. iu 

28 Eigenvalues of infinite multiplicity are included among the “‘limit points’’. 
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of f(x) exists and belongs to L,(M), where 1.i.m. es dx means the limit in the 


mean of the function [ ...dx as Noo. Also ae generalized Fourier coefficient 
K (N) 


(8.2) f.=J Pr() (x) dx = UF, Gn) 
E 
of f(x) exists, and s lin PB<oo. 
i 
ii) We have the following expansion formulas : 


(8.3) f(x) = (2a)~flim. f p(x, 2) f(h) dk +1 i m. Shee 


(8.4) (APR? =S ((ePak + s /,|2  (PARSEVAL’s equality), 
M 


(B55) hig) f F(R) &( k) dk + Si tfn (generalized PARSEVAL’S equality) . 
M wat 


n (8.5) g€L,(E), and in (8.3) lim. f...dk means 7 Lime if nears 


M a a PSO SEG) 

iii) Conversely let F(k) € L. 2(M) be given. Then we can construct an L.-function 
f(x) through formula (8.3) with f(k)= F(k), i =0 such that (8.1), (8.2), (8.4) and 
(8.5) hold good. 

iv) Let E,, be the resolution of the identity associated with H, and let P=I— Ey. 
Then the part in PL,(E) of H is unitarily equivalent to Hy. 

v) [€D(H)=D(A,) if and only if | ki? i( (K)€L,(M) we 2 fin lfn P <oo. We 
have the following representation of H: 


~ 


(8.6) Hi (x)= (2) (x, ) FB) dhe + Lim. Stig fy Pn (2): 


The proof of i), 11) and v) except for the complete domain characterization 
of H will be given in the following section; it seems appropriate to remark here 
that the proof depends on some properties of the spectrum of H (Theorem 4) 
but not on the absolute continuity of the positive part of the spectrum (Theorem 6), 
which will be established in the course of the proof of ii). We shall also make 
use of Lemma 4.4, which is almost equivalent to the non-existence of the positive 
eigenvalues of H, but not in the strict sense, for we have not shown that every 
eigenvector associated with any positive eigenvalues should satisfy the homo- 
geneous integral equation (4.7) and that it should be a B-function. Pertaining 
to the positive part of the spectrum of H, we have 


Theorem 6. There exists no positive eigenvalue of H. Moreover, the spectrum 
on the positive real line is absolutely continuous. 

The proof of Theorem 5, iii), iv) and the domain characterization of H will 
be given in connection with the proof, by the aid of the time-dependent theory, 
that the S-matrix is unitary. 

It is implied by (8.6) that the Schroedinger operator H=—A-+V admits 


of a diagonal representation, or H is diagonalizable, in terms of the eigenvectors 
7, (x) and the eigenfunctions p(x, k). 
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§ 9. Proof of the completeness of {p,,(a)} and {ep (x, k)}. We shall divide 
the discussions that follow into two parts. One concerns the projection P and 
the other E,. So to speak, we shall first take up the expansion problem for 
the function P/(x) and then extend our research to the function Eyf (x). 


Let us consider the conjugate Fourier transform g(x, k; #) of H(x, y; x) with 
the defining equation 


(9.1) ROC NCEE ah fehl) ek ydy 


for each x and Imx>0, Imx?=£0, noting that H(x, -; x) L,(E) (see Theorem 2 
and Lemma 2.1). Actually, however, “‘l.i.m.’’ in the above definition is un- 
necessary. This will be made clear by the following 


Lemma 9.1. // Imx>0 and Imx?=:0, then H(x, y; x) 7s absolutely integrable 
im x or y and g(x, k;x) ts a bounded function of x and k for xCE and kCM. 


Proof. According to the kernel equation (1.1) for H(x, y; x)=G(x, v; x?) we 
have (4=a-+71) 


[\Hy:alaxsc [ g 
E E 


—bl#—y| 


|x—y| 


x—3| 


—b|x—z2! 
dx 4 Cf |g) H(z, y;x)|dz PS dx. 
E E 


A further estimation is possible by (1.5) and SCHWARz’s inequality, which gives 
[\H(x ys *)|axSC+Cllg|[2. 95 OSC, 
E 


where C depends on x but not on y and is finite if Imx>0 and Imx?=-0 (see 
Theorem 2). By the symmetry stated in Lemma 2.1 H(x, y; x) is absolutely 
integrable in y, too, and hence it follows from (9.1) that g(x, k; x) is bounded 
in x andk for xC€E andkemM. 


Set 
(9.4)’ TNC tire (22)? (|k|2 — x?) g(x, hk; x) = esha) W) (x, k; x). 
Then we have 


Lemma 9.2. h® (x, k; x) is a B-function of x for ke M, Imx>0 and Imx? +0 
and satisfies the integral equation 


ina—y| 
(9.2) W(x, k; x) = p(x, k; x) — fH q(y) AY (y, k; x) dy, 


4n J |x—y| 
E 
where 


(9.3) p(x, iw) =— tf 


4n J |¥—y| 
E 


ix|x—y| 


q(y)e* dy 


is a B-function of x fork€ M and Imx=0. Moreover, we can extend the definition 
of h(x, k;x) to the case where Imx=0 and Rezx>0. h(x, k; x) 1s bounded and uni- 
formly continuous in x,k and x for xc E,kCM and Imz2Z0,0<aSRe Hs B00 
In particular, h(x, k;|k|)= (x, k). 
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Proof. Let us first show that g(x, k; x) satisfies the equation 


th-x 


1x |x—y| 
2g 41 é ay 5 
9.4) gle bie) = Qa) t oe af ala) aly Bama ys 
E 


From the kernel equation (1.1) we have for any /CL,(E) 


tulx—y| 
f Hays) fy) ay : [= f(y) dy — 
F 


~ aa J Tea 
E 


Lf a 90) 42 [ Heys) fo) ay, 
E 


E 


where we have performed interchange of the order of integrations, which is 
permitted since the last integral is absolutely convergent because f | | dz S\laysc 
E E 


liall {I (z, <5 |] Aldll SCllall l|fl| (by Theorem 2). Introducing the Fourier trans- 
form of f(x) by 
(9.5) fo (h)= (2az tam. | e-"* h(a) ax 


E 


and making use of PARSEVAL’s equality, we get therefore 


ee tk-x« ns 
fel Piste) fate) ie (2n)-4 [ re fol) Ph pee 
4 et %lx—2| % 
70 a were g(2, k; %) dz| fo (Fk) dk 
M 


J |% 
EB 


where we have used the fact that the conjugate Fourier transform of (42|x—y|)1 x 

e’*|*—% as function of y is (27)! (|k|2—x2)4 e'*:* and have interchanged the order 

of integrations in the last integral, noting that |x—z|1e™!*~‘l g(z)EL, and 

that f |g(z, 2; %)| |fo(k)|@k S C||A(z, -; x)|| ||f|| SC|l¢|| (by Theorem 2). Since 
M 


fy€L,(M) is arbitrary, (9.6) yields (9.4). Then the equation (9.2) for h(x, k; x) 
follows at once from (9.1)’ and (9.4). 

From the remark after Lemma 4.1 we can see that p(-,k;)@€ B and that 
p(x, k; x) is bounded for «EE, REM and Imx=0. Similarly, the second term 
on the right of (9.2) is also in B if Imz>0 and Imx?=-0 (see Lemma 9.1). Thus 
AY (., k; x) € B for Imx>0, Imx?=+0, so that we can rewrite (9.2) ash (., k; x)= 
(i —T,) 1p(-, k; x). Noting that p(-, k;|k|)=p(-, k)CB, this proves the lemma 
in view of Lemma 4.8 and Theorem 3, if we show that f(-, k; x) € B is strongly 
continuous uniformly in ke M and x for Imx=0. The uniform continuity in 
k can be proved by an argument similar to, but simpler than, the one used for 
proving the uniform continuity in x and k of w(x, k) in §6. Thus it remains 
to prove the uniform continuity in x of (-, k; x). 


Given any 0>0, we can choose an R=R(6) so that for |x]=R we have 


(9.7) lpi kewl f 2 or ay < 0/2 
; 
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/v virtue of Lemma 3.1. For |x| <R 


r—y| 


1, |x ix |x—y] 
Pe Rim) — ple kindle | —eitele9 


la(y)| ay 


|~—y| 


a vl et %2|¥— Pa eel a 
sc [la ayer cf [0 a 
[xy] 
K(R, 
a 
Again we can choose an R,= R, (6) so large that J’ may be bounded by 6/2 (see 


the estimation of J; in the proof of Lemma 2.2). In order to estimate J, we 
first note that 


tees ele ae iighes el Slain, 


where L denotes a straight line which starts from x, and ends in x,. Application 
of this inequality to J gives 


JSC\|x,— | fla) dy. 
K(R;) 


Taking x, and x, so near to each other that /<0/2, we are led, noting (9.7) and 
the above estimates for / and /’, to the result that for any given 6>0 we have 


(9.8) |p (%, h3 2%) — P(x, k:-%,)| <0 (x€E, REM), 


if |%,—*%,| 1s sufficiently small. (9.8) proves the required uniform continuity in 
x of p(-, Rk; x). 
Now we shall enter into the expansion problem for (Pf) (x), where we first 
assume that /(x)€ Co (E). As the first step we prove PARSEVAL’S equality. 
Since g(x, k; x) is the conjugate Fourier transform of H(x, y; x), it follows 
from PARSEVAL’S equality that for Imzx>0 and Imx?=-0 


WH x: 2) Hey x) dz= I ea kix) erm) dk 
E 


M 


eer. ol : { DIG Bed tau Seg 
9) aaecea a) | Lape ~ Tapa” 


M 
Khi{ 6k) bys had: 
Let us multiply both sides of (9.9) by /(x) and f(y) and integrate over E with 
respect to x and y. Then the left-hand side gives, on multiplying by (x? — x’), 
(9.10) (x? — 2) (Raf, Raf) = 6? — 1) (RaRaf, f) = ((Re— Ralf), 


where we have made use of the first resolvent equation®® R,—R,=(A—J') RK, R, 
and freely interchanged the order of integrations in the integral considered, 
noting that it is absolutely convergent: 


aa | (-. x; NAC, vs ll [FOO |F 00) dy dx <oo. 


29 See HILLE & PHILLIPS [12], p. 126. 
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Here let us introduce the function 
(9.11) DAR td) — (20) EN eg (f(x) €Co°(E)). 


Before proceeding further we note that @(k; \Vu+te) is bounded for ke M, 
w€ [a, BJ (0<a<f) and OSeSéq (€)>0), as is easily seen from Lemma 9.2 
and that ®(k, x) >} (k) as x->|k|, as isseen from (8.1) and the fact that f(x) €Co(E), 
and hence the integral is well defined. 

Let us put x2=y+7e. Since f{(x) Co (E) and the integrand of the right 


side of (9.9) 
2i€ ———, 


Tee ayer 1% Bs Yet te) bly, bs You + te) 
is absolutely integrable in k and bounded in x and y, we can freely interchange 
the order of integration in forming the above-mentioned multiple integral on 
the right-hand side of (9.9). This enables us to obtain in view of (9.10), for e>0, 


(0.12) (Ree = Ranta pea x |P(k; Wut te)Pak. 
M 


Now we can avail ourselves of the fundamental relation 


+ [(Ep+ pd bf) — (Ext Ex) = 54, btn [8 wie Rie) bf) de® 


which, incorporated with (9.12), yields for « and 6 such that O0<a<f 


pai a lan dbf] 


1); ince 
~ ity f fale C= Serie a Sa Agave 


ie, fof, aac OO Ve ee) dy, 


IU a 


where we have used the above-mentioned boundedness of ®(k, x) and Fusini’s 


B 
theorem. Let us consider the inner integral f...du. In view of the equation 
a 


B 
eee = f= i Ne a | ki? 
| var oa — tant SAE 
we have for ¢€ [0, &] (&> 0) 
B 
f...dul|<F(k) =C if |k|SVb+1, 


G . = 
ees ik ake \/B +1, 


%° (a, B] denotes the closed interval «<<. 
$1 See ¢.g. STONE [81], p. 183. 
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where the C’s do not depend on ¢. F(k) is integrable, and hence by the bounded 
convergence theorem we can interchange im and the k-integration. Since P(k; x) 
el0 


tends to / (R) for x>|k| as has been seen above, we obtain for {(x) €C$°(E) 


(9.13) 2 (Ep + Eso) in) ee 3) plas i. Pak; 


ya<|k|<yB 
here we have made use of the well-known relation 22 


B j 
0 NG ieee ie) elm 8 <a] 


fore peaches Wipe 
in | @aarrale ie (a 0M ea@< p, 


in which f(u, €) is a continuous function of w and ¢ for wE[a, 8] and e€[0, é| 
(> 0). Letting «> in (9.13), we see that ((E,; —E,_ 9) f, f) =0 for any f(x) ECy°(E), 
which is dense in L,(£) in the L,-norm. This shows that E,=E,_, and f is not 
an eigenvalue of H, and since f is positive and otherwise arbitrary, we conclude 
that no positive eigenvalues exist. This proves the first half of Theorem 6. 
Hence we have from (9.13) 

0.14) (E,—EJ- N= f \il@Pae 

ya<|kl<yp 
for f(x) CCO(E). Letting «+0 and B->oo, we arrive at PARSEVAL’S equality 


(9.15) PAP =J f(k) 2 ak 


for f(x) CCo°(E). We have so far assumed that /(x) €Co°(E). But the extension 
to the general case f(x) €L,(£) can be made in a standard manner*’. Namely, 
we can show that i(k) exists and lies in L,(M) if we take the limit in the mean, 
as specified in Theorem 5, in forming /(k). Thus (9.14) and (9.15) hold true for 
f(x) €L,(E). 

Now let us step forward to the second problem, 7.e., the expansion of (Ef) (x). 
According to Theorem 4 the subspace E,L,(E) is spanned by the eigenvectors 
fp, (x)} belonging to the eigenvalues w,,<0, or {@,(x)} form a complete ortho- 
normal system in £)L,(E£). Consequently we have 


r =2 h i 


In view of (9.15) and (9.16) we are led to the required relation 


(9.16) |Eof 


4) IIe =[2o/IP + IPAPR= Jf lfmpae + > hl 


Derivation of (8.5) from (8.4) is almost obvious*. 
Now let us prove (8.3). It is obvious from the completeness in Ey L,(E) of 


{p,,(x)} that 
(9.17) (Eof) (x) =1.i.m. x fn Pn (X)- 


32 See e.g. TITCHMARSH [33], p. 31. 
33 See e.g. TITCHMARSH [34], pp. 55— 56. 
34 See any books on Fourier integrals, e.g. TircHMARSH [33], Chapter III. 
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It remains to show that 

(9.18) Cee ie eye — (2%)~*Li.m. f w(x, kh) f(k) dk. 
M 


We can obtain from (9.15), if g(x) €Co (E), 


((Eg—En) fe)= f F(R) BQ) ak 


ya<|kl<yp 
a sf i(k) dk fv x, k) g(x) dx 
27)? 
ya<|kl<yp E 
=(2n)-?f|  f ple A) fle) dk) g(a) ax, 
E ya<|kI<yB 


since p(x, k) is bounded and /(R) is integrable for |x <|k|< V6 and x€ E. Co (E) 
being dense in L,(£) in the L,-norm, we have 
(9.19) (Ep—E,) f)()=2a)4 f p(w kf) ak 
ya<|kl<yB 

If we let a0 and Boo, then the left side of (9.19) converges strongly to 
(Pf) (x), and hence the limit in the mean of the right also exists and is equal 
to (Pf)(x). This is nothing but the desired relation (9.18). We have proved 
i) and ii) of Theorem 5. 

Now we show the diagonal representation (8.6) of H without complete 
characterization of the domain. For this purpose it suffices to show 


(9.20) (H rm sel i(k TH dk + Said 


for {(x)€ D(H)=D/(H,) and g(x)€L,(E). We see from (9.14) that 


(Ef, @) =D fbn (u <0) 
(9.21) ae 


=Lhe t+ SHEAR  (uZ0):; 
hence (9.20) holds if the relation 


(Hf, g) i I wa(E, fF, g) 
is taken into account. 


Finally let us prove the last half of Theorem 6. Taking g=f in (9.21), we 
can see that 0,(u)=(E,f, f) is absolutely continuous in « for ~>0. This implies 
nothing but the absolute continuity of the spectrum of H on the positive real 
line. 


III. S-matrices 


§ 10. Wave operators W,. As has been mentioned in the Introduction, we find 
it convenient to introduce as wave operators W, not only for the mathematical 
treatment of the S-matrix S, but also for the ay of Theorem 5, iii). 
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It is known*® from the time-dependent theory of the S-matrices that, under 
the assumption that g(x) €L,(E), the wave operators W,, exist: 


(10.1) We = s-lim Ui(ey-28 3 U(t) == plea Aa iho 


t—> +00 


WV are isometric operators satisfying 
(10.2) PW,=W., 


which implies that the ranges of W, are subsets of PL,(E). The S-matrix is 
defined by 


(10.3) Sa AW: 
and it is generally believed to be unitary, which is, however, not self-evident. 
S is unitary if and only if the ranges of W, are identical. We shall prove that 
S is unitary by showing that the ranges of W, are identical and coincide with 
IZING ae 

To this end we introduce the operators U, by 


(10.4) (Uf). (4) = @a)7# lism. (x, 2) i, (kh) dk. 
M 


(10.5) Opi UT, 
where /,(k) is the Fourier transform of f(x) € L,(E) (see (9.5)), and prove 

Theorem 7. U,, formally defined by (10.4) and (10.5), are tsometric operators 
with domains L,(E) and ranges PL,(E), and we have Us re S defined by 
(10.3) 7s unitary. a 

We shall prove the above theorem in the next section and then complete 
the proof of Theorem 5. Here we note that the statement ii1) of Theorem 5 is 
equivalent to the following : Let Z be an isometric operator from PL,(£) to 
L,(M) defined by Z/(k)=/(R) ({€ PL,(E)). Then Z is unitary. (The isometry 
of Z is already shown in (8.4).) 

§ 11. Proof of Theorem 7. We first show that U., are well defined operators. 
To this end we need only to show that this is the case with U_, since U, is related 
to U_ through (10.5). 

Let us consider the integral 


(11.1) L(g) =S fo BM Ak= ay *f fo | Sole) g(a) da] ak, 


where g(x) €Co°(E), iy (R) €Co°(M) and the carrier of i(k) does not contain the 
origin of M. The totality of such functions of k will be denoted by Cq°(M)’. 
In (11.1) the x- and k-integrations are, actually, extended over compact domains 
and g(x, k) is bounded in k in a compact domain from which the origin is ex- 
cluded (see Theorem 3). This fact allows us to interchange the order of inte- 


gration, yielding : friar 
(11.2) L(g) = (27)-¥ f | f 9 (%,&) folt) 48] g(a) dx. 


35 See Cook [7], Karo [1/7] and Kuropa [20]. 
36 “‘s_lim’’ denotes strong limit. 
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On the other hand it follows from (11.1) that 
Lis (Sloe ar flame ae) [irl Pell SMAlliisll 


Thus L,(g) is a densely defined, bounded linear functional of geCo(£). But 
such a functional can be extended to the whole L,(£), and we shall also denote 
the extension by L,(g). Then by means of Riesz’s theorem*’ there exists a 
unique element /* of L,(£) such that L;(g)=(/*, g), which defines an operator 
U by Uf=f* and ||Uf||<||f||- Now in view of (11.2), again restricting g(x) 
to Co(E), we have 


(11.3) (U/) (x) =22)*f p(x #) folh) ak ae. (fo(k) €Co(M)). 


A standard argument shows that (11.3) can be extended to every f€ L,(E) by 
writing l.i.m. {...dk instead of {...dk. Comparison of (11.3) with (10.4) gives 
M 


the result that U is nothing but the required operator U_. Thus we have proved 
that U_ is an everywhere defined, bounded operator (D(U_)=L,(E)) and its 
norm ||U_ || does not exceed 1. 

Next let us determine the adjoint U* of U_. This can be done most con- 
veniently by using the relation (U*/, g)=(f, U_g); the result is easily seen to be 


ae 


(11.4) (OA) (ai) a= (Soe * laa, ais fel ea eee 


In the derivation of (11.4) we have used reasoning quite analogous to that for 
getting (11.3). 
U_ has the property that 
(11.5) UT SeaePt 
which follows immediately from the defining equations (10.4), (11.4) of U_, U* 
and (8.3). Also we have 
(11.6) UP TCH. OS 
This can be seen from the following equations: For f(x) € D(H) =D(H,)) 


(ULE ix)== (2a)-tLim. fe? lay) tae 
M 


= (2m)~! Lim. fei*-*|R 27 (A) dk: 


M 


(H, U* f) (x) = (2m)-! Lim. f ef" |hl*(U* f) fo (BR) dk 
M 
= (22)-! Lim. f e** | ef (A) dk, 
M 


where we have made use of the diagonal representation of H (8.6) and the fact 
that Ho is a special case of H with V=0 and has no eigenvalues. 

Now we can complete the proof of Theorem 7. Application of U* from the 
left to both sides of (10.1) gives 


(14.7) U2W_= lim U*U() — (lim=s-lim), 


7 See e.g. Riesz & Nacy [29] and Yosipa [39]. 
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since U* is a bounded operator and lim U(t) exists. On the other hand, we see, 


t—> —oo 


by differentiating and again integrating, that 


t 
(U* Ul) fg) —(U* fg) =i f (US bt Vet Mf, g) dt 
: 0 
t 
=a ee iis View oe] '2) at 
0 
where /(x) €D(H,)=D(H) and g(x) €L,(£) and we have made use of the fact 
that U* e!# = e¢''#sU* | which is a consequence of (11.6). In view of (11.7) and 
(11.8), we obtain 


t 
(11.9) (CEs 72) = (CS 7p, 2) rc lismeslile? Cea eas 
> 9 


(11.8) 


The integrand of (11.9) is calculated as follows: 


(e'#t0U* Ver sf 2) =?) ah (UF Veet), (R) By (k) dh 
p A eit RP (Veit He iyo (Rk) &(k) dk 


=jv e (MoKIRN) #)” (k) Bo (k) dk 


= (2m)-*f | fg (%H) gta) (eB (x) dx] Bo(h) ak, 


(11.10) 


where l.i.m. is not needed before {, because g(x, k) is bounded in x, g(x) €L,(E) 
padi (en oF) A) (4) CLs (2). ¢ 

Returning to (11.9), the limit for t+ — co of the f-integral can be replaced 
by the Abelian limit, 7.e. 


lim f (et U* Ve-ttHe fevdt = lim f (es UEVe of, 2) at 
§ 2105 § 


t+ —O 9 


since the existence of the ordinary limit is known. Then (11.9) and (11.10) give 


ie ig y a eet dtx 
0 
x f| fee ) (IN f(x) dx] Fo(B) de 
M E 


If we assume here that 2 (k) CCq°(M)’, then we can interchange the ¢ and R- 
integrations, for g(x, k) is bounded for x€E and ke D=carrier of go(k) (see 
Theorem 3) while 

{lata e tt j(2)] ds [lll 


(DEW fe) = (OE ey te tim | aati ff (ee oy, Fe yas, 
M 


where we have put 


F(x, k) = (x, k) 9(x)€La(E) (for each fixed kM). 
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Hence 


(U*W_f,e)=(U* 2) — an lim [ 2 () dk x 
ea aaa 


x Lam | (fe ete Aa Tope ae one 


t— —c 
= (U2 f.e)+—* lim 3 Bo (®) (Ro, pie fs F(-, B)) ah 


(22)! e|0 
=i CE Hg) si lim | (f, Ro,jerie (> ®) G(-)) Bo (A) aR. 


2n 3 e|0 A 
M 


(11.11) 


Here we further assume that f(x) CCo°(E) CD(H,). Now the function 


(Ro, init -ie P (+, R) q(-)) (*)= a fri p(y, k) dy 


converges to 


; 2ilallx—yl 
ta | emp 1 9-H d= — (Tne) 


as ¢/0 uniformly with respect to x€E and k contained in the compact carrier 
of g,(k) by virtue of Theorem 3 and Lemma 4.6. Thus the convergence of the 
integrand of (41.11) for e€) 0 is uniform, and we can interchange i and hes ak 
so as to obtain 


{ oe RI |z—91 


WWI.) = WF d+ 1, [aan [Fax t [Fa oly, aay 
M E E 


laa st 


Banc (k) aah hs findx (by (62) 


(22)? 


Peat et kl |*—9| b)d lq 
y (x, k) 4 al eran ae a x 
E 


(11.12) 


= | o(h) B®) dk = (fe). 


Co’ (M)’ is dense in L,(M), and hence the totality of such g(x) that g, (k) €C9°(M)’ 
is dense in Ly(E). Cy°(E) is of course dense in L,(E£). We can, therefore, obtain 
from (11.12) the following relation: 


(11.13) OS Ware I, 
If we operate with U_ on both sides of (11.13) and take account of (10.2) and 
(11.5), we get 

(11.14) i= U (UP We = PWo=W. 


and, on substituting this relation into (11.13), 


(44545) U* U_=W*W_=I. 
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The relations (11.5) and (11.15) show that the operator U_ is an isometric operator 
with domain D(U_)=L,(E) and range R(U_)= PL,(E)*8. Quite a similar result 
follows for U, and W,.. Thus we can see that S=W* W_=U+* U_and Sis unitary. 
Now the proof of Theorem 7 is complete. 

In connection with eigenfunction expansion it remains to prove that the 
isontetric operator Z is unitary, to characterize the domain of H, and to establish 
the unitary equivalence between the part of H in PL,(E) and H,. To these 
ends we shall define the operators Y and Z’ by 


(Y f) (2) = fo (A) (from L,(E) onto L,(M)); 
(277) (4) = (2or) at hima fry, 'k) fk) ak (from L,(M) into LX(E)). 


Then we have U_=Z'Y and U*= YZ, and we have 


(11.16) LOL WBZ VIZ) AT 
On the other hand (8.3), (8.4) and (8.5) imply that 
(11.17) B= ssiZe 


(11.16) and (11.17) show that Z is a unitary operator from PL,(E) onto L,(M), 
and that Z’=Z*, the adjoint operator of Z (Z* is a unitary operator from L, (M) 
onto PL,(E)). 


Now it follows from (8.6) that 
(Z PHA) (k) =|AP/(), 


and it is known that H, represents the operation of multiplying by |&|? in L,(M). 
Thus we have proved iv) of Theorem 5. 

Finally let us prove the first assertion of Theorem 5, v). Since the necessity 
is obvious from the diagonal representation (8.6) of H, it remains to show the 
sufficiency. We see from the unitarity of Z and (8.6) that if F(k) CL,(M) and 
|k|?F(k) CL(M) then Z*F (x) € PD(H) (D(H). Thus we need only to prove that 
for any sequence F, such that >) |F,|?< oo and )) ui|F,|?< 00, li.m. D) F, 9, € 
E,D(H) D(H). But this follows easily from (8.6) and the fact that the corre- 
spondence between F,L,(E£) and /,*°, which assigns every fC EyL,(E) to f, Cle, 
is one-to-one and isometric, where /,; denotes the Hilbert space of all sequences 


{F,} such that >) |, |?< oo. 


IV. Concluding remarks 


§ 12. Higher-dimensional and two-dimensional cases. Our theory, which has 
so far been developed for the 3-dimensional Schroedinger operator, cannot be 
extended directly to the case of higher dimension. One of the main difficulties 
is that Theorems 1 and 2, which serve as a bridge between the resolvent R, and 
the kernel G(x, y; 4), fail to hold, because the kernel G(x, y; A), even though 
it exists, is not of the Carleman type in the case of more than three dimensions. 


38 R (A) means the range of A. 
39 J’ denotes the identity operator in L,(M). 
40 If the dimension m of E,L,(E) is finite then /, is the totality of finite sequences 


of length m. 
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A remedy for it has been suggested by GARDING [9]; it is to consider the operator | 
(H'—A)4, where ¢ is a certain positive integer, instead of R,=(H—A)7. The 
operator (H‘—4)+ is, in fact, known to have a kernel of the Carleman type. 
Of course, this will require many an alteration in the details of our theory. 

On the contrary, in the 2-dimensional case our theory suffers no essential 
modification, for then the resolvent kernels can be shown to be of the Carleman 
type. But it should be remarked that results for the 2-dimensional case do not 
follow as special cases from the 3-dimensional treatment, because the potential 
function g(x), which diminishes at infinity in the 2-dimensional Euclidian space 
E2, does not always diminishes at infinity in E=£%, which is an extension of 
E. The conditions imposed upon q(x) are to be replaced by 


(A’) q(x) ts a real-valued function which 1s locally Hoelder continuous except 
at a finite number of singularities. Furthermore, q(x) ts square integrable 
(q(x) €L,(E?)) and behaves like O(\x|~*") (h>O0) at infinity, t.e. there exist 
positive numbers h, Cy and Ry such that 


lo(ls Cole [—te* for xh Ro. 


Not all the details of the theory need be modified if we replace the function 
(40|x — y|)te'*l*—9l, which is the resolvent kernel of Ry,2, by“ 47H) (x|x—y)), 
where H() denotes the Hankel function of first kind of order 0. 

In concluding we remark that if we had assumed in addition to (A) that 
q(x) €L,(E) or more explicitly g(x) =O(|x|~*~") at infinity (in the 3-dimensional 
case), our arguments could have been much simplified. 
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Bounds in the Cauchy Problem for the Laplace Equation 


ee ANE 


Communicated by C. LOEWNER 


1. Introduction 

It is well known that the Cauchy problem for the Laplace equation is not 
well posed in the sense of HADAMARD%, 7.e., the solution does not depend in a 
continuous way on the Cauchy data. Because of this fact and perhaps because 
of the feeling that such problems were unnatural, mathematicians have until 
recently neglected the study of the Cauchy problem for elliptic equations. How- 
ever, during the past few years considerable attention has been given to rather 
general problems of this type, especially to the study of uniqueness of solution 
[1—14, 16, 18, 19|**. One of the reasons for this recent interest is the fact that 
although such problems are not well posed, they nevertheless actually correspond 
to concrete physical situations. The physical problem is usually one in which 
a portion of the boundary is inaccessible for measurement of boundary data. 

It has been demonstrated by LAURENTIEV [2] that if one prescribes in addi- 
tion to the Cauchy data a supplementary condition, namely that the harmonic 
function be uniformly bounded by some constant &/ in its region of definition &, 
then the Cauchy problem for the Laplace equation has a stable solution. The 
problem is then well posed. It will be noted in Section 2 that various other 
conditions are sufficient to insure stability. 

In this paper we present a method for obtaining pointwise bounds for the 
solution to the Cauchy problem for the Laplace equation in u-dimensions. If 
the Cauchy data is known exactly the error term may be made small by the 
Rayleigh-Ritz technique. However, it usually happens that the Cauchy data 
is obtained by measurement and is not known precisely. On the other hand, 
an upper bound for the absolute value of the error in measurement of the data 
can usually be determined. We show that knowledge of this error and of an 
upper bound for the absolute value of the harmonic function in Y leads to a priort 
bounds for the value of the harmonic function at any point in @. 

The method proposed here is in a certain sense a generalization of the work 
of LAURENTIEV [12]. In fact some of the results of LAURENTIEV may be used 
to give pointwise estimates near the portion of the boundary on which the data 
is given, but only if part of the data is homogeneous and if @ satisfies certain 


* See, for instance, J. HADAMARD, Le Probleme de Cauchy, Paris 1932, or refer- 
ence [17]. ; 

*k See also P. R. GARABEDIAN & LIEBERSTEIN, H. M., On the numerical calculation 

of detached bow shock waves in hypersonic flow. J. Aero. Sci. 25, 109—118 (1958). 
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geometrical requirements. Even if these conditions are met, LAURENTIEV’S 
method still fails to yield pointwise estimates at every point in %. 

The results of this paper may be extended to more general types of equations 
or systems of equations, e.g., polyharmonic equations, equations with variable 
coefficients, etc. We restrict consideration here to the Laplace equation, because 
it is the most important case physically, and because certain simplifications can 
be made which are not possible in the treatment of the more general problem. 
The more general problems will be dealt with in a subsequent paper. 


2. Notation and Formulation of Problem 


We consider an n-dimensional domain @ bounded by a closed surface @. 
On a portion » of @ Cauchy data is prescribed. We denote the remainder of the 
surface (6 — 2) by 2, and take the origin of a spherical coordinate system at 
a point P,. This point is chosen in such a way that there exists a sphere K(z) 
of radius y and centre at P, which has the property that the intersection of @ 
with A(rv) contains no points of ». (If two or more disjoint regions are formed 
by the intersection of K(rv) and J, it is sufficient that the boundary of any one 
of these portions contains no points of >.) We denote by &, the intersection of 
K(r) and @ (the portion whose boundary contains no points of » if two or more 
disjoint regions are formed), and require also that FR not be a point of Z,. We 
denote by R the radius of the largest sphere satisfying the conditions prescribed 
above and by 7) the radius of the smallest such sphere. We can always locate 
P, in such a way that K(7,) is tangent to Z at some point of 2. (One need only 
use a portion of 2 if necessary or convenient.) We assume this to be the case 
and prescribe further that A be so chosen that 7)>0. (The point P, should 
not lie on 2.) In most instances P, will not be a point of Z but will lie outside 
of @ near the portion 2. It is obvious that we have considerable freedom in 
the choice of &. The optimum choice will depend on the geometry of 2 and 
upon the Cauchy data. Finally, we denote by w(r7) the portion of the surface 
of K(v) which les interior to @. 

Let w be a function uniformly bounded in @ satisfying 


(2.4) Aw=0, A= 
i=1 °? 


therein and satisfying on the portion 2 of the boundary the conditions 


f (w—f2do<a, 
(2.2) z 

f (w,;—F) (w,;— F) do<az 
where do denotes an element of surface area on 2, the comma denotes partial 
differentiation, and the repeated index denotes summation over 7 from 1 to ». 
The quantities f and F; are the measured values of w and its gradient on J, 
and «, %, are known constants. 

In order to obtain pointwise estimates for w at any point in J we first derive 

a suitable a priori estimate for an arbitrary function u which satisfies (2.4) in, 
Then by choosing uv in an appropriate way we obtain the desired bounds for w. 
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We introduce the following notation: 


fi eee do, 
(2.3) z 
= fu ,u ,do 
and 2 
(2.4) A,(u) = fu ju ,dx 
Dr 


where dx denotes an element of volume of Y,. We assume further that 
(2.5) |u| su 

in Y. We shall derive an a priori estimate for u of the following type: 
(2.6) HA Pani eeeseeiiogl 2 te On Ke Ai el? 


where Ay, A,, and K, are determined constants and 6 is a known quantity 
B= 0 <1. 

Before deriving (2.6) we point out how such an a priori estimate may be 
used to obtain pointwise estimates for the function w satisfying (2.1) and (2.2). Let 


(2.7) uU=WwW—@ 


in (2.6), where is an arbitrary harmonic function in @*. Then 


|w(P) —9 (PP SK MOOK, f [ww —/) —(p— P40 + 


+, J [(@,;—F) — (p,—F)) (@,;— F) — (9, :— B)] do 
S 2K, M*-1K, lay + f (p —f)?4o| + 


+ Koa + J (ps —F) (9,.—F) dol}. 


Thus if «, and «, are small and @ can be chosen to approximate the data of w, 
then it is possible to obtain close upper and lower bounds for w(P). The usual 
procedure for approximating the data is to choose 


(2.8) 


N 
(2.9) y ay bP, 


(where the g,; are N independent harmonic functions) and use the Rayleigh-Ritz 
technique to minimize the right-hand side of (2.8). Thus, once we have obtained 
(2.6), pointwise estimates for w(P) follow immediately. We turn now to the 
derivation of (2.6). 

3. Preliminary Estimates 


For >2 we define for any function wu a functional in the following way: 


(3.1) Vir, u) = log F(r,u) + Cr—?™?) 
where 
(3.2) F(r,u) =f o'-"A,(u) do + Kya+ Kee 


Yo 


* It is clear that if w is uniformly bounded by some constant M, in and pisa 
prescribed function, the constant MM in (2.5), with w defined by (2.7), can be computed 
easily. 


38 L. E. PAYNE: 


and ce (V2+1) 7 ier ras 3 
ACEP : 
r nH 
(3-3) l= Fao 
na 
a= 32) 


We show that with the choice of constants given by (3.3), V(7, u) is a convex 


function of ~~). To prove this we form 
» &V(r, u) 
[F (7, u) |? (07— MD) 


Soe aE u) uu ds — (4, (u)}*} + 2C[F(r, w)]2. 


o(7) 


(3.4) 


Let &’ denote the portion of 2 lying inside A(r); then since w is harmonic 


F(r,u) = ea ge Fagus) thet 2») — 0-9) 4» SH dg + 


ve) 1 y= CB (n—1) 90 42 
EE apa ical) 2 — a= Cs BI get ae 
et fwdstsfe An Yao 
@(r) Day! 
where 6/0” denotes the outward normal derivative on 2’, and o measures the 
distance from the origin P. Clearly 


ve n—D 7 8) 
[ee rerun 


(3.6) eke 
AUGe ent ee 


to. (n—1) 


(n—2) 


—(n—3) 
Sates 
2(n—2) |? 


frwass IK, ! 


a(r) 


é4 [Ket 5 


Also by GREEN’s identity and the Schwarz inequality we have 


+ > Je furas +R [uimaas [Rte + R e.|+ & &. 
a(r) (7) 
We need also the inequality 
du\2 = 
(3.8) [ [leradsuP — (32) |@s2—J2e, 
x) 


where grad, denotes the tangential projection of the gradient. Inequality (3.8) 
follows directly from the RELLIcH identity [20], which for the region Y, reduces to 


al 


(3.9) w(r) 


Jgrad, w)?— (FE) |ds=2 f ograd,g- gradu S do — 


Sopigeees || grad, w|? 2 (se) | do+(n —2) A,(w). 
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This identity has been used previously for obtaining bounds in Dirichlet and 
Neumann problems [1/5]. We now insert (3.6), (3.7) and (3.8) into (3.4), ie., 


(PG, «Pp oSmagn 2a ei) ae {5 f weds f [leraa ul? all ds + 


y a(r) w(r) 
le —(n—1) * —(n—3) pie 
+ (a — ene Ky Ae 1 fu su sds — 
(3.10) oe 
—S)Ra futastr fu U as] [Rte+ Re] —e ee 
2 Boat A 2. 6 7 
(7) o(r) 


2 

y—"M—I) e 2 igs ED Va 

F2¢| sf ds + (Ky oman | (K, ss) 
o(r) 

If we replace the second integral on the right by — |/2 ¢, and make use of (3.3), 

we find that the right-hand side is non-negative. It follows then that V(r, w) 


is a convex function of r~"~?), 


It is clear that the values for C, K,, and Kk, given in (3.3) will not in general 
be the best possible values. The best choice will depend on the geometry of 2, 
the values of ¢, and ¢,, the number of dimensions and the location of Pj. For 
instance, instead of the left-hand portion of the inequality (3.6) one could have 
written 
3.11) F(z, u) — 


y—"n—-) 


eos ees | a 
[fu ae Ky 2(n—2) fae | 2(n—2) ao 
a(r) 
where « is an arbitrary positive constant. Similarly (3.7) could have been re- 
placed by 


[4(P< f weds [ (S*)ds+ 
(3.12) (1) o(r) 
+4|(BR)2 


o( 


weds +BRfuiujds [(y R2) Hy Re] + & 5 
r) a(r) 


( 


where f and y are arbitrary positive constants. Thus, C, A, and kK, could have 
been made to depend on the three arbitrary parameters a, B, and y. 


In two dimensions we define V (7, ) as 


(3.13) Vir, u) = log fet4,(uw)de+ Ki ene | + C’(log 7)? 


where C’, Kj, and Ky are given by 
(3.14) ee 


It is easily verified that V(r, w) given by (3.13) is then a convex function of log7. 
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We assume again that 7>2 (the same procedure applied to (3.13) will yield 
similar results for »=2). From the convexity of V(r, u) for 797K it follows 
that 

exp[iC x7") Rye) 


(3.45) <{exp [C Roe 2) F(R, u yyi- °fexp[C 2(n—2) Fy. u)}? 
where 
(3.16) b= (22) et yA 4), 


But from (3.1), (3.3) and (3.6) we see that 


UES Uke ye it Gee | [im—2) R?-+275) ss. Me 
: Diten Bee 2 fap : 


647) FR,» 3|S ae + a 


where S(R) is the area of w(R). Note that when upper bounds for ¢ and é, 
are known then M, is completely determined. Since these upper bounds will 
be computed%, let us write 


(3.18) | 6 C,M?,  <C,M 


where C, and C, are known constants. Then instead of Mf we write 


(3.19) M? < B, M? 
where 

=P USi(R} (n+2) C, [(m— 2)? R? +272] 
B20) P= spt + 2(n—2) G1 a2) 


From (3.2) it is obvious that 


(3.21) F(7,, u) = Ky 4 -+ Keo és. 
Thus (3.15) may be rewritten as 

(3.22) Fv, 0) = BM? (Ky ee Ke) 
where 


(3.23) B= Bis exp{C [ie 0) Roe 2) eg gee aby 
Using (3.6), we then have 
(2224 very = ey RACES ACh Me Sy fer eM gy OB Mee Se ere ene 
Note that a somewhat simpler expression arises if we know instead of (2.5) 
a bound for A(w), 7.e., 
(3.25) A(u) = fu ,widxs Ll. 
g 


This condition is obviously sufficient to insure stability. We have then in- 
stead of (3.17) 


, < L? wits a Ps c 
(3.20) F(R, u)< ie (gee ReaD Rie eee raal e 


* See the last paragraph of Section 4. 
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Thus, one could write instead of (3.24) the expression 


P) tte Ura s A ott? 7") e, 
DIE o(7) , . = , r Pa 
Cy Les) (A, e+ Ky &9)° +C3(Kya+ Ky 2) 


where C; and Cj are easily determined. 
It is clear also that stability and pointwise bounds will follow if we use 
a bound of the following type: 
(3.28) J |u| | grad u| dv< Lj 
g 
instead of (2.5) or (3.25). 


4. Pointwise Estimates 
We derive first an a priori estimate for the value of u at a point P in @,. 
Let G(P, Q) be the harmonic GREEN’s function for the sphere K(7) with singu- 
larity at P. Then by GREEN’s formula we have 


(4.1) = sid Oe ae ot do — fu ride: 


If the point P is not a stant point of Z,, equation (4.1) may be differentiated 
any number of times with respect to the variables of the point P, the order of 
differentiation and integration being interchangeable on the right-hand side. 
Using SCHWARz’s inequality we obtain 


eR ee fats hy Oe eC Ray. et) 65] x 


w(7) 
ye f (ey ds+ rn ch & i do +r) R- aac 


o(r) » 3, 


(4.2) 
x 


If we replace the first bracket term on the right using (3.24), we find 


| «( oe) ee OBONT 2 ois 2K 1&1 + Ky &9)°X 


ro f(s 2) ds +r ak 
én 
(7) 


Inequality (4.3) thus gives an a priori estimate for the value of w at a point P 
in Y, which is of the desired type (2.6). It depends on the uniform bound M 
and on the Cauchy data in such a way that if the data tends to zero the value 
of u(P) tends to zero. 

We denote by &; the variables of the point P. Then by differentiation of (4.1) 
we obtain 


(4.3) 
x 


ii (Sy dot+trn*R? it Grado] 


yy 


7 : ou 
(4.4) uw (P)=— fuse Gidst+ [ 6, do— [« 4, Gide 
w(7) Da De 
where the comma in this case denotes differentiation with respect to the variable 
E,. If P is not a boundary point of Y,, we again apply the Schwarz inequality 


and obtain a bound for any derivative of w at the point P. Since we shall need 
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it later, we write down the estimate for the gradient of w at P. 
[u (P)u (P)]S2B M*0-9) (Kye, Ky @5)° x 


(4.5) , ; 
call cairn ce con 


alr) x 


Woke G+ R*6,6,] da) 
On 2 On i= eee x 


Suppose now that we wish to obtain bounds for the value of w at a point 
P in J—@Q,. We proceed in the following way: Let Y* be a subdomain of GQ 
whose boundary @* coincides with @ except for a small band 4+y2e2n—y 
where y is a small positive constant and 7>1>7. If w*(n) is defined relative 
to €* as w(%) was relative to @, we prescribe that Y* be so constructed that 
w*(r,) and @ have no points in common. In practice the sphere A(7,) will first 
be decided upon, and then @* will be determined. We now integrate (4.3) and 
(4.5) over the surface w*(7%), 7.¢., 


m= fds <2B M?-9(K, e, + Ky &) x 
w*(7;) 


ela fe ts arn a eee eee 


w* (71) o( w*(r,) 2 


(4.6) 


and 
aie fu,uj,dsS2BM*"—9(K, e+ Ky &)°x 
* (11) 


(4.7) x (r 


4 a 
Eee as* 
on on a 


ae fi ee RAGL Ee Jods}. 


on on 


Note that instead of (4.6) we could use (3.24) for y=” and the fact that 


(4.8) | dst foPds. 
wo*(r) wr) 

The object of introducing the subdomain Y* was merely to avoid integrating 

all the way to the boundary of Y, in (4.6) and (4.7). For then the integrals on 
the right would not have converged. Actually we do not use Y* or €* but only 
*(r,). Consequently another way of stating what we have done is to say that 
we formed (4.6) and (4.7) by integrating (4.3) and (4.5) over a portion w*(7%) 
of the surface of the sphere K(%) (where 7>7,>79) with the stipulation that 
every point of w*(7,) should be in Y, at a distance greater than some number 
f from @. 

It is obvious from (4.6) and (4.7) that we may regard the surface w*(7,) as 
a Cauchy surface, that is, one on which inequalities for the Cauchy data are 
known. Our object is to obtain an estimate for the value of « at a point P in 

— 9* with boundary @*. Thus we consider the following problem: 


(4.9) Au=0 in O*— g* 


with upper bounds for 7, and 72 known on w*(7,). We proceed as before, letting 
«*(7) be our new surface 2* and letting @* —2* be the new »*. We choose a new 
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origin Ry* to have the same properties with respect to the new »* that Py possessed 
for the original ». We define new quantities R* and r¢ in the obvious way, 
and, proceeding as in Sections 3 and 4, we are able to obtain bounds for the value 
of wv and its derivatives at any point P* in GF. (By GF we mean the intersection 
of J* — G* and K(r*) with center at R*.) These bounds will be given by (4.3) 
and (4.5) with the coordinates changed to starred coordinates, the constants 
B, Kk, and k, replaced by constants B*, Ki and K, and e, and ¢e, replaced by 
m, and 7,. For example 


| «(P*)|? SSD IME oa) Kaediyate Kowa). xX 
n— oGen? G* 4 
for f (Case + var [|Z P+ ce (cr9 co 
w*(r*) z*’ 


with 7, and 7, defined by (4.6) and (4.7). It is obvious that (4.10) is of the form 


(4.10) 


(4.11) |\w(P*)\\?= K*MPA0— 999 (Kye, + Ky &9)?? 


where the K* is a determined quantity. 

If we wish to obtain bounds for the value of « at a point which is neither 
in J, nor in YF, we again define a new Cauchy surface, and in fact continue to 
define new such Cauchy surfaces until we arrive at the point at which bounds 
are desired. It is obvious that there is considerable arbitrariness in this procedure, 
and the optimum estimates will depend on the geometry of 2’, the values of « 
and é, etc. 

One can also obtain in an obvious way pointwise estimates for any derivatives 
of wu in ZY. It is likewise possible to obtain pointwise estimates for w and its 
derivatives if the region is unbounded. No conditions need be imposed on the 
behavior of w at infinity. 

We note that in obtaining pointwise estimates for the solution of (2.1) and 
(2.2) we subsequently set w= w — , choose the @ as indicated in (2.9) and minimize 
with respect to the coefficients. Regardless of the number of intermediate Cauchy 
surfaces it may have been necessary to define, the constants 0; in (2.9) must 
always be chosen to minimize the expression 


(412) Kyla +S (p —A8do] + Kalta + f (9,1 — Fi) (p,.— B) do 


occurring in (2.8). Thus one starts by choosing a large enough number NV of 
harmonic functions g; to make the expression (4.12) sufficiently small. The 
constants C, and C, of (3.18) are then easily computed. The determination of 
the constant K, in (2.6) is then straightforward. 


5. Practical Application of Results 


It is apparent from the procedure described in the preceding section that 
the precision with which we can approximate the solution of (2.1) and (2.2) at 
a point P with a given set of functions g; depends very much on the number 
of Cauchy surfaces which we must define in order to arrive at P. Thus if 2 
is sufficiently small or if the domain J is sufficiently thin or irregular (for example, 
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a region in the shape of a dumbbell with Cauchy data given on a portion of one 
of the knobs), it may not be possible to get practical bounds for w at a point 
sufficiently far removed from 2. In many physical problems, however, the point 
P may be used to determine X. If we wish to determine certain steady state 
atmospheric conditions at a point P in the atmosphere, we take measurements 
of Cauchy data over a portion 2 of the surface of the earth which is so large 
that the point P will lie inside the original Y,. A number of interesting problems 
in electrostatics and magnetostatics, as well as gravitational potential problems, 
may be handled in a similar way. It should perhaps be pointed out that the 
gravitational potential v satisfies a Poisson equation, but if we approximate 
with functions p satisfying the same equation, then the quantity wu=v—y will 
be a harmonic function. 

The results of this paper may likewise be used to investigate the field (flow, 
electrostatic, magnetostatic, etc.) which gives rise to a given free surface. 

Also in certain boundary value problems in which the point P lies close to 
the boundary it may often be more convenient to measure Cauchy data quite 
accurately over a portion of the boundary near P and use the results of this 
paper, rather than to consider the standard boundary value problem. This is 
particularly true if the region is extremely large or if portions of the boundary 
are irregular. 

It may happen that isolated singularities (or singular surfaces of dimension 
less than »—2 in general) occur in Y. In theory one can obtain bounds for 
u(P) at any point P outside a small neighborhood of the point at which the 
singularity occurs. In practise one may not be able to obtain close bounds at 
each such point P unless the constants a and a, in (2.2) are extremely small 
and the number N in (2.9) is extremely large. 

In other problems the portion of the boundary X may not be known. One 
may wish to choose it to be a level surface for a larger (possibly infinite) domain 
as is the case in certain design problems. The method developed here may be 
applied to such a problem. 

If it were possible to measure the Cauchy data exactly, and if »' were suffi- 
ciently smooth and we knew uniform bounds for a sufficient number of derivatives 
_ of wu, then the procedure could be simplified considerably. We could actually 
find a bound similar to (3.24) for f|gradw|?ds, and it would not be necessary 

o(r) 
to define Y*. However, it is Fae ae unlikely that this will ever be the case 
in a real physical problem. Also if we could obtain bounds of the type (2.2) 
for higher derivatives of w on 2, and if » were sufficiently smooth, then we 


could again find a bound similar to (3.24) for f |gradu|*ds. It would again 
w(r) 

be unnecessary to define Y*. However, one must expect poor accuracy in the 

measurement of higher derivatives of w on 2, so that the net result would usually 

be a worse bound for w(P) than would be obtained by the method suggested 

in this paper. 


This research was supported in part by the United States Air Force through the 
Air Force Office of Scientific Research of the Air Research and Development Command 
under Contract No. AF 18(600)-573 and in part by the National Science Foundation 
through a postdoctoral fellowship. 


id 


Bounds for the Cauchy Problem 45 


References 


ARONSZAJN, N.: A unique continuation theorem for solution of elliptic partial 
differential equations or inequalities of second order. Math. Pure Appl. 36, 
235— 249 (1957). 

CALDERON, A. P.: Uniqueness in the Cauchy problem for partial differential 
equations. Amer. J. Math. 80, 16—36 (1958). 

CaARLEMAN, T.: Sur un probleme de unicité pour des systémes d’équations aux 
dérivées partielles a deux variables indépendantes. Ark. math. Astr. Fys. 
B 26, 1—9 (1939). 

Corvbes, H. O.: Uber die Bestimmtheit der Lésungen elliptischer Differential- 
gleichungen durch Anfangsvorgaben. Nach. Akad. Wiss. Géttingen, Math.- 
Phys. Ila, No. 11, 239—258 (1956). 

Dovuettis, A.: Uniqueness in Cauchy problems for elliptic systems of equations. 
Comm. Pure Appl. Math. 6, 291 —298 (1953). 

Foras, C., G. Gusst & V. PorNaRu: Sur le probleme de Cauchy pour le type 
elliptique a deux variables. Acad. Rep. Pop. Rom. Bul. Mat. Fiz. 7, 97—103 


(1955). 


| Hartman, P. & A. WiInTtNER: On the local behavior of solutions of non-parabolic 


partial differential equations. III. Approximations by spherical harmonics. 
Amer. J. Math. 77, 453—474 (1955). 

Hernz, E.: Uber die Eindeutigkeit beim Cauchyschen Anfangswertproblem einer 
elliptischen Differentialgleichung zweiter Ordnung. Nach. Akad. Wiss. Gét- 
tingen, Math.-Phys. Ila No. 1, 1—12 (1955). 

HORMANDER, L.: On the uniqueness of the Cauchy problem. Math. Scand. 
(in print). 


} Joun, F.: A note on “improper’’ problems in partial differential equations. 


Comm. Pure Appl. Math. 8, 494—495 (1955). 


|] Lanprs, E. M.: Certain properties of equations of elliptic type. [In Russian. ] 


Dokl. Akad. Nauk SSSR. 107, 640—643 (1956). 

LaurREnTIEV, M. M.: On the Cauchy problem for the Laplace equation. Izvest. 
Akad. Nauk. SSSR., Ser. Mat. 20, 819—842 (1956). 

MULLER, C.: On the behaviour of the solutions of the differential equation 
Au=F (x, u) in the neighborhood of a point. Comm. Pure Appl. Math. 7, 
505—515 (1954). 

NIRENBERG, L.: Uniqueness in Cauchy problems for differential equations with 
constant leading coefficients. Comm. Pure Appl. Math. 10, 89—105 (1957). 

Payne, L. E., & H. F. WEINBERGER: New bounds in harmonic and biharmonic 
problems. J. Math. Phys. 4, 291—307 (1955). 

Preperson, R. N.: On the unique continuation formula for certain second and 
fourth order elliptic equations. Comm. Pure Appl. Math. 11, 67—80 (1958). 

Petrowsky, I. C.: Lectures on Partial Differential Equations. Cambridge: 
University Press 1954. 

Pucci, C.: Sui problemi di Cauchy non ,,ben posti‘’. Rend. Accad. Naz. Lincei 

(8) 18, 473-477 (1955). 

Pucci, C.: Discussione del problema di Cauchy per le equazioni di tipo ellitico. 
Annali Mat. Pura Appl. 46, 131—153 (1958). 

Reriicu, F.: Darstellung der Eigenwerte von Au-+ Au durch ein Randintegral. 
Math. Z. 46, 635—646 (1940). 


King’s College, Newcastle-upon-Tyne 
and 
University of Maryland 


(Received January 15, 1960) 


Uber die Randkomponenten 
ebener harmonischer Vektorfelder 


ERICH MARTENSEN und KURT VON SENGBUSCH 


Vorgelegt von H. GORTLER 


§ 1. Einleitung 


Bei verschiedenen physikalischen Anwendungen der Potentialtheorie spielen 
harmonische Vektorfelder! eine Rolle, die auf und auBerhalb einer bestimmten 
Berandung in der (x, y)-Ebene erklart und im Unendlichen beschrankt sind. Die 
iibliche Methode, solche Vektorfelder aus gewissen Randbedingungen zu berech- 
nen, besteht in der Darstellung des Feldes als (negativer) Gradient eines skalaren 
oder als Rotation eines (bei ebenen Problemen nur eine Komponente besitzenden) 
Vektorpotentials. 


Uber die Einwertigkeit solcher Potentiale ist zunachst noch nichts ausgesagt. 
Man wei8 jedoch, da8 ein harmonisches Feld im Falle verschwindender Zirkulation 
(Linienintegral der Tangentialkomponente) ein einwertiges skalares und im Falle 
verschwindender Ergiebigkeit (Linienintegral der Normalkomponente) ein ein- 
wertiges Vektorpotential besitzt. Bekannte Beispiele hierfiir sind das (zirkula- 
tionsfreie) elektrostatische Feld eines Zylinders mit einer von z unabhangigen 
Ladungsdichte und das (ergiebigkeitsfreie) magnetostatische Feld eines strom- 
durchflossenen Zylinders mit von z unabhangiger Stromdichte [Z|. Die mathe- 
matische Physik kennt jedoch auch Beispiele harmonischer Vektorfelder, die im 
allgemeinen weder ergiebigkeits- noch zirkulationsfrei sind, bei denen also beide 
Potentiale im allgemeinen notwendig mehrwertig sind. Dazu zahlt das Aus- 
blasungs- und Absaugungsproblem der Aerodynamik [2]. Man sucht darin die 
inkompressible Potentialstromung zu bestimmen, die sich auf und auBerhalb 
eines 1m Parallelstrom befindlichen Tragfliigelprofils unendlicher Streckung aus- 
bildet, wenn auf dem Profil eine ,,Ergiebigkeitsverteilung‘‘, d.h. die Normal- 
komponente vorgeschrieben ist. Hier werden beide Potentiale? in dem Falle 
mehrwertig, wenn auBer der Ergiebigkeit auch die den Auftrieb erzeugende 
Zirkulation von Null verschieden ist. 

‘ Unter einem harmonischen Vektorfeld wird ein divergenz- und rotationsfreies 
Feld verstanden. Es folgt hieraus (in kartesischen Koordinaten) zwar, da® jede 
Komponente eines solchen Feldes auch eine harmonische Funktion ist, d.h. der 
Differentialgleichung von Laplace genitigt, jedoch ist nicht jeder Vektor ein harmo- 
nischer, dessen Komponenten harmonische Funktionen sind. 


» Die einzige Komponente des Vektorpotentials wird in der Strémungslehre als 
Stromfunktion bezeichnet. 
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Von besonderer Bedeutung ist die GréBe des betrachteten Vektorfeldes auf 
der Berandung selbst, z.B. die elektrische Feldstarke auf der Oberflache eines 
Leiters oder der Vektor der Strémungsgeschwindigkeit auf der Oberflache eines 
Tragfliigels, und zwar teils, weil von der physikalischen Problemstellung her 
oftmals nur die Randkomponenten des Feldes interessant sind und teils, weil das 
Feld im AuBengebiet auf Grund bekannter Satze der Potentialtheorie durch die 
Randwerte bestimmt ist. Die Untersuchung harmonischer Vektorfelder liBt sich 
daher zuriickfiithren auf das aquivalente Problem, die ,,Vertraglichkeiten‘‘ der 
Randkomponenten (Normal- und Tangentialkomponente) zu bestimmen, die, wie 
da man weif, nicht voneinander unabhangig sind. 


In der vorlegenden Arbeit werden ebene harmonische Vektorfelder im AuBen- 
gebiet beliebiger Berandungen, die im wesentlichen nur gewissen Stetigkeits- 
und Differenzierbarkeitsvoraussetzungen geniigen miissen, betrachtet. Ergiebig- 
keit und Zirkulation diirfen grundsatzlich von Null verschieden sein. Es wird 
der Zusammenhang zwischen der Normal- und der Tangentialkomponente auf der 
Berandung untersucht, und im Falle der Vorgabe einer Komponente werden 
entsprechende Existenz- und Eindeutigkeitssatze formuliert. Im Hinblick auf 
praktische Anwendungen wird ein numerisches Berechnungsverfahren angegeben. 
SchlieBlich wird gezeigt, daB sich diese Uberlegungen ohne Schwierigkeiten auf 
die in der Hydrodynamik vorkommenden ,,Gitterprofile“ tibertragen lassen. 

Die erhaltenen Ergebnisse kann man auf folgende einfache Weise von der 
Anschauung her verstehen. Sei eine beliebige Ergiebigkeitsverteilung auf einem 
Kreis vorgeschrieben, so laBt sich die Gesamtheit aller im AuBengebiet ein- 
schleBlich des Unendlichen méglichen, inkompressiblen Potentialstr6mungen 
bekanntlich aus vier ,,Grundstr6mungen“ linear superponieren; diese sind im 
einzelnen eine im Unendlichen verschwindende ,,reine Zirkulationsstrémung“, die 
den Kreis als Stromlinie enthalt, ferner je eine zirkulationsfreie, den Kreis als 
Verzweigungsstromlinie enthaltende Strémung, die im Unendlichen bzw. nur eine 
x- oder y-Komponente besitzt, und schlieBlich eine im Unendlichen verschwin- 
dende, zirkulationsfreie Strémung, die den Kreis in der Weise durchsetzt, daB 
die vorgegebene Normalkomponente angenommen wird. Auf dem Kreis selbst 
sind die zugehérigen Tangentialkomponenten dieser vier Anteile hinlanglich 
bekannt bzw. mit Hilfe des Poissonschen Integrals geschlossen darstellbar [3]. 
Es wird nun gezeigt, daB die fiir den Kreis bestehenden Zusammenhdange auch fiir 
beliebige Berandungen unverandert giiltig bleiben. Allerdings sind die praktisch 
interessierenden Tangentialkomponenten der genannten vier Grundstrémungen 
nicht mehr in geschlossener Form angebbar, wohl aber als Lésungen von Inte- 
gralgleichungen zweiter Art mit stetigem Kern ohne Schwierigkeiten numerisch 
zu gewinnen. 

§ 2. Die betrachtete Funktionenklasse 

In der (x, y)-Ebene sei eine stetig gekriimmte, doppelpunktfreie, geschlossene 
Kurve & gegeben. Es ist immer ohne Einschrankung méglich, fiir © einen Kurven- 
parameter y derart anzugeben, daB x(p), y(g) fir alle y zweimal stetig diffe- 
renzierbar sind, die Periode 2 besitzen und © mit wachsendem 0S yS2z genau 
einmal im Gegenuhrzeigersinn durchlaufen und da8 auBerdem die ersten Ab- 
leitungen x(q), y(~) nirgends gleichzeitig verschwinden. 
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Gegenstand der Betrachtung sind alle auf und auBerhalb © erklarten Vektor- 
funktionen v(x, y) = (v,(, ¥), V(X, y)) mit folgenden Eigenschaften: 


y ist auf und auBerhalb © stetig; 
» ist auBerhalb © stetig differenzierbar und harmonisch, d.h. auBerhalb © 


gelten die Differentialgleichungen 


: OVz Ovy 
= - ee 5 
(1) div t aE —. By 
= OVy OV x i iy Ds 
(2) Tot, o = a a Os 


» konvergiert im Unendlichen gleichmaBig gegen eine vorgegebene vektorielle 
Konstante 0.= (Uc, x) Voo,y), d-h. fiir alle OS OS 27 existieren die beiden Grenz- 


werte 
(3) jim» (R cos O, Rsin 0) = 


im Sinne gleichmaBiger Konvergenz. 


Sei v, die Normalkomponente und v, die Tangentialkomponente von 0, ge- 
nommen auf ©?2, so verstehen wir unter der Ergiebigkeit von © das Integral 


(4) Bisel aid ss 
€ 
wahrend die Zirkulation um © gegeben ist durch 


(5) Bl pope? 
G 


Stellt insbesondere » eine inkompressible Potentialstré6mung mit der Dichte 0 
dar, so sind Ergiebigkeit und Zirkulation maBgebend fiir die auf das Profil € 
ausgetibte Gesamtkraft & (Dimension: Kraft/Lange); es gilt namlich als Folge 
des Impulssatzes 


(6) KR =—oEv, +o! [vq, f 3 


Demnach bewirkt negative Ergiebigkeit einen Widerstand, d.h. eine Kraft- 
komponente in Richtung von », und negative Zirkulation einen Auftrieb, d.h. 
eine Kraftkomponente in der aus »,, durch Drehung um + 2/2 hervorgehenden 
Richtung. Fiir E=0 geht (6) iiber in den Satz von Kutra- JouKOWSKI. 


‘ Obwohl hier nur ebene Vektorfelder » betrachtet werden, ist es gelegentlich 
von Vorteil, » durch eine identisch verschwindende z-Komponente formal zu einem 
dreidimensionalen Vektor zu erganzen; insbesondere besitzt dann der (dreidimen- 
sionale) Vektor rot » die denticch verschwindenden Komponenten rot, » und rot, v. 

* Wenn nicht ausdriicklich etwas Gegenteiliges gesagt wird, soll n der Finhers 
vektor i in Richtung der 4uBeren Normalen und t der im Gencuihircicersnm orientierte 

Tengentialembeitsveltor einer geschlossenen Kurve sein; wir verstehen unter der 
Normal- und Tonmentialkomponcute v, und v, bzw. die Skalarprodukte (1, ») und 
(t, ») und unter 4/énm und @/ét bzw. die Differentiationen in Richtung von un und t. 
Die Einheitstangente t geht aus der Einheitsnormalen n durch Drehune um + 2/2 
pee ae 

* i, J, £ sind die Einheitsvektoren in x-, y- und z-Richtung. 
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§ 3. Darstellung harmonischer Vektorfelder durch ihre Randwerte 


Lemma. Es sei & ein p-fach zusammenhiingendes Gebiet der (x, y)-Ebene mit 
stiickweise glatten Berandungen ©;, i1=1...p. Einheitsnormale n und Einheits- 
tangente t der Berandungen seien so orientiert, daB n aus © herausgerichtet ist und 
© ber Umlaujung in Richtung von t zur Linken liegt. Fiiy alle Aufpunkte (x, y) 
aus % und alle Integrationspunkte (&,) == (x, y) aus © und G;, i=1... p, set 
F(x, ¥;&,n)=f(x—&, y—n) als stetige Funktion erklart; fir alle (é, n) += (%, ¥) 
aus & set F ferner zweimal stetig differenzierbar und harmonisch. Ferner soll 


| F’=F— In - fiir alle (E,n) aus einer Umgebung von (x, y) stetig differenzierbar 


sein, wenn v den Abstand zwischen Aufpunkt und Integrationspunkt bedeutet. Dann 
laBt sich ein in & und auf S;, 1=1... p, stetiges und in & stetig differenzierbares 
und harmonisches Vektorfeld »=(v,, vy) in ® gemaB 


(7) py = — grad @ + rot (Vf) 


zerlegen in den negativen Gradienten eines (einwertigen) skalaren Potentials ® und 
in die Rotation eines (einwertigen) Vektorpotentials Y £; dabei sind 


Pp 
1 
(8) ee f Fas, 
aig 
p 
(9) a Greg), [Fas 
aig 


in S beliebig oft stetig differenzierbare und harmonische Funktionen und v,, =(n, v), 
v,= (t, v) die Randkomponenten von v. 


Beweis. Uberfiithrt man das Gebiet & im Falle #>1 durch #—1 Schnitte in 
ein einfach zusammenhangendes, so ist die Behauptung fiir £>1 sicher richtig, 
wenn sie fiir /=1 gilt, da die tiber die Schnittberandungen erstreckten Integrale 
gemeinsam keinen Beitrag liefern. Es gentigt daher, den Beweis fiir ein einfach 
zusammenhangendes Gebiet & mit stiickweise glatter Berandung © zu erbringen. 
Aus der einmaligen stetigen Differenzierbarkeit von » in © und der Divergenz- 
und Rotationsfreiheit folgt die Existenz von stetigen Ableitungen beliebig hoher 
Ordnung aus der Tatsache, da8 w=v,—zv, eine analytische Funktion ist, da 
div »=0 und rot,v=0 mit den Cauchy-Riemannschen Differentialgleichungen 
identisch sind. AuBerdem sind Real- und Imaginarteil von w und damit die 
Komponenten von » in © harmonische Funktionen. Der in & liegende Aufpunkt 
(x, y) werde von einem hinreichend kleinen Kreis & umschlossen. Dann sei S* 
eine in ®& verlaufende, stiickweise glatte, doppelpunktfreie, geschlossene Kurve, 
die ® umschlieBt, ohne mit ® einen gemeinsamen Punkt zu besitzen?. Durch 
Anwendung des Greenschen Satzes auf das von G* und St berandete zweifach 
zusammenhangende Gebiet &* folgt 


(10) [ [AF — FAv) 48 dq ={(v 5-—F Sas [= Fo as. 
se 3 : 


1 Beziiglich der Orientierung von Normale und Tangente auf ©, ® und G* gilt 
das in § 2 Gesagte. 
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Zusammenfassung von (1) und (2) in 
grad v, = [f, grad v, | 


und skalare Multiplikation mit n und t ergibt 


OV, OVy OUz _ OUy 
Dii~. wha wheihent nat ae Yeu 
oder 
ov olf, v] 
(11) Op. 


Aus (10), (41), A4F=0 und Av=0 folgt nach Ausfiihrung einer partiellen Inte- 
gration 


[oH +t )ac— foH san F ano 
S* & 


Mittels der (durch skalare Multiplikation mit n und t unmittelbar zu beweisenden) 
Identitat 
OF 


rao grad F — v, [grad F, f} 


OF 
va +16») 


und der auf ®t giiltigen Gleichungen 


oe ca eames a” 
an y an’ at at 
ergibt sich 
i) OF OF’ 
[as fle | ar) oS 
R K 
=— bi aaind Beat [ x, grad F ds, t ; 
S* S* 


Strebt der Radius des Kreises St gegen Null, so konvergieren die Integrale der 
linken Seite bzw. gegen 27» und Null. Denkt man jetzt den Gradienten auf 
den Aufpunkt statt auf den Integrationspunkt angewandt (es tritt dann wegen 
F(x, y;&,")=f(x—&, y—m) ein Vorzeichenwechsel ein) und jeweils vor das 
Integralzeichen gezogen, so bekommt man 


2 tes ova di ira s.— igrad fu,F ds, f ' 
G* S* 


Die rechte Seite bleibt ungedndert, wenn man den Integrationsweg S* in S 
tiberfiihrt. Mit Hilfe der Identitat 


(42) rot (Yt) = [grad ¥% f] 
ergibt sich dann die Behauptung (7) bis (9). 


Satz 1. Ein Vektorfeld ») mit den in §2 genannten Eigenschaften ist diberall 
auBerhalb © vollstindig bestimmt durch seine Werte auf © und im Unendlichen 
und wird dargestellt durch 


(13) b= b.. — grad ® + rot (Vf) 
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mit 
(14) Cae i), In ds, 
PEGE iB 
€ 
(15) Y= [ yntds. 
2 26 a 
€ 


Beweis. Umgibt man den (auBerhalb © liegenden) Aufpunkt (x, y) mit einem 
hinreichend groBen Kreis © vom Radius R derart, daB © ganz im Inneren von 


S liegt, so gilt, wenn man F= In — setzt, auf Grund des Lemmas fiir das von 


S und © berandete, zweifach zusammenhangende Gebiet 


b = — grad (®* + ®) + rot ((W* + W) f) 


mit 


25 
€ 


p+ =— [v,n ds, @ sal fenntas, 
OGG - Y v 
S 


y+ — 1 fonds, i [ ein + as. 
PEGE: & Y PH | Y 
S 
Es geniigt dann bei Beachtung von (12) offenbar, die Giiltigkeit der Gleichung 
(16) — grad ®* + [grad Y*, f] =v, 


im Aufpunkt (x, y) zu bestatigen. Zieht man die Gradienten unter die Integral- 
zeichen, so gilt im Punkte (x, y) 


grad p* = — 35 Uy, F ds, 
S 
Se oldas 41 ue 

grad Y ; Up ds; 

folglich : 
, 1 1 
(17) —grad ©* + [grad Y*, = 1 ih (v,n+01)ds=5"5 fo & 
S S 


Die beiden fiir alle O<@<2z im Sinne gleichmabiger Konvergenz existierenden 
Grenzwerte 


(18) jim» (x + Ros, y + Rsin ©) = d.0 


sind eine unmittelbare Folge der Voraussetzung (3): Sei namlich k’>0 der 
Radius eines Kreises um den Ursprung derart, daB zu vorgegebenem e>0 fir 
alle R= R’ und alle OSX O<22 


|) (R cos. O, BR sit O) — 0,5, ,|"<, 


|v,(R cos O, Rsin@) — ve,|<e 
4* 
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a 


wird, und R= R’ der Radius eines Kreises um (x, y), der den zuerst genannnten 
Kreis umschlieBt, so gilt auch fiir alle R2R” und alle OSOS27 


|u,(x + Roos @, y+RsinO) ——,,|<e, 
|v, (x + Reos@, y+ Rsin@) — v~,,| <e. 
Es folgt daher mit (48) - 
Lint meee Hime t_[ 2 +Re0s@, y +Rsin@) dO = deo, 


R->-co 2%R R>o 2% 
0 


und der in (17) durchgefiihrte Grenziibergang R—oo liefert (16). 

Im Hinblick auf das Folgende betrachten wir in irgendeinem auBerhalb © 
gelegenen Punkt (x, y) zwei Einheitsvektoren n* und t* mit der MaBgabe, daB 
t* aus n* durch Drehung um + 2/2 hervorgeht. Dann folgt aus (13) in Verbindung 
mit (12) durch skalare Multiplikation 


o@D Ce 

(19) is: b) = (n*, Veo) Ree ot* ’ 
o@D ov 

(20) (Pn) (it alee eee 


Speziell erhalt man fiir n*=t und t*=j die Komponentenzerlegung von (13) ~ 


o@ OV 
21) ae ie Fog ee 

o@D as 
(22) Vy = Vo0,y ET ey = “tie . 


§ 4. Aquivalenzsatz 
Mit der vollstandigen Beherrschung des Vektorfeldes » durch seine Randwerte 
auf © und im Unendlichen ist freilich nicht gesagt, daB diese Daten beliebig 
vorgebbar sind. Vielmehr stellen wir gerade die Aufgabe, zu festem v,, alle mit- 
einander ,,vertraglichen‘’ Randkomponenten v,, v,; auf © zu ermitteln. Dies 
fiihrt im folgenden auf die Auflosung von linearen Integralgleichungen. 


Satz 2. Es sezen v, und v, zwer auf © stetige Funktionen. Dann ist das Bestehen 
irgendeiner der beiden Integralgleichungen 


1 Seg eee, 4° *o 1 
(23) vs + 5g | nds =2(n, Yoo) +2 3 f In—-ds, 
G G 
1 de O 
(24) oe oll v, In —ds = 2(t, 04) er} [ ein tds 
G € 


notwendig und hinreichend fiir die Existenz eines Vektorfeldes», das auger den 
m § 2 genannten Eigenschaften auf © die Normalkomponente v, und die Tangential- 
komponente v, besitat. Zu jeder Kombination v,,v,, die (23) oder (24) gehorcht, 
gibt es genau ein derartiges Feld»; dieses wird durch die Gleichungen (13) bs (15) 
dargestellt. 
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Beweis. Sei &* eine auBere aquidistante Nachbarkurve! zu € und seien n* 
und t* Einheitsnormale und -tangente auf €*. Dann gilt fiir das auf €* genommene 
Potential einer stetigen einfachen Belegung 9 auf © bei Annaherung von ©* an 
die AuBenseite von ©, die mit ©, bezeichnet werde, 


; 7d 1 Bal 4 
(25) climax fem ds= 2 [ ointds—ze, 
€ € 
; 6 1 Sad 4 Y 
. dim ap {ent as= zy f eln ds; 
€ € 


dabei sind die auf den rechten Seiten auftretenden Differentiationen unter den 
Integralzeichen ausgefiihrt zu denken und die Integrale selbst als Cauchysche 
Hauptwerte zu verstehen [3, 4,4]. Beachtet man, daB ® und YW in (14) und (15) 
als Potentiale der einfachen Belegungen v,/22 und v,/27 dargestellt sind, so 
gehen (19) und (20) in die mit (23) und (24) identischen Gleichungen 


2a On D, Qn Ot 


Tipe 1 
it )5) - [em as ie C frum ds, 
€ € 


UV (0.05) == sa} In : ds = a folntas anes 
€ € 

tiber, wenn sich der Aufpunkt (x, y) unbegrenzt dem Profil nahert. Fir den 
Beweis der Umkehrung erklaren wir durch (13) bis (45) ein iiberall auBerhalb © 
stetig differenzierbares und harmonisches Vektorfeld », wobei die v, und v, 
auf © stetig sind und irgendeiner der beiden Integralgleichungen (23) und (24) 
gehorchen. Da die Grenzwerte (25) und (26) fiir stetig gekriimmte Profile und 
stetige Belegungen 0 im Sinne gleichmabiger Konvergenz existieren [5], k6nnen 
wir v tiberall auBerhalb © und auf ©, als stetiges Vektorfeld erklaren. Die Ab- 
leitungen von ® und WY verschwinden im Unendlichen, so daB » dort gegen 0, 
strebt. Wir zeigen, daB diese Konvergenz eine gleichmaBige ist. Sei (€,7) ein 
laufender Punkt auf © und U der Umfang von ©, so gibt es Schranken M>0 
und M*>0, mit denen auf € 


M 
Psibea aes at 
Plea sae 
Ge Gals Chal Oe 
wird. Zu vorgegebenem ¢>0 gibt es dann immer ein R’>M derart, daf fiir 
alle R= R’ 
M* R+M _, 
R R—-M 


1 Kine zu einer gegebenen Kurve © aquidistante Nachbarkurve ©* wird definiert 
durch ihren Ortsvektor rept én 


wobei ry und n bzw. Ortsvektor und Normale von © sind und 6 eine von Null ver- 
schiedene, dem Betrage nach hinreichend kleine Zahl ist. Aus der geometrisch un- 
mittelbar einleuchtenden Tatsache, da entsprechende Punkte von © und @* einen 
gemeinsamen Kriimmungsmittelpunkt besitzen, folgt, daB sich die Kriimmungsradien 
genau um 6 unterscheiden und da somit €* eine stetige Kriimmung besitzt, wenn 
dies fiir © der Fall ist. 
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gilt. Fiir einen beliebigen Aufpunkt (R cos 9, R sin @) mit R2 k’ und 0OSOS22n 
folgt dann aus (14) 


op * atal y Vn eek COS Ow mr 
ax 0 ee | " R2—2R (Ecos O+ynsin O)+&+7? : 
€ 
oD ete poe R+M 3% 
éx|—~ 2x J. U R2—RM 
€ 
M* RAM € 


~ onR R-M ~ 27° 


da sich die Ableitungen 0D/dy, 6W/ox, 0W/0y entsprechend abschatzen lassen, 
liefern (12) und (13) 


|v — »,.| <|grad ®| + |grad ¥| < ores 


also auch Ran a haat 


|v, — V6,y| S| 0 —D|< é. 


Damit hat » die in § 2 genannten Eigenschaften. Es bleibt zu zeigen, daB » auf 
©, die Komponenten v, und v, besitzt. Bei Annaherung des Aufpunktes an © 
liefern (19) und (20) bei Beachtung von (25) und (26) 


: hy = 1 bee 1 0 “ae 
(27) lim_ (n*, 0) = (1, 0. 7] % n-dst+ +3. 5 [un 4s, 
€ 
i f2590 1 4 oe sir iY 
(28) lim (t, 0) = (t, ¥..) —5> 5, fo, Inds — fo In ds+ 4 
€ C 


Geniigen nun v, und v, der Integralgleichung (23), so liefert (27) 
h bi) aaa ee 
(20) lim (n*, v) = 94; 
ist dagegen (24) erfiillt, so 14Bt sich mit (28) der SchluB 
l panel ee 
(30) oo 
ziehen. Denken wir uns » durch (13) bis (15) auch im Innengebiet 8 und auf 
der Innenseite ©; von © als stetiges Vektorfeld erklart und beachten dabei auBer 
(25) und (26) noch die analogen Beziehungen fiir die Innenseite [3, 4, 5] 


(31) lim, sax [ elntds= 2 femtast+ze, 
€ € 
é 1 ee 
(32) HE ts fom tae, gln——ds, 
€ € 


wobei &** eine innere aquidistante Nachbarkurve zu © mit der Normalen n** 
und der Tangente t** ist, so erhalten wir mit Hilfe von (19) und (20) den 
Sprung, den » beim Durchgang durch € erleidet: 


] * a Seeeiare esse 
(33) ee ee Maree a an 
(34) lim (t*, ») — im (PERO) eet 


C*¥—> Cy > Cz 
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Je nachdem, ob nun (29) oder (30) erfiillt ist, liefert entweder (33) 


(35) lim, (1, v) =0 
oder (34) 
6) clim, (t**, 0) =0. 


Liege (%9, Yo) auf ©; und (x, y) in B oder auf ©,, so definiert das vom Wege unab- 
hangige Integral 
| x,y 
(37) u=fv,dx+u,dy 


%o, Vo 


eine in 8 und auf ©; stetig differenzierbare Funktion w(x, y) mit 
grad u=b, 


die tiberdies in 8 beliebig oft stetig differenzierbar und harmonisch ist und im 
Falle der Giiltigkeit von (35) die Eigenschaft 
(38) lim ,°“. =0, 


C* > ¢, On** 
besitzt ; anderenfalls ergibt (36) in Verbindung mit (37) 
(39) lm u=0O 


C** > C; 
(man braucht dazu nur den Integrationsweg von (37) entlang der Berandung G, 
zu fiihren). Anwendung des Greenschen Satzes auf das von ©** umschlossene 


Gebiet liefert 


[ [(@4u + grad? x) CRC aH Pm uot ds 
: or 
und damit durch Grenziibergang ©** — (; in jedem der beiden Falle (38) und (39) 


ffordxdy=o. 
B 


Also verschwindet v in 8 und auf ©, und (33) und (34) besagen, daB das Vektor- 
feld des AuBengebiets auf ©, die vorgegebenen Komponenten v,, und »v, besitzt. 
Zum Beweis der Eindeutigkeit nehmen wir an, es seien v, und bv, zwei Vektor- 
felder mit den geforderten Eigenschaften. Da diese Felder die gleichen Rand- 
werte v,, 2;, Ve. besitzen, sind sie beide durch (13) bis (15) darstellbar und somit 
identisch. 

Ist eine der beiden Komponenten v, und v, vorgeschrieben, so kann man 
zur Bestimmung der anderen grundsatzlich zwischen (23) und (24), d.h. zwischen 


einer Integralgleichung erster und einer solchen zweiter Art wahlen. Auf Grund 


der bekannten Tatsache, daB der Kern - 


ist, hat man hiermit die Wahl zwischen einer Fredholmschen Integralgleichung 
zweiter Art mit stetigem Kern und einer Integralgleichung erster Art mit singu- 
larem Kern. Indem wir uns anschlieBend fiir den ersten Weg entscheiden, d.h. 
v, bei vorgegebenem v, aus (23) und v, bei vorgegebenem v, aus (24) berechnen, 
stehen uns fiir Fragen der Theorie die Fredholmschen Satze und im Hinblick 
auf praktische Anwendungen einfache numerische Verfahren zur Verfiigung. 


1 ; re) fears s 
In A stetig, dagegen aR In =F singular 
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§ 5. Allgemeine Lésungen der Integralgleichungen zweiter Art 
Die Integralgleichungen (23) und (24) sind von der Form? 


10d 1 

(40) o+ tf fomds=t. 
€ 

Man weiB, daB die zu (40) transponierte homogene Integralgleichung als allgemeine 
Losung eine beliebige Konstante besitzt [4,6]. Infolgedessen hat (40) nach den 
Fredholmschen Satzen [3, 4,6] dann und nur dann eine auf © stetige Lésung o, 
wenn 7 stetig ist und zur allgemeinen Lésung der transponierten homogenen 
Gleichung orthogonal ist, d.h. der Bedingung 


(41) lards==0 
G 
genugt. 

Um diese Bedingungen im Falle der Fredholmschen Integralgleichungen (23) 
und (24) nachzupriifen, gentigt es offenbar, da (n,v,.) und (t,v,.) auf © stetig 
sind und im Mittel verschwinden, die zu einer auf © stetigen Funktion @ ,,kon- 
jugierte“’ Funktion (Cauchyscher Hauptwert!) 


ERO 1 
(42) kat S [en as 
€ 


auf ihr Verschwinden im Mittel hin zu untersuchen. Zu diesem Zweck werde 
eine gewisse (auBere) Nachbarschaft von © mit einer Schar aquidistanter, stetig 
gekriimmter Kurven C¥ schlicht tiberdeckt. Der Index 6 gebe dabei den Ab- 
stand zwischen CF und € an. Die Orthogonaltrajektorien dieser Schar, die gerade 
Linien sind, gestatten es, auf den CF einen Kurvenparameter s einzufiihren, 
der der im Gegenuhrzeigersinn wachsenden Bogenlange an der entsprechenden 
Stelle von © entspricht, jedoch im allgemeinen nicht mit der (ebenfalls im 
Gegenuhrzeigersinn wachsenden) Bogenlange s* auf CX identisch ist. Damit haben 
wir in einer (AuBeren) Umgebung von © ein krummliniges, orthogonales Kurven- 
netz (0, s) erklart, das € als Parameterlinie 6=0 enthalt. Sei d/d¢* die tangentiale 
Differentiation auf ©%*, so ist 


Pie ge ae: ; 1 
= oe | on ds 


€ 

eine stetige Funktion von 6 und s fiir 6>0. Da itberdies k* gemaB (26) fiir 
60 gleichmabig fiir alle s gegen den durch (42) gegebenen Grenzwert k strebt, 
ist dieser selbst auf © stetig. Bezeichnet ds/ds* das Verhaltnis der Bogenlangen 
von © und GC} an Stellen mit gleichem Kurvenparameter s, so folgt 


nih. asia { 
(43) h* = ae fen as. 


Cte 
¢ 


1 Auf diese Integralgleichung wird man z.B. auch bei der Behandlung des inneren 
Neumannschen Problems gefiihrt, wo man nach in 8 harmonischen Funktionen sucht, 
die auf © die Normalableitung at besitzen. Die Lésungen dieses Problems werden 


durch alle Potentiale von einfachen Belegungen o dargestellt, die der Integralglei- 
chung (40) gentigen [4, 5, 6]. 
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Nun ist ds/ds* eine stetig differenzierbare Funktion von 6 und s, da alle Para- 
meterlinien des (0, s)-Netzes stetig gekriimmt sind. Integrieren wir (43) fiir festes 
O>O0 iiber das Intervall O0<s<U, wobei U der Umfang von © ist, so wird 
daher nach einer partiellen Integration 


U U 
: 1 unas 1 
(44) [ ktas= s { tae} ( [elm 5 as)as. 
0 0 


€ 


Hier konvergieren auBer k* auch die Faktoren des rechts stehenden Integranden, 
insbesondere das (gleichmaBig stetige) Potential einer einfachen Belegung [5], 
gleichmabig fiir alle O<s<U, wenn 6-0 strebt; da ds/ds* auf © den Wert 1 
besitzt, gilt 

a Ws 


hen = 
se os ds* uy 


so daB der Integrand rechts in (44) gleichmaBig gegen Null konvergiert. Wir 
diirfen daher in (44) die Grenziibergange 6->0 mit den Integrationen vertauschen 
und erhalten somit das gewiinschte Resultat 


(45) Neds 
€ 


Die Existenz einer Lésung fiir die Integralgleichungen (23) und (24) bzw. 
fiir (40) ist damit gesichert. Die allgemeine Lésung setzt sich dann bekanntlich 
aus der allgemeinen Lésung der homogenen und einer speziellen Lésung der 
inhomogenen Gleichung zusammen. Die Mannigfaltigkeit der Lésungen 6 der 
homogenen Gleichung 
(46) ee fain ds=0 

ae ODP s Y 


entspricht nach den Fredholmschen Satzen der der transponierten homogenen 
Gleichung, sie geht also aus einer nichttrivialen Lésung von (46) durch Multi- 
plikation mit einem beliebigen konstanten Faktor hervor. Weiter gilt, daB eine 
im Mittel verschwindende Lésung von (46) notwendig die triviale ist. Infolge- 
dessen gibt es genau eine durch 


(47) fods=1 
Ly 


normierte Losung 6 von (46)!. Aus dem gleichen Grunde existiert genau eine 
spezielle Lésung o der inhomogenen Gleichung (40), die der Nebenbedingung 


(48) fods=0 
€ 
gehorcht. Die allgemeine Lésung von (40) lautet damit 


(49) o = const eG. 


1 Tm folgenden verstehen wir unter 6 die einzige Losung von (46) und (47); die 
allgemeine Lésung von (46) lautet dann const 6. 
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Bezeichnen k, und k, die durch (42) definierten Konjugierten zu v, und v 
und bedeutet ein Akzent die Differentiation nach der im Gegenuhrzeigersinn 
wachsenden Bogenlange s, so lauten (23) und (24) 


qe il iis 1 
(50) a cared fo, In dS 25, 3 Y= 2 op gh ea, 
€ 
1 4) ] 1 dees , , ZB 
(51) U;- re an Us er S = LU 0 x % + 2% 0.4 Y ~ mpg 


€ 


Die im Mittel verschwindenden speziellen Lésungen der inhomogenen Gleichungen 
(50) und (54) lassen sich jeweils in drei im Mittel verschwindende Bestandteile 
zerlegen, die bzw. von den drei Summanden der rechten Seiten herrtihren. Da 
die allgemeine Lésung die Form (49) besitzt, werden alle Normalkomponenten 
im Falle einer vorgegebenen Tangentialkomponente dargestellt durch 


0 1 2 3 
(52) U, = Ev, + Ve0,2Un F Vo0,y Un + Un 


ne _ (i ee Pelee : c Y 
wobei die ,,Grundlésungen“‘ des Problems v,,, v,,, v,,, U,, eindeutig bestimmt sind 
durch 


0 1 0 1 0 
€ € 
1 1 (4) 1 4 1 1 
i, ++ 5 fo, in ds =2y, fP.ds=o, 
(53) ‘ : 
2 2 
Vy = Bs ES fe In ar ag 25"; Uy ds= ? 
wu On Y Ay 
€ € 
3 es, S 3 
Pra nds=k,, [¥,ds=o; 
€ € 


entsprechend lauten im Falle einer vorgegebenen Normalkomponente alle hiermit 
vertraglichen Tangentialkomponenten 


(54) Dare Ten + Voo x2; Ar tess dv, aie v, 


mit den Definitionsgleichungen fiir die Grundlésungen 


0 ic 0 1 0 

vO) — = — == — 

+ Gq f Hints 0, feds i 
€ € 

1 Ae 1 1 ; 1 

U; => ——— — —— " ) = 

i eter fein Vas Don. [#ds=o, 

(55) 2 € 

2 AiG: 2 

U.-- —— Uv _ = i == 

tn On | By tos: haat hue [eas : 
€ € 

3 hg) 

w+ls fiintds=—hy, [i,as= 
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_ DaB die willkiirlichen Konstanten, deren Vorgabe das Problem jeweils eindeutig 


macht, in den allgemeinen Lésungen (52) und (54) die Bedeutung der Ergiebigkeit 


_£ und der Zirkulation I’ erhalten, folgt durch Integration dieser Gleichungen 
_tiber © unter Beachtung von (4) und (5) sowie der Nebenbedingungen, denen 
die Grundlésungen gemaB (53) und (55) unterliegen }. 


Die Grundlésungen lassen sich folgendermaBen anschaulich beschreiben: 


0 d ieee Ss ; a oe : 
ee v,: ,reines Ergiebigkeitsfeld‘‘, das © als Potentiallinie, bzw. ,,reines Zir- 


kulationsfeld‘‘, das © als Feldlinie enthalt, im Unendlichen verschwindend; 


v,, vy: ergiebigkeitsfreies Feld mit € als Potentiallinie bzw. zirkulationsfreies 
Feld mit © als Feldlinie, beide im Unendlichen gegen (1, 0) gehend; 
Dp, Oy: dto., nur im Unendlichen gegen (0, 1) gehend; 


vy, Uy: ergiebigkeitsfreies Feld mit vorgeschriebener Tangentialkomponente 
bzw. zirkulationsfreies Feld mit vorgeschriebener Normalkomponente, beide im 
Unendlichen verschwindend. 


§ 6. Transformation auf einen beliebigen Kurvenparameter 


Wir denken uns © durch eine Parameterdarstellung x(q), v(m) mit den in § 2 
genannten Eigenschaften gegeben. Nach Einfiithrung neuer Funktionen ,, @,, 


a a 
ee 415, Ons 0 0»... 3; durch 


(62) 
56 Un = ISL , Us a ) 
Be j#+R > +H 
# + 
Hl. Re = neue 5 Re ees 5 
7 Paes ee eas), 
+ + 
jos fog 
8 Z, == Ss aebeb 6 v = ae 5 — Or 3 5 
ies [REDE eo PRs 
+ + 


transformieren sich (53) und (55) in 


2 2a 


pee 0 0 
CD) —5— | K(9, y) Wn (y) dy =0, i wo, (y) dp =1, 
0 0 
Qa 27 
os 1 ° fia 
on (9) — == [ Kp.) np) dy =25(9), | On(p) dp = 0, 
(59) Ads ea 
> 2 : a 
u() — ste { K(p.¥) only) dy=—24(9), [ Only) dy=0, 
eee es 
t 3 3 
dO —a [Kev ow dv=ay), — f anlp)de=o, 
0 0 


} : 3 NEE Oe ; 
1 Die Dimensionen der Grundlésungen sind nicht einheitlich; v,, v, stellen eine 
Tee op oe. ; ’ n° & , : 
reziproke Lange dar, v,, 04, Un, U4 sind dimensionslos und v,, v, besitzen dieselbe 
Dimension wie |v). 
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bzw. 2m Ae 
OP) — 5  [ Ke. y) oo (y) dp = 0, | ele) d= 
yr os 
lg) — ah, [ Kp.) duly) dp =24(9), | ole) dp =o, 
(60) ve a 
on(®) — 5— [ Kg.) axty) dy = 29 (9), fo (p) dp=0, 
os ree 
ole) —- [Ky valy)dy=—malp), [ole 4p=0, 
0 0 


mit dem gemeinsamen Kern 


_ 5 ¥(9) [x(y)—*y)]—4 (9) YP H)) 
Enola [x (p) — * () 2+ Ly (gy) —y (y)]? 


Dieser Kern ist mit 2m periodisch und auBerdem stetig; denn es existiert der 
Grenzwert 


i K AGING) = IAP) b= Onde ee 

2h) yl, Pe Oe = eae amaRet « i 

Die Konjugierten k,, k, transformieren sich gemaB (42), (56) und (57) in 

20 
(63) x(~) = — ss L(y, yp) oy) dy, 
0 
wobei w durch w, bzw. w, und x durch z,, bzw. x, ersetzt zu denken ist und 
re Sea: ete (p)I+y et g)—y (p)) 

(64) (Q, Y) [x (@ wy) |? a lv (@ )-—y (y ye 
ein mit 22 periodischer, fir p=g+ ee 4; +2... ssingularer Kern ist: 


Das als Cauchyscher Hauptwert zu verstehende Integral (63) lat sich einfacher 
auswerten, wenn man von L(, y) einen stetigen Anteil H(g, w) abspaltet: 


(65) L(y, y) =— ctg 2—* + Hig, 9); 


die Stetigkeit von H(q, wy) folgt aus der Existenz des Grenzwertes 


(66) lim: ile pee I Pe eee ee 


Ypspt2kn Le (gy) ]?+ [v (p) |? 


Elimination von L(g, wy) aus (63) und (65) ergibt 
1% 22 
; | = a 
(67) xl) = 35 f oly) ete *% dy — 1 f Alp, y) oly) ay; 
0 0 
damit ist die Bestimmung der x,, x, im wesentlichen zuriickgefithrt auf die 


bekannte Aufgabe, zu einer stetigen und mit 2 periodischen Funktion f(p) die 
zugehorige ,,konjugiert harmonische‘‘ Funktion 


(68) sy) =s5 [ Hyetg *5% dy 


Harmonische Vektorfelder 61 


_ aufzufinden. Das _,,Kotangensintegral‘ (68) beherrscht man sowohl in der 
Theorie (die Fourier-Entwicklung fiir g(p) erhalt man aus der fiir /(m), indem 
man die Konstante fortlaBt und cos wy durch —sin w@ und sin wp durch cos wg, 

== 1,253 5.., ersetzt) als auch in numerischer Hinsicht, worauf in § 7 nadher 
{ eingegangen werden soll. 


Die stetigen Kerne K(g, y) und H(g,y) lassen sich auch folgendermaBen 
darstellen: 


> LS y(¢)—y (y) 
(69) K(g,y) = oy arc tg CS 
1g [x (p) — * (w)]?+ Ly (yp) —y (y)]? 
Sri 2 — 
D 


woraus die fiir das anschlieBend beschriebene numerische Verfahren wichtigen 
Identitaten 


27 
(71) JK.) dp = 20, 
20 
(72) tAior yd p— 0 
folgen. 


Ist speziell © ein Kreis (x= % +a cos(p — Yo), Y= Vo +4 sin(y — gp)), So wird 
K(y, y)=1 und H(g, y)=0. Dann liefern (56), (59), (60) und (67) die bekannte 
Tatsache, daB v,, (q) gleich der harmonisch Konjugierten zu v,(~) und 0, (9) gleich 
der harmonisch Konjugierten zu — v, (@) ist. 


§ 7. Diskrete Approximation der Grundlésungen 
Wir gehen von einer Parameterdarstellung x(), y(g) fiir © mit den in § 2 
genannten Eigenschaften aus. Es geniigt jedoch, daB die Funktionen x(q), v(q) 
an 2N diskreten Stellen auf © bekannt sind, die den aquidistanten Parameter- 
werten 
_ wn 


(73) Pu IN [p= Ore. (2N — 1), 


_ entsprechen (V=1). Ist eine solche Parameterdarstellung nicht explizit bekannt, 
sondern © z.B. graphisch gegeben, so kommt es offenbar nur darauf an, die Lage 
dieser ,,Stiitzstellen“ geeignet zu wahlen. Es sei betont, daB die Stiitzstellen 
beliebig und somit ganz den praktischen Erfordernissen der numerischen Rech- 
nung entsprechend auf © verteilt werden kénnen; die Verteilung braucht nur 
von einer gewissen EbenmaBigkeit! zu sein, damit die anschlieBend beschriebenen 
numerischen Operationen (z.B. die Differentiation) sinnvoll bleiben (Fig. 1). Fiir 
bestimmte, praktisch interessante Profiltypen kann man derart geeignete Profil- 
koordinaten auf einfache Weise unmittelbar geometrisch konstruieren [7]. 


1 Als Kriterium hierfiir kann z.B. dienen, da8 das durch trigonometrische Inter- 


polation der Profilkoordinaten ¥(9,), V(P), M=0.-.. (2N—1), entstehende Profil von 
dem gegebenen (etwa im Rahmen der Zeichengenauigkeit) nicht unterschieden werden 


kann. 
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Zunichst ist zu einer an den Stellen (73) gegebenen und mit 22 periodischen 
Funktion f(p) die durch (68) erklarte harmonisch Konjugierte g(p) und die 
Ableitung f(g) naherungsweise numerisch zu bestimmen+. Mit Hilfe der Methode 
der trigonometrischen Interpolation fand GARRICK [8] 


2N-—1 
(74) = 2X Chute = ON), 2 Ni Aes 
mit 
(75) Ga eee Y= US, M2 =A), 


Entsprechend ergibt sich fiir die Ableitung [7] 


oe ea 
(76) = XP y hase MOAN = Thy 
mit 
(77) p= fue Ce es Fea eth 


2N-1 

(78) ds B=, 
u=0 
2N STs 

(79) igo 
u=0 


Zur Berechnung der x,,, x, formen wir (67) zunachst mit Hilfe von (72) um in 


x(9) = 3 [ oly) ctg *5% dp + tf Aly, 9) @(9) — H(,v) oy) lay; 


wertet man die auftretenden Integrale durch (74) bzw. durch die Rechteckformel 3 
aus, so ergibt sich 
2N—1 2N-1 


(80) #,= Sy Coe. (Hy joes Fe Ta. B=0...(2N —4), 
v= 


wodurch zweierlei erreicht wird: erstens werden die Diagonalelemente (66) des ° 
Kernes H(q, y), deren numerische Berechnung wegen der vorkommenden zweiten 
Ableitungen mit erheblichen Fehlerquellen verbunden ware, nicht bendtigt, und 
zweitens gilt wegen (74) und (78) 

2N-1 


(81) 2 kOe 
u=0 


‘Ist /(y) bzw. f/(p, y) in den vorkommenden Veranderlichen mit 22 periodisch, 
so wird im folgenden /(9,) = fy baw, 1(Pu» Pr) =fuy gesetzt. Naherungen zu f(g) bzw. 
/(p, y) werden mit /(p) bzw. f(g, yw) bezeichnet. 

* Fir w 22N ist f, durch f, = fyu—gn erklart. 

* Fur diese im Falle periodischer Integranden bekanntlich sehr genaue Quadratur- 
formel wurde kiirzlich von HAMMERLIN [9] eine Fehlerabschatzung angegeben. 
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Mit Hilfe von (65), (73) und (75) leiten wir einen Ausdruck fiir (80) her, in dem 


nur noch die Werte des (einfacher zu berechnenden) Kernes L(g, y) vorkommen. 
Es wird 


2aN-1, (—1)- ( Fate Lo ee 
i Prin: aay p= ELA eal = 
H, a ACU ah a Ds (LG), — lay, ©) 
ae, re 
2N—1 
mule? (u—Y) 7 (yv—p)a 
ae an 24 (>, ctg , w, ctg on 4, 
ve 
1 2N—1 Aer (v bw) a { 2N-1 ecb 
aN, & feat) O, cte- ae GaN 2; (Lew han, Oy) e 
Yo +u 
2=0...(2N —4); 
hieraus folgt mit Hilfe von (76) und (77) 
S , N52 
(82) %,, = a Seay CE — LO.) w=0...(2N —1) 
v=0 
Vo Ll 


Hierbei ist zu beachten, daB On nur formale Bedeutung zukommt; da namlich 
w(g) nur als stetig vorausgesetzt ist, braucht die Ableitung @(p) gar nicht zu 
existieren. Der praktische Vorteil, die Bezeichnung On dennoch zu verwenden, 
besteht darin, da8 man zu den bekannten Daten x,, y,, w, nur noch die Ab- 
leitung (76) bestimmen muB, die Formel (74) ftir die Berechnung der konjugiert 
harmonischen Funktion jedoch nicht mehr bendtigt. Die w, bzw. die w,,, oder 
w,, sind bei gegebenem v,,,, oder v,,, selbst’ Naherungen und aus (56) zu berech- 
nen, wo man %, und y, durch 5 und Vy zu ersetzen hat. 


Wir haben nunmehr Naherungen fiir die inhomogenen Seiten der Integral- 
gleichungen (59) und (60) zur Verfiigung. Diese Integralgleichungen fiihren nach 
einer einfachen Umformung mit (71) und nach anschlieBender Anwendung der 
Rechteckformel auf Systeme mit zunachst 2N-+1 linearen Gleichungen fiir 2N 
Unbekannte [10]; sie lauten im Falle einer vorgegebenen Tangentialkomponente 


{ 2N=L 53, oO ~ 0 { paler A { 
> N (A, wOnw ES, 5Ony) =O, oN bd ane eA 
v=0 u=0 
yeu 
eee age aot tees 2 Pert xh 
aN ; (A, wpOnw es Dy ‘) = 2, , aN On yw — O; 
v=0 u=0 
vty 
(83) 1 ya a eae 2 eo) 2 . { 2N—19 
ys (ey On ih AGF yn y) =e 2% 4 ’ 4 N On “ 0 ) 
2N ON ee 
v=0 L=0 
peu 
1 2N-—1 ~ 3 or 3 % 1 aha le: 
AT > (Iori: Onn pe ES plo Ord PE ws aN Ds Dhan, 
2N a aN ao 
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und im Falle einer vorgegebenen Normalkomponente 


pee ene ee 1 2 1 
_ a ny al ike ny — ee (60) = —_—_ 
aN oan (KG, Oty uv Oty) OF 2 a tu Qn 
v=0 u=0 
VEU 
Fi let ae - ; 4 tH 
=e! (Om Kyy Oty) = 2%, a Wi,= 0, 
2N é 
v=0 p=0 
ye 
va fo yee OB nee, x Wigs. 
ya (Koj Oi — Buy Op) —2V ps Si Wi,—O0, 
v=0 u=0 
va 
1 2N—1 2, 3 as 3 = 1 pole | 
ON >) (LS On, oa, ey) = Hy ws 2N O14.= 05 
v=0 u=0 
VEL 


wobei stets w=0...(2N—1). Die Systeme (83) und (84) enthalten sdmtlich 
eine lineare Abhangigkeit und sind daher nur scheinbar iiberbestimmt: Die 
Summation iiber alle w=0...(2N—1) ergibt namlich wegen (79) und (81) 
identisch Null. Infolgedessen kann man ohne Einschrankung z.B. jeweils die 
Gleichung fiir ~=0 in (83) und (84) unterdriicken. Damit sind in jedem Fall 
vier Gleichungssysteme, bestehend aus 2N Gleichungen fiir 2N Unbekannte, 
zu lésen; diese Gleichungssysteme unterscheiden sich nicht durch ihre Koeffi- 
zientenmatrix, sondern nur durch ihre inhomogene Seite. 

Die in (83) und (84) auftretenden Naherungswerte fiir die Kerne K(q, y) 
und L(g, w) ergeben sich aus (61) und (64) zu 


(BS) Re alread ta el pt SL ot ere eee 


-- 
(86) ie g=="2 Hy (Fy ~ ee Hs? $0 ..(2N =e ete 


Mit den Lésungen von (83) bzw. (84) beherrscht man schlieBlich durch (58) 
die gesuchten Naherungen fiir die Grundlésungen des Problems 


% % 
On ¢ Ory 


<—, %%,=— ) ROSOes FoR Ot eN— ie 


a 
(87) Un ae ——» = ~ 
VaR oR V#R+9 
f + 


§ 8. Numerisches Beispiel 


Zu einer vorgegebenen Konstanten a>0 sind z*=-a Fixpunkte der nach 
JouKowsk1 benannten konformen Abbildung 


(88) i (2 as <). 


z 


Ein ,, Joukowski-Profil im engeren Sinne“ erhalt man aus (88) bei Abbildung 
eines Kreises der z*-Ebene, der den Punkt z*=a auf dem Rande und z*=——a 
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im Inneren enthalt, in die z-Ebene. Ein solches Profil besitzt jedoch bei z=a 
eine ,,scharfe Hinterkante“, ist also nicht stetig gekriimmt und entspricht somit 
nicht den in § 2 genannten Voraussetzungen. Diese werden dagegen von einem 


_,,verallgemeinerten Joukowski-Profil“ erfiillt, das durch Abbildung eines Kreises 
_entsteht, der beide Fixpunkte im Inneren enthalt, z.B.: 


yk — en O28 (0,1 4 1,6) 1,2 cos P+0,4 sin Pp \ 


( 1,6 le 
(89) meet et 
| yt = a{o,4 ents 6) 222 || . 
/ 1, 


Dieser Kreis mit dem Mittelpunkt (— 0,24, 0,4a) und dem Radius (0,1 +1/1,6) a 
wird mit wachsendem 0X gyX2z einmal im Gegenuhrzeigersinn durchlaufen; 
der kiirzeste Abstand zum Fixpunkt (a, 0) wird fiir p=0 angenommen und 


_ betragt 0,1@. In der z-Ebene bekommt man daher fiir p=0 eine ,,abgerundete 
_ Hinterkante“. 


Mit einer Konstanten c>0 laute die Normalkomponente auf dem Kreis (89) 


Cc 


* —— 
(90) R= 


(1 — 0,7 cos p — 0,1 cos 2y — 1,8sing + 0,4sin 29), 


so daB die Ergiebigkeit E*=(0,1+]/1,6) ac betragt. Da man die Verhiltnisse 
beim Kreis geschlossen beherrscht!, ergeben sich die Grundlosungen in der 2*- 
Ebene zu 


java cere 

*  2na(01+)1,6) ’ 

ux a 2 1,2. sin Grr On COS p 

1,6 
41 
(9 ) je ay 1,2 cos P+0,4 sin p 
/'1,6 
3x 


Vv, = — (— 0,7 sin » — 0,1 sin 29 + 1,8 cos my — 0,4 cos 29). 


Durch Einsetzen von (89) in (88) und anschlieBende Zerlegung in Real- und 
Imaginarteil bekommen wir die Parameterdarstellung des Joukowski-Profils © 
(Fig. 1) 


a R24 R24 2 
x (~) “5 2 ye? > 


(92) 


ye y*24 y*2— GQ? 


y(@) = > ‘i HDL yD 7 


Bei der Transformation der Grundlésungen in die z-Ebene ist zu beachten, daB 


Cet. , Kop Nine ar 
fiir v* die zu 1 normierte Zirkulation, dagegen fiir v;" und v," die Verhaltnisse 


3 
im Unendlichen erhalten bleiben miissen. Die Transformation von wv; und 7; 


ist zundchst nur bis auf einen willkiirlichen konstanten Faktor bestimmt; sie 


1 Vgl. §1 und $6. 
Arch. Rational Mech. Anal., Vol. 5 5 
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LET LESSEE EESEEESSES SSIES ann OS 


CO 


= 
COON AMF WN = 


laBt sich jedoch im Falle nicht verschwindender Ergiebigkeit durch die Forderung 
E = E* eindeutig festlegen. Die Normalkomponente (90) transformiert sich 
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Tabelle. Koordinaten, Evgiebigkeitsverteilung und Grundlésungen 
SS ee ——— ee 


1,003 72944 
0,997 067 54 
0,958 53214 
0,89407013 


0,80795896 | 


0,703 42879 
0,583 218 64 
0,44995579 
0,306 37044 
0,15538841 
0,000 141 89 
—0,156071 67 
—0,309861 66 
— 0,457 83315 
—0,596681 77 
—0,72328691 
—0,83479631 
— 0,928 69613 
— 1,002 86000 
—1,05557070 
—1,085 506 54 
— 1,091 683 62 
— 1,073 34466 
—1,02978710 
—0,96013491 
— 0,863 092 24 
— 0,73679688 
—0,579045 31 
— 0,388 36315 
—0,16641745 
9,078 480 31 
0,328 96981 
0,56101920 
0,75197631 
0,888 894 92 
0,97062747 


| 


- Vu 
a 


— 0,002 632 34 
0,026 59062 
OO7 7312412 
0,144021 83 
0,22041844 
0,300 37568 
0,378 38693 
0,449 738 20 
0,510 56031 
0,557834 64 


0,58938314 | 


0,603 85276 


0,600 695 38 | 


0,580141 38 


0,543 16478 | 


0,491 43820 
0,427 27697 
0,353 57225 


0,273 71329 | 
0,191498 96 | 


0,11103651 
0,036 62270 
— 0,027 405 60 
— 0,076876 54 


—0,10809100 | 
—0,118308 53 | 


— 0,106 50910 
—0,074490 52 
—0,02813491 
0,021 903 66 
0,061 63036 
0,078 63289 


Cn 
c 


0,364 55230 
0,03287987 
— 0,044 35321 
—= 0,075: 73595 
— 0,096 618 90 
—0,114 54827 
—0,13140922 
—0,146954 53 
—0,15969286 
—0,167 23415 
— 0,166 52219 
—0,15405499 
—0,12610049 
—0,07887795 
— 0,008 654.03 
0,088 313 99 
0,216071 23 
0,37957648 
0,585 59860 
0,843 15222 
115031251 
95133833077 
1,78595509 
1,762647 69 
1,46203978 
1,11191840 
0,826006 90 
0,61497448 
0,464 14673 
0,358489 34 
0,28646980 
0,23972911 


0,069009 11 0,21258381 

0,04085216 0,20229516 

0,010185 59 0,21162528 

— 0,007 400 32 0,263 55763 
Pr : tt 


| 
| 


0 
AU y 


Wey: 
0,58889 
0,34975 
0,26447 
0,223 52 
0,200 99 
0,187 88 
0,18030 
0,17635 
0,17503 
0,17581 
0,17844 
0,18285 
0,18917 
0,197 68 
0,208 85 
0,223 46 
0,242 66 
0,268 15 
0,30210 
0,346 13 
0,395 74 
0,426 24 
0,398 $7 
0,324 08 
0,249 54 
0,193 66 
0315537 
0,13038 
OnaS.o8 
0,10975 
0,113 14 
0,128 77 
0,16516 
0,24842 
0,492.06 


Fig. 1. Verallgemeinertes Joukowski-Profil mit 2N = 36 Konturpunkten 


0 
avs, 


1,33470 
0,58808 
0,35004 
0,264.85 
0,223 26 
0,201 11 
0,187 84 
0,180 33 
0,176 34 
0,17504 
0,17581 
0,178 44 
0,18285 
0,18917 
0,19768 
0,208 85 
0,223 46 
0,242 66 
0,26815 
0,30211 
0,34613 
0,395 74 
0,426 24 
0,398 57 
0,324.08 
0,249 54 
0,193 66 
0,15537 
0,13038 
0,11563 
0,10975 
0,113 14 
0,128 78 
0,16517 
0,248 44 
0,49209 
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fiiy ein verallgemeinertes Joukowski-Profil (2N=36) 


1 A 2 2 3 5 
v v ) tu tu 
tu tu Utu Ut z ca 
3,621 68 3,62595 10,865 03 10,858 66 2,55187 2,55444 
0,74081 0,75476 4,995 66 4,984 08 0,997 51 0,998 73 
— 0,081 92 —0,081 40 2,998 32 2,99975 0,51620 0,514 76 
—0,45471 —0,46279 2,22204 2,22745 0,332.91 0,33247 
—0,70458 —0,700 66 ‘leases, 1,78013 O32 y/ 0,23289 
—0,90231 | —0,90367 1,468 69 1,46955 0,162 34 0,161 92 
—1,06902 — 1,068 62 4520555 1,205 36 0,104 42 0,104 63 
= HP Vo —1,21149 0,961 20 0,961 31 0,049 22 0,049 11 
—1,32995 —1,32982 O20 45 0,72014 —0,00845 — 0,008 40 
—1,42406 —1,42417 0,474 69 0,47471 —0,071 67 (OKO AN FAG) 
—1,49145 —1,49142 O;22 14157, 0,22118 — Oa 2Z02 — O14 202) 
—1,52969 —1,52972 0 C2 572) —0,041 78 —0,21986 —0,21986 
— 1,53630 —1,53632 —0,31438 — 0,31437 — 0,304 39 — 0,304 39 
— 1.50878 —1,50880 — 0,596 30 — 0,59629 — 0,393 84 — 0,393 84 
—1,44445 —1,44446 —0,88742 —0,88740 —0,48561 —0,485 61 
—1,34010 —1,34011 — 1,188 33 —1,188 32 —0,57646 —0,57646 
—=41519127 —1,19128 —1,50112 = {550414 —0,66277 —0,66277 
— 0,990 92 — 0,990 94 —1,83001 —1,82999 —0,74055 —0,74055 
—0,72721 — 0,727 23 — 2,181 63 —2,181 62 —0,805 20 —0,805 20 
—0,38004 —0,38006 —2,56282 —2,56280 —0,84992 —0,84992 
0,081 07 0,081 05 — 2,967 26 — 2,967 24 —0,861 11 —0,861 11 
0,68040 0,680 37 —3,32495 — 3,32492 — 0,807 76 —0,80777 
1,34358 1,34355 — 3,399 56 — 3,399 53 —0,63814 —0,63814 
1,78926 1,78923 —2,91238 —2,91235 —0,35350 — 0,353 50 
1,843 98 1,84395 —2,07947 —2,07945 —0,07979 —0,07979 
1,676 34 1,67631 —1,33031 1550 29 0,096 83 0,096 83 
1,46043 1,46041 —0,79080 —0,79079 0,18991 0,18991 
1,264 04 1,264 02 — 0,421 35 —0,421 33 O23 3127; 0,233 26 
1,10601 1,10599 —0,16401 —0,16399 0,251 27 OV2Z5 1027 
0,991 24 0,991 21 0,027 08 0,02710 0,25947 0,25947 
0,92207 0,922 04 0,188 69 0,188 71 0,268 55 0,268 55 
0,902 38 0,902 34 0,35664 0,35668 0,287 42 0,287 42 
0,94095 0,940 87 0,57809 0,57816 0,326 54 0,326 53 
1,05977 1,05960 0,93975 0,93991 0,404 75 0,404 73 
1,32434 | 1,323 84 1,668 81 1,669 24 0,57260 0,57258 
2,009 34 2,007 66 3,71079 3,711 96 1,04270 1,04287 
somit in 
dex 
(93) Vy, (¢) dz 6’ 
wahrend die Grundlésungen (91) in 
0 0x | dz* 
v;(@) aie dz. 6’ 
1 ly | dz* dz 
(9) = dz las Con 
(94) 2g) <be | det] | 4 
NG Y | dz \e| dz* \c0’ 
3 3 g* 
v, (9) = deale 
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iibergehen; dabei ist wegen (88) 


GF eit eel Mice bet 
(95) ae Ne 2 =) Vb ya)? — 4.0? x?” 
az 1 
(96) Gz \5 MD 
einzusetzen. 


Dieses Beispiel wurde fiir 2N=36 Konturpunkte mit der elektronischen 
Rechenmaschine G 2 (mittlere Rechengeschwindigkeit: 20 Operationen pro Se- 
kunde) des Max-Planck-Instituts fiir Physik und Astrophysik in Miinchen nach 
dem in §7 beschriebenen Verfahren numerisch berechnet. Die Ergebnisse 
(Tabelle) zeigen im Vergleich mit der strengen Lésung (94) eine Genauigkeit 
bis zu fiinf Dezimalstellen; in einer Umgebung von y=0, wo ungenauere Werte 
wegen der starken Kriimmung von vornherein zu erwarten waren, liefert das 
Verfahren immerhin noch eine den praktischen Erfordernissen geniigende Ge- 
nauigkeit. Die G 2 bendtigte fiir dieses Beispiel eine Rechenzeit von 3 Std 25 min; 
daran ist die Auflésung des linearen Gleichungssystems (36 Unbekannte, 4 in- 
homogene Seiten), die nach dem Eliminationsverfahren von GAUuss- JORDAN 
durchgefiihrt wurde, mit 1 Std 45 min beteiligt. 


§ 9. Erweiterung auf Gitterprofile 

In der Str6mungslehre interessieren die inkompressiblen Potentialstromungen 
auf und auBerhalb einer gitterf6rmigen Anordnung von Profilen [7]. Ein der- 
artiges ,,Gitterprofil bestehe aus unendlich vielen kongruenten, einfach zusam- 
menhangenden Gebieten %8,, k=O, +1, +2..., mit doppelpunktfreien, stetig 
gekriimmten Berandungen ©,, k=0, +1, +2..., die mit gleichem Abstand t>0 
(,,Gitterteilung‘‘) in y-Richtung angeordnet sind. Die ©,, R=0, +1, +2..., 
sollen sich gegenseitig weder bertihren noch iiberschneiden. Die Menge der auBer- 
halb des Gitterprofils gelegenen Punkte sei % genannt (Fig. 2). Fiir Cy existiere 
eine Parameterdarstellung x(g), v(y) mit den in § 2 genannten Eigenschaften. 
Dann ist eine beliebige Berandung ©, gegeben durch x(g), y(y) +t, R=O, +1, 
+2.... Die Konstanten 


(97) a, =Max{x(y)}, a_=Min{x()} 


legen die seitliche Begrenzung des Gitterprofils fest. 


In &% und auf den ©,, k=O, +1, +2..., sei ein Klasse von Vektorfeldern 
v(x, y) mit folgenden Eigenschaften Gegenstand der Betrachtung: 


y ist in M und auf den ©,, k=0, +1, 4 


2 on, SOLES 
» ist in YW stetig differenzierbar und harmonisch; 
y ist in Mund auf den ©,, k=O, +1, +2..., hinsichtlich y periodisch mit ¢: 


(98) v(x, y+?) =b(x, y); 
es existieren die Grenzwerte 
1D 90% 
(99) limi eat (2,9), =0 


*—> +00 
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fiir alle 7 eines beliebigen Intervalles y — Bs Si Hee im Sinne gleichmaBiger 
Konvergenz. ip - 

Ergiebigkeit und Zirkulation werden analog zu (4) und (5) definiert, wobei 
nur jeweils tiber eine beliebige Berandung, z.B. tiber Gy, integriert wird. 


Satz 31. Der tiber irgendeinen Periodenstreifen genommene Mittelwert von » 
ast vechts bzw. links vom Gitter gleich einer vektoriellen Konstanten 


Vio = (Os Ud coy) bzw. DE eo UE, ys Unico, 9) 


d.h. es gilt fiir alle a, >a, und a'_<a_ und alle y 


(100) : if Dla, 1) C70 ye 


Beweis. Betrachten wir ein rechts vom Gitter gelegenes einfach zusammen- 
hangendes Gebiet 8 mit doppelpunktfreier, stiickweise glatter Berandung G. 
Oberer und unterer Teil von © sollen kongruent sein und in entsprechenden 
Punkten mit gleicher Abszisse den Abstand ?¢ besitzen. Rechter und linker Teil 
von © seien gegeben durch die Geradenstiicke 


II\ 


t 
nSw++ baw. £=%, y—+<SnSyyt 


2 2 Ores 


lo 


E = x, At am 


wobei (%», Yo) und (%,, y,) zwei Punkte mit x,>%)>a, sind. Auf Grund des 
GauBschen und Stokesschen Integralsatzes gilt dann: 


Ws Yo 4 
Pa ree teas 1 iiieo ee 
+ff divudxdy=~ | V, dS =—— i Vy (%,%) dy, | v,(%g,) dn = ) 
: é a Ve 
t t 
‘ ie 2 Yor 5 
+ [[ rot.vdxdy= 7 | “ws = | Vy(%1,4) an — [ % (20,1) dn =0. 
y é Ane ae 
moe AE 


Vektorielle Zusammenfassung ergibt die Existenz eines », ., mit der Eigenschaft 
(100). Der Beweis fiir die Existenz eines vektoriellen Mittelwertes »_,, links 
vom Gitter verlauft analog. 


Satz 4. Die Differenz der Mittelwerte (100) ist durch die Ergiebigkeit und die 
Zirkulation fiir das Profil ©) bestimmt: 


(101) Ds oo — 9» =—I+- 


awe 
i t 


1 Dieser Satz wird in [7] unter der (weniger allgemeinen) Voraussetzung der 
Existenz einer einwertigen Stromfunktion bewiesen. 


70 Erich MARTENSEN und KuRT VON SENGBUSCH: 


Beweis. Durch Anwendung des Gaufschen und Stokesschen Satzes auf das 
zweifach zusammenhangende Gebiet &, mit den Berandungen ©, und ©, (Fig. 2) 
erhalt man 


Wt y+ 4 
J {I divudxdy = | v, (4,9) dn — [ va.) dn—+ fm ds=0, 
Gx phd jue Ge 
ty y+4 
1 1 , 1 , 1 
{ ff rotsndxdy = 7 i Vy (a4,4) an — - it Vy (a_,%) dy — + [ 4ds=0. 
Op Ws 5 get Cx 


Fig. 2. Gitterprofil 


Diese Gleichungen sind wegen (4), (5), (98) und (100) identisch mit der Kom- 
ponentenzerlegung von (101). 


Satz 5. In U wird v dargestellt durch 


(102) D =D. — grad @ + rot (YF), 
WOKE? Voo= (Voo, x» Voo, y) das vektorielle Mittel 
(103) oo = (04 co + D- oo) 
ist und 
(104) ca [ %,F ds, 
Pye o 
G 
(105) a [uFds 
2 0te 
©, 


harmonische und hinsichtlich y periodische Funktionen sind mit 


(106) Pye ee 


\2 (coy oe ae 2atyeD) 
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Beweis. Ein beliebiger in & gelegener Aufpunkt (x, y) laBt sich immer als 
innerer Punkt eines zweifach zusammenhangenden Gebietes (, mit den Beran- 
dungen ©, und ©, auffassen (Fig. 2). Dort erfiillt (106) die Voraussetzungen 


zum Lemma aus §3; insbesondere ist /’=F — ln + in einer Umgebung von 
(x, y) analytisch [7]. Es gilt daher 


(107) » = — grad (®* + ®) + rot ((W* + ¥) f) 
mit 
o*=— 1 |v, Fads, O= 3) [o,F ds, 
Dae. 20, 
Sk C;, 
fe ETT HRs sake F 
pr — sz TOUT «tigen fukas. 
ic Cy 


Wir untersuchen den Ausdruck 
(108) 1p = — grad ®* + [grad ¥*, f]. 


Bei der Auswertung der Integrale auf ©, liefern die Anteile des oberen und 
unteren Wegabschnittes infolge der Periodizitatseigenschaften von » und F zu- 
sammen keinen Beitrag. Wir kénnen daher (108) gemaB w=, -+tv_ in zwei 
Anteile zerlegen, die den Integrationen iiber die beiden Geradenstticke rechts 
und links vom Gitter entsprechen. Zieht man die Gradienten jeweils unter die 
-Integralzeichen, so findet man durch elementare Rechnung 


t / 
, ety Gin amen 
meee | 40d) aes Ae tis 
2t J Cos 22(x—a’,) ioe 22(y—n) 
Y4+— 2 t t 
sin ol 
ANU ae Ai) eee oe pea: 
+. Gof gees a.) ee 2m = ) 
Bei Beriicksichtigung von (100) folgt 
t 
ones 22a 
(109) Wp=Z% eZ ‘r (e ty (a',,9) P(x, ys @,0) + 
v4 


mit 
2a (x—a’,) 
2 é Fite cos Zeer 
UX 
/ at i 
(110) P(x, yy, a+, n) Che ; 22 (x—a',) 4n(x—a,) ? 
ACS ets Sa Fagee 
sin 2n(y—n) 

Qax She 

eer t 
45 on ge yee: == — : a 
ett) O(%, 95 40M) ta) on (yy) eet 
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* / ss vt / t iE 
Diese Ausdriicke sind wegen a’, > x fiirallea, 2a, undalle y, — = Sn a 
gleichmaBig beschrankt : 


22 (x—a’,) 


af 2u%x es: 
[PR VAa ia at en 
(1—e 
2a% 1 
OC a Cee ee 


Im Zusammenhang mit der Voraussetzung (99) konvergiert der Integrand in 
(109) fiir alle y, — 5 Sy SS) + gleichmaBig gegen Null, wenn a’, + 00 
strebt. Folglich existieren die Grenzwerte 
Lin | sy, Da 
al, —> + 00 
Nach einer entsprechenden Uberlegung fiir w_ ergibt sich 
lim. op lim (oF = F105 + Diss) = Dose 


a, —> +--00 @;.— + 00 
a’_—+>—oo a_—— oo 


Zusammen mit (12), (107) und (108) folgt die Behauptung (102) bis (106) bis auf 
den geringfiigigen Unterschied, daB die Integrale fiir ® und WY zunachst noch 
iiber ©, zu erstrecken sind. Auf Grund der Periodizitatseigenschaften andern 
diese Integrale ihre Werte jedoch nicht, wenn sie tiber Gy statt iiber ©, erstreckt 
werden. 


Satz 6. Im Unendlichen konvergiert v gleichmaBig gegen die konstanten Mittel- 


werte im Endlichen », ., a.h. fiir jedes Intervall y — « <ynSy+ + existieren die 
Grenzwerte - 


(112) lim 0(%,7) =V+0 


xX—>+ 
um Sinne gleichmapiger Konvergenz. 


Beweis. Wegen (101) und (103) gilt allgemein 
(113) Disco = Doo t it Fj. 
Bei Benutzung von (4), (5) und (12) folgt aus (102), (104), (105) und (113) 
(114) Die eee 5 | = UV, (= grad F + i + v, - grad F + t, tl} as, 


C, 


wobei der Gradient auf den Aufpunkt (x, vy) bezogen ist. Bezeichnet (€, 7) den 
Integrationspunkt auf ©), so wird mit (106), (110) und (411) 


2U% 


l a z a ; 
(115) , gradl +i=2e + {i P(éE,n; x,y) +j{Q(En; x, y)}. 
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| Mit einer Zahl a’. >a, erhalten wir aus (97), (110) und (411) fiir alle (&, 7) auf 
©), sowie fiir alle xa’, und alle y eines beliebigen abgeschlossenen Intervalls 
die gleichmaBige Beschrankung von P(&,; x, y) und Q(&, 7; x, y): 


2n(a,—a) 


27a, 
Peta Xx, y)| =e} ye 27 
ag é 
oma, 4 
OS) OR re rae 


Damit folgt aus (115) die Existenz des Grenzwertes 


(116) lim (= grad F + i) =0 

xX—>-+ co \ 7 
im Sinne gleichmaBiger Konvergenz fiir alle (€, 7) auf ©, und alle y eines beliebigen 
abgeschlossenen Intervalles. Durch Vertauschung des Grenziibergangs x —>-+ co 
mit der Integration tiber ©, liefern (114) und (116) die Behauptung (112) fiir 
x—>-+ oo. Entsprechend verlauft der Beweis fiir x—— ov. 


Satz 7. Es seten v,, und v, zwei auf Gy, stetige Funktionen und v,, eine vektorielle 
Konstante. Dann ist das Bestehen irgendeiner der beiden Integralgleichungen 
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(117) vy +2 af 0,F ds = 2(tt, Vgc) ++ z{ v,F ds, 
@ C, 

(118) yp [uFas=21, Deo) —— & [ %,Fas, 
G, , 


wober F die Bedeutung (106) hat, notwendig und hinreichend fiir die Existenz eines 
Vektorfeldes » in X und auf ©,, R=0, £1, +2..., das aufer den oben genannten 
Eigenschaften auf ©, die Normalkomponente v, und die Tangentialkomponente v, 
besitzt und dessen Mittelwerte v. . der Gleichung (103) gentigen. Zu jeder Kombi- 
NALLON Vy, Vz, Veg, Ate (117) oder (118) geniigt, gibt es genau ein derartiges Vektor- 
feld v; dieses wird durch die Gleichungen (102) bis (106) dargestellt. 


Beweis. Man kann zeigen, daB die Grenzwerte (25), (26), (31) und (32) unver- 


andert gelten, wenn man In durch F (definiert durch (106)) und € durch ©, 


ersetzt; dariiber hinaus kénnen auch die in §5 benutzten Existenz- und Ein- 
deutigkeitssdtze sinngema8 auf Gitterprofile tibertragen werden [7]. Der Beweis 
verlauft daher weitgehend analog zu dem fiir Satz 2, und wir beschranken uns 
auf die Darstellung der Besonderheiten. Dazu zahlt zunachst die Periodizitats- 
forderung (98), die wegen der Periodizitat von (106) gesichert ist. Durch Ein- 
setzen von (104) und (105) in (102) erhalt man nach einer Umformung 


p=t2ee | v,ds + or [ ds — 


(110) Ge G 
aii fa [ v, (< grad F4 iJds | a ie grad F +i, fds. 
6, G, 
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Es wurde bereits gezeigt, daB - grad +i fiir alle (€,7) auf Cy und alle y eines 


abgeschlossenen Intervalls gleichmaBig gegen Null konvergiert, wenn %—>-- 00 
2% 


geht, so daB nach Multiplikation von (119) mit e ¢* auch der Grenzwert (99) 
fiir x >-+co im Sinne gleichmaBiger Konvergenz existiert und den Wert Null 
hat. Das Entsprechende laBt sich fiir «——oo beweisen. Es bleibt schlieBlich 
zu zeigen, daB die Mittelwerte v,. von (119) Gleichung (103) gentigen. Dazu 
beachten wir, daB fiir (119) als Element der zugrundegelegten Funktionenklasse 
der Satz 6 gilt, so da8 wir die Mittelwerte p,,., durch die Grenziibergange 


x—>--co bekommen. Da auBerdem = grad F — i fiir alle (&, 7) auf © und alle y 


eines abgeschlossenen Intervalles gleichmaBig gegen Null konvergiert, wenn 
x—>— oo strebt, gestattet (119) die Grenziibergange > + co auszufiihren: 


5 


* ae 5 fey, 
Lo= Doky | v, aS Ses fads. 


C, G, 


Damit ist (103) erfiillt. 
Fiihrt man den Kurvenparameter y auf ©, ein, so unterscheiden sich die 
erhaltenen Integralgleichungen von den in § 6 betrachteten nur durch die Kerne 


am 2 (v) Sin anlage 4 () sin 2 — 
K(9,y) =7.- ~ ar! i 7 
(120) K(p,y) => 2x [x (9) —*(y)] 2x[y(9)—y (y)] 
Co] ; cos oe Y) 
on (0) Sin = Be ifolt oink 2a rel 
(121) L(y, y) = t .  2n[x(y)—x(y)] 22/v(o)— (y)] 
Col aarp pee er ae e Lah 3 


Definiert man H(y, y) wiederum durch (65), so kommen die zu (69) und (70) 
analogen Beziehungen 


[Vv (~) — y (p)] 
t 


[4% (p) — «(y)] 


tg 


: ie 
20 \ = an 
(122) K(g, yp) =2 a@ arc tg 


Ig 


; Cof 2nla(p—*y)] _ 6, 2a lylo—y ty] 
(123) FLX Gy) fee a he 
Op 7 — 


zustande, womit auch (71) und (72) Giiltigkeit behalten. Die Grenzwerte (62) 
und (66) gelten ebenfalls unverandert. Damit kann das in § 7 beschriebene 


numerische Verfahren wortlich iibernommen werden, nur sind jetzt K wy und 6 
bn 
Naherungen von (120) und (121). 


v 


Durch Anwendung des Impulssatzes auf die Berandung ©, des Gitterprofils 
erhalt man fiir die Gesamtkraft gleichfalls den Ausdruck (6), allerdings mit der 
veranderten Bedeutung (103) fiir ».. 


[9] 
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Poincaré and the Shortcomings 
of the Hamilton-Jacobi Method for Classical 
or Ouantized Mechanics 


L. BRILLOUIN 


Communicated by C, TRUESDELL 


Abstract 


A great theorem was proven by H. PorncarE in celestial mechanics. It states 
that, in the most general problems of mechanics, the total energy of the system is 
the only well behaved first integral of the system, while other so-called integrals ~ 
cannot be represented by uniform and convergent series. This very important result 
can be explained and visualized by comparison with standard methods of discussion, 
as, for example, the Hamitton-JaAcopr procedure. The discussion shows that there 
are serious limitations to the use of this procedure, which collapses in the most general 
problems (PorncarE theorem) and can be used only for “‘almost separated”’ variables. 

The PorncarE theorem appears to provide the distinction between determinism 
in mechanics and statistical mechanics according to BOLTZMANN. 


1. Poincaré’s Science and Hypothesis 


Famous old authors are always worth rereading. Let us open POINCARE’S 
book ‘‘Science and Hypothesis’! to Chapter VIII, where the author discusses 
the definition of “‘energy’’; we find on page 158 (French edition, or page 121, 
English edition) a very curious problem stated in the following way: 

“In the determinist hypothesis, the state of the universe is determined by 
an extremely great number m of parameters which I shall call %,, %,..., %,. 
As soon as the values of these parameters at any instant are known, their 
derivatives with respect to the time are likewise known and consequently the 
values of these same parameters at a preceding or subsequent instant can be 
calculated. In other words, these m parameters satisfy ” differential equations 
of the first order. 


“These equations admit of m—1 integrals, and consequently there are n — 1 
functions of %,, %2,..., %,, Which remain constant. If the we say there is some- 
thing which remains constant, we only utter a tautology. We should even be 
puzzled to say which among all our integrals should retain the name of energy.” 


' Porncarsé, H.: La Science et l’Hypothése. Paris: Flammarion 1902; English 
translation, The Foundations of Science, Science and Hypothesis, The Value of Science, 
Science and Method. New York: The Science Press 1929. 
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The author presents there, in a paradoxical way, a very profound problem. 
In a private letter (December 1958), SCHRODINGER writes: ‘“‘POINCARE is sort 
of pulling his reader’s leg, in order to explain in the simplest way some very 
new and relevant points of view ... There is not just one constant of the motion 


| but (2 —1). Which one of them is the energy? or rather, what function of those 


(x—1)? I believe these were points very worth making, and perhaps worth 
thinking about even now.”’ 

We intend to discuss the preceding statement, and we shall compare it with 
a very important theorem discovered by Porncarké? himself in his work on 
celestial mechanics, which contains a very pertinent answer to the problem. 
This will lead us into a discussion of the limits of accuracy in classical mechanics, 
and of the foundations of statistical mechanics. The problem was ‘originally 
intended by PoINcaRE to apply to a general physical system, but his statement 
certainly cannot be used for a quantized structure since it does not include any 
possibility of quantum jumps. Hence we shall restrict the discussion to a classical 
mechanical system. We thus examine the following problem: 

We consider a conservative mechanical system, with no dissipation (no friction, 
no passive resistance) nor regeneration (no external energy source, no negative 
(active) resistance). For such a model we easily define potential and kinetic 
energy, hence total energy too. This is the way all problems of classical mechanics 
are usually stated. PoINcARE himself, when discussing a celestial problem, has 
no hesitation how to write down the total energy of the system. 

There are N coordinates %,, %),..., Xy and N corresponding momenta #,,..., Py. 
This represents the m parameters of POINCARE: 


(1) n= QIN : 


The initial state of the system at t=0 is defined when we give the initial values 
%1(0),..., Xy (0), £1 (0),..., y(O) of the 2N parameters. The mathematician 
considers these initial values as ‘“‘given’’, but the physicist asks ‘“‘given by whom’? 
We must enquire how these initial values are measured and investigate the 
accuracy. This involves a complete change of point of view, and here we quote 
a remark of HADAMARD: “‘The question is not whether you make a very small 
change in the initial conditions, the real problem is whether this small initial 
change may not result in a very large change in the final results.’’ In other words, 
the problem concerns the stability of the trajectories. 

From t=0 on, the system moves on according to the laws of mechanics, and 
we may compute the coordinates x;(f) and momenta #;,(/) at time ¢. These 2N 
quantities are defined in terms of one variable ¢ and (2N —1) constants, to be 
computed from the initial conditions. These 2N —1 constants are the ‘‘first 
integrals’ considered by PoiNcaRE in his paradoxical statement in Science and 


Hypothesis. 
2. Poincaré’s great theorem on celestial mechanics 


Let us now open a previous book? of PoincarE himself (1892), the first volume 
of his ‘‘ Méthodes nouvelles de la Mécanique céleste’’, and read Chapter V, entitled: 


2 PoINcaRy, H.: Méthodes Nouvelles de la Mécanique Céleste. Paris: Gauthier- 
Villars 1892—1899; reprinted by Dover Publications, New York. 
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Non-existence of uniform integrals. The whole chapter is devoted to the dis- 
cussion of one fundamental theorem: 

The canonical equations of celestial mechanics do not admut (except for some 
exceptional cases to be discussed separately) any analytical and uniform integral 
besides the energy integral. 

This great Poincaré theorem contains the answer to the previous problem. 
It definitely proves that the (7 — 1) integrals previously introduced by POINCARE 
have very different properties. The total energy is the only expression which 
is represented by a “‘well behaved” mathematical function, and all other so-called 
integrals are not analytical; they may have a very strange behavior, with dis- 
continuities, and cannot be compared with the energy. The integrals representing 
momentum or moment of momentum are.comprised in the other ‘‘exceptional 
cases’’ discussed by POINCARE. 

For instance, PoINCARE applies his theorem to the three body problem (vol. I, 
p. 253). Assuming the center of mass to remain at rest, there are 4 actual first 
integrals, corresponding to moment of momentum and energy; POINCARE’S 
theorem proves that no other well behaved first integral can be found, and it 
thus puts an end to one century of unsuccessful attempts. 

Why did Poincaré choose not to mention this great theorem in his later | 
book? The reason must be that, in Science and Hypothesis, he had in mind a 
much more general problem and was thinking of the behavior of any arbitrary 
physical system; but we already said that his statement certainly does not hold 
for quantum theory and that the problem must be restricted to classical mechanics 
only. 

Our aim is to compare this theorem with the original statement of POINCARE, 
and in doing so, we have to present a short summary of the standard methods 
of Analytical Dynamics. We choose to follow an excellent textbook written by 
MAx Born? on Atom Mechanics, where the classical methods are very precisely 
discussed for the problems of the “‘old’’ quantized mechanics, for which they 
represented the skeleton of the theory. 

This text is also chosen because of a curious coincidence: after discussing very 
precisely the Hamilton-Jacobi method, and using it for computation on perturbed 
systems, M. Born considers the great Poincaré theorem, which apparently contra- 
dicts the “‘general”’ perturbation methods of M. Born and restricts their validity 
to the exceptional cases where the theorem does not apply. Born just quotes 
two theorems of BrRuNs and PoINcARE (p. 292). Both theorems refer to the 
convergence of expansions obtained by the perturbation theory. The theorem 
of Bruns (1884) states that points of absolute convergence or of divergence 
are equally dense. The perturbation procedure introduces non-analytic functions 
on one side and assumes the property of analyticity at another point: it thus 
contains internal contradictions. PorNcark&’s theorem also refers to perturbation 
methods in celestial mechanics and goes further than Bruns. BorN mentions 
these grave difficulties and does not elaborate further (he must have been too 


3 Born, Max: Vorlesungen iiber Atommechanik, erster Band, collection ,,Struktur 
der Materie“, vol. 2. Berlin: Springer 1925; English translation, M. Born, The 
Mechanics of the Atom. London: Bell 1927. Also C. Lanczos, Variational Principles 
of Mechanics. University of Toronto Press 1949. 
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busy at the time with the ‘‘matrix mechanics’’ he was inventing, with HEISEN- 
BERG). He nevertheless adds some comments: “‘It thus was never possible to 
prove the stability of the solar system ... The perturbation methods used in 
celestial mechanics are not convergent, but only semi-convergent (POINCARE, 
vol. II, Chap. VIII). We thus see that it is impossible, for purely theoretical 
reasons, to prove the absolute stability of atonric systems.” Let us add that these 
semi-classical methods completely failed in the case of the helium atom. 

For atomic structures, the solution was found with the invention of wave 
mechanics, which considers the old-fashioned quantized mechanics only as a 
first approximation of limited validity, as proved by the well known B. W. K. 
method. 

We shall first remind the reader of the fundamental methods of classical 
mechanics, in order to be in a position to discuss the importance of POINCARE’S 
theorem and to examine the limitations it introduces upon the validity of classical 
methods. 


3. The methods of analytical dynamics for separated variables 


As stated in the preceding section, we refer the reader to BorNn’s textbook 
for all details, and we shall be satisfied with a short summary of the general 
situation; quotations will be distinguished by a B (for Born). 


If we have just one variable (N=1), there is only one, 
(2)  2N—1=1, 


first integral, which will directly define the energy. Assuming stability, the 
motion is periodic, and we introduce an angular variable w, of which both x and 
p are periodic functions (period 22) 

(3) w=vi+ WwW; 


during one period, w increases by one unit. 


The corresponding momentum is the action variable 


(4) J=$pdx 
taken for a full period. The energy depends only upon / 
(5) E=E(j), 
and the equation of motion yields 
dw,» 0E 
.) eS ed 


where » is the classical frequency. The action J was chosen for quantization, 
and the old quantized mechanics used the assumption 


(7) Wii 1 He n an integer, 
with a quantum frequency 

AE NE 
(8) Ye Ayah. 
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The similarity of equations (6) and (8) was the basis of Bour's correspondence 
principle. On the other hand, it was proven that any slow change of a parameter 
in the energy formula would leave J (and hence m) unchanged. This is the prin- 
ciple of adiabatic invariance (B., Chap. II, §§1—12) which applies as long as 
the frequency v does not become zero. 

These formulas can be easily generalized for a stable system with separable 
variables*. This is the case when the independent variables x, ..., %,% can be 
chosen in such a way that each pair xx, px obeys a separate equation of motion 


and is periodic with a frequency », [the period being t= _}. We have (3), 


Vr 
(4), (6), (7) as before (each one being given a subscript k), and the total energy, 
instead of equation (5), is given by an expression 


(9) E=E(ji, Je,---s Jn) 


containing only the action variables J. Such a stable system is called multi- 
periodic. If we consider a quantity represented by a function of the variables 
X,, Pp, 1t can be rewritten as a function of the angular variables w, 


F(W,, Wa, -«+, Wn)» 


and it exhibits period 1 in each of the wx variables. These variables themselves 
are linear functions of time (equation (3)) and the function F can be expanded 
in a multiple Fourier series 


(10) F = > C erm (m-w) = S (6 ap hE Mauss 6 eek Fo) chat 
m 


m ™ 


where m represents a set of integers m,, mg, ..., M™);, and (m-w) or (m-y) are 
scalar products 
(m-w) = > m, + W, (mv) = >) (my, + v,) . 
k k 


The system is non-degenerate if no relation 
(11) (m-v)=0,  ™m integers, 


can be found; otherwise it is called degenerate. In such a case, the number of 
independent angular variables is reduced. If we have m relations (11), the number 
of independent frequencies is N —n. The system is completely degenerate when 
there remains only one angular variable and one frequency y. This is the case 
for instance for Kepler ellipses or for Lissajoux figures. 


The condition (11) for non-degeneracy is very restrictive. It means that all 
frequencies are incommensurable, so that there exists no zero combination 


* Born, M.: Chapter II, paragraphs 14—17. Beware of the fact that Born uses 
different notations 


| Born here 
f : | 
requency |  % Vp 
period | Wp Th 


integer Pe aa Mp 
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_ frequency among them. On the other hand, some quasi-periods T, (POINCARE) 
can be found, that satisfy relations 


(12) Tq = My Ty + Ey = My T. + Eg -+* = My Ty + En 


with errors &, é,..., €y as small as desired. If all the ¢ were zero for a certain 
choice of the m,’s, the t would really be a common period and the system would 
_ completely degenerate, since equation (12) could be written 


(13) BIS ed eee ee © 
q 


, 


M, Ms, mn T. 


which represents N —1 rational relations between the original periods vx. 


In the old quantized mechanics, the J action variables corresponding to the 
w’s were quantized, according to equation (7). These J’s were adiabatic invariants 
for any slow change of a parameter in the equations, so long as the non-degeneracy 
was maintained and no condition (11) would occur between the frequencies slowly 
modified by the change of the parameter. This condition of no degeneracy is so 
restrictive as to make the imvariance very seldom applicable. 


It may however happen that the change of the parameter does not establish 
any coupling between the separate variables, instead maintaining complete inde- 
pendence of their equations. If such were the case, the occasional degeneracy 
would no longer be an obstacle. 


Coming back to our equations with separated variables, we note that the 
system has 2N —1 integrals, namely the N action variables Je (which define 
the energy) and the N —1 phase differences (w, 9 — w»). 


These are the constants considered in Section 1. Their respective importance 
and their general behavior will be discussed later. 


4. Non separable variables. Hamilton-Jacobi procedure 


When variables cannot be separated, the problem becomes much more in- 
volved. One may try a formal solution, based on Hamilton-Jacobi equations, 
or one will wish to be more practical and use a method of approximation. 


We restrict our discussion to the case of stable problems, when no particle 
may go to infinity. 

The Hamilton-Jacobi equations yield a possibility of finding angular variables 
w, and the corresponding action variables J, which formally satisfy equations 
very similar to those of the preceding section. These variables are often called 
Delaunay variables (see Lanczos). A typical feature results from the fact that 
this procedure now involves a “‘contact transformation” instead of a “‘point 
transformation”’. The Poincaré theorem of Section 2 will play a distinctive role 
in emphasizing the difficulties of this general theory. The classical application 
of the Hamilton-Jacobi method is given by Born (B., §§5—8, 15), and the 
discussion proves that difficulties arise whenever a degeneracy condition (11) 
obtains. This means that the definition of the w, J variables may be materially 
altered at points of degeneracy. 
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Such a restriction was easily eliminated in a problem with separated variables, 
because there was no coupling possible between different variables. Here one 
must admit that this sort of coupling does exist in general, and that occasional 
degeneracy must cause trouble. 

At any rate, the Hamilton-Jacobi procedure has never been a way to solve 
a practical problem. It is very useful for general theorems and discussions, but 
not for actual solution. 


5. Successive approximations 


The only available method for the solution of a special problem is the method 
of successive approximations (B., Sections 18, 34—35 and 40—44). In intro- 
ducing the method, it is necessary to start with the case of degeneracy (B., 
Section 18). We assume an unperturbed energy function given by the Hamiltonian 
H, and a small perturbation 2H, (A small), 


(14) H=H,+AH,. 


The problem is supposed to have been already solved for the unperturbed 
Hamiltonian H,, and it may yield a partly degenerate solution. The small 
perturbation is usually enough to remove the degeneracy. For instance, the 
initial problem, being » times degenerate, depends only upon N —w angular 
variables (see Section 3, equation (11)) and obtains that many actual frequencies, 
while ” frequencies are zero. When the degeneracy is lifted, the N —» original 
frequencies are modified by an amount proportional to A, while the 7 additional 
frequencies take small non-zero values, of the order of A. Each perturbation 
function H, will react differently on the m degenerate degrees of freedom, and 
result in a different choice for these additional w and J variables. Hence we 
see directly how a point of degeneracy will be a point for possible sudden changes 
in the definition of angular and action variables. 


Let us consider a non-degenerate system characterized by a perturbation 
term AH,, and the corresponding w, J variables. Then, let us decrease 4 to 
zero and produce degeneracy; finally we introduce another perturbation 4’ Hj, 
which will result in a different choice of angular variables w’ and action variables 
J’. The original J’s were adiabatic invariants for a slow change of A; for A=O 
we have a sudden change in the J values, and the new /J’s are adiabatic invariants 
for a slow variation of 1’. 


A well known example in atomic models is represented by an atom first 
perturbed by a magnetic field (ZEEMAN effect, B., §34) and then by an electric 
field (STARK, B., § 35). 

To be more specific, let us assume a system with 2 degrees of freedom which 
both yield the same frequency 7 in the unperturbed problem. This may be 
a Kepler ellipse where y and # have the same frequency v9. The same ellipse 
may also be represented in x, y coordinates, and they have again the same 
frequency 7%). The first perturbation AH, may separate 7 and %, giving them 
frequencies y, and vg, both slightly different from 9; the difference »,—v, re- 
presents the frequency which is zero in the degenerate case and becomes a small 
quantity of the order of A in the perturbed problem. The second perturbation 
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i’ H' may be a field in the x direction, acting differently on x or y and producing 


frequencies y, and y,. In both cases, action variables are easily found. The 
transformation AH,’ H; exhibits 


continuity of the frequencies See = a 
(15) 
discontinuity in the action variables J,, J,== J,, J,. 


The formulas (6) for the frequencies are always valid despite the discontinuous 
change in the J variables: 


oH, oH oH. 
2 tA =, +A 
Soy, as, : ej, 
oH 
Vo =e Vo mila A By 
At) oH” 
Ve — Vo ip ae ay, 
x 
Py=YytHN a 3 
y 


We have oversimplified the problem in order to show more clearly the fundamental 
mechanism involved. Let us generalize and assume a degeneracy of the type 


(17) My V1 — My Vo =0, M,>0O m,> 0 integers, 


hence with the periods 


TT = M Us TT. = MT 
1 1 1 

Tv — — = 
v My, Vy My Vo 


The action integral J at degeneracy is related to J, J, before degeneracy by the 
formula 


(18) T= (dy t pdm) =— hts Ie 


This relation shows the J discontinuity. 


6. Approximations for non-degenerate systems 


The preceding problem represents actually the general case. When we start 
from a non-degenerate system and modify even very slightly any parameter, 
we immediately pass through a very large number of degenerate situations. 
Non-degeneracy requires a set of incommensurable frequencies. We already 
noticed, in Section 3, that adiabatic invariance of the action variables could 
be proven only so long as non-degeneracy was obtained or separation of variables 
maintained, preventing any coupling between different groups of variables. In 
the general case, we cannot assume complete separation, and the slightest per- 
turbation destroys the non-degeneracy, hence the definition of angle and action 
variables becomes hopeless. It may have to be changed at each point of degen- 
eracy, and these points lie infinitely close to each other. This essential difficulty 
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is easily recognized in the general discussion sketched by Born (B., § 41). Assum- 
ing a Hamiltonian of the type 

(19) H=H,+AH,4+77H,... Asmall 


where H, is the unperturbed Hamiltonian and H,, Hy, ... are perturbation func- 
tions, the resolution requires the formation of a transformation function 


(20) S. ==, Sot: ASy iA aieae 
and the successive terms S,, S, are found to be of the general type (B., p. 291) 


(21) S,=>) x 7 Am 5 ane) (with no constant term), 
Mt (M:V% 


a formula where (m - v9) and (m - w) represent scalar products defined for equation 
(10), while the A,, coefficient is computed from H,, H, ... perturbation functions. 
If the system is originally non-degenerate, then no denominator (m - v9) can be 
identically zero, and BoRN assumes that the expansion is convergent. The case 
is not that simple since a great many of them will be as nearly zero as desired 
(or rather, as not desired). 

These extremely small denominators may make the sum divergent; and we 
must add that the smallest perturbation would make some of these (m - v) actually 
zero and S, infinite. The perturbed action variables J are defined in terms of 
the S function and thus open to any kind of irregular behavior and discontinuities. 
The frequencies are again continuous. This situation will help us understand 
the exact meaning of the great Poincaré theorem. 

In the preceding discussion, it was assumed that degeneracy or non-degeneracy 
was directly controlled by the nature of the Hamiltonian. Let us remember 
that initial conditions of the motion also play a decisive role; for instance, the 
frequency of a physical pendulum is a function of amplitude; if we have a pendu- 
lum oscillating along x and y, with different non-linear restoring forces in both 
directions, we easily see that the frequencies », and y, can be made commensurable 
or incommensurable by very small changes of the initial conditions. Hence initial 
conditions also control degeneracy. This point will be of importance for later 
discussion. 


7. Poincaré’s great theorem again 


We must now come back to PoINcARE’s great theorem, which we quoted 
in Section 2 without examining carefully its meaning and importance. The 
importance of this theorem is well recognized by astronomers, but mathematicians 
or physicists often overlook it; its implications reach very far and shatter 
many current ideas in classical mechanics. Let us first explain the exact statement 
of Poincare. He considers a mechanical system obeying the canonical equations 


(22) dx; = OF dys Pee Foe OF 
dt OY; 7 dt Ox; f 
where %,,...,%;,.-.., %y are the variables and y,,..., y;,..., yy are the corre- 


sponding moments, while F(x;, y;) is the total energy of the system. The problem 
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of Porncar® is the three-body problem of celestial mechanics, and the energy 
F is supposed to be expanded in power series of a certain parameter 


) (23) NO TNO a ee 


This represents a typical case of perturbation, and the yw expansion of equation 
(23) corresponds to the A expansion of our previous equation (19), since H or F 
both represent the total energy. What is different in PoIncAR#’s case is the 
choice of the %;, y; variables. 

In the unperturbed problem, it is assumed that /, depends upon the N 
_ variables x; but not upon the y;’s. The perturbation functions F,, Ff... depend 
upon both x; and y,; variables. Furthermore, the whole problem is multiperiodic 
(hence stable), and here is how it happens: 


The mechanical equations (22) give, in the unperturbed problem, 
(24) a =k er constant 3 y,=— a t+ Yo, 


and the original coordinates of the system (PorncarE, loc. cit., pp. 22—32) are 
multiperiodic functions of the y;’s. Hence the y; (momentum variables) play the 
role of our previous w,; (angular variables) while the x; variables replace the 
previous action variables. Finally, the situation is exactly similar to the previous 
problem, with just an interchange of variables and momenta, and such an inter- 
change is of no consequence whatsoever. 


Accepting y,; as representing an angular variable w;, we see that 
(25) eres 


is the corresponding frequency. (There is, however, a different use of the 27 
factors in the treatments of BORN or POINCARE, but this again is immaterial.) 

The system is multiperiodic with respect to the y; variables, and this is also 
true for the /,, £, perturbation terms. 

With these assumptions, PoINCARE explains the precise meaning of his 
theorem: 

Let g be an analytic and uniform function of x,y, which must also be 
periodic with respect to the y variables. Let this function g be supposed to be 
analytic and uniform for all values of the y’s, for small values of w, and for x’s 
contained in a certain domain D. Under such conditions, the function g can be 
expanded 
(26) P= Poth Ate yet- 


with 9, %, Y2, -.. uniform in x, y and periodic in the y’s. PoINcARE then states: 
Except for some special cases to be discussed later, such a function cannot be an 
integral of the system of equations (22). Only the total energy F has such properties. 

This theorem invalidates the system of approximations developed by Born, 
who tried to find expansions of type (26) for his action variables J in a perturbed 
problem by using a transformation function S also expanded in a similar fashion. 
The S and J’s cannot be well behaved analytic functions. This is the very 
remarkable consequence of POINCARE’s theorem. 
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We do not intend to repeat here the details of the proof given by POINCARE, 
but we shall stress a, few points of importance, since they explain clearly that 
all the trouble comes from points of degeneracy, which also give difficulties in 
Born’s discussion. This old theorem of PorNcARE proves that these difficulties 
are much more serious than BorN imagined and that they build up an almost 
impassable roadblock. 


8. The role of degeneracy conditions in Poincaré’s theorem 

The clue, in Pomncaré’s discussion, is to consider, for each stage of approxi- 
mation, the behavior of the function y at the points of degeneracy. These points 
are defined by conditions 
(27) >» oe = >'m,v,;=0, m, integers, 
according to (25). And Porncaré£ notes that, in any domain of the x,’s, even 
the smallest domain you can imagine, there is always an infinity of such points 
of commensurability (or degeneracy; see, for instance, POINCARE, p. 239, bottom 
of the page). 

The perturbation function & is expanded in multiple Fourier series (POINCARE, 
Pp. 236—238), 
(28) 1 ei a 2g ida ko 


a formula similar to our equation (10). The general case refers to a perturbation 
fF, that may contain all the B,, terms. In such a case, the proof shows that the 
function gm cannot exist. When some of the B coefficients are zero, the discussion 
runs differently (POINCARE, pp. 240—245) since there is a possibility of secular 
perturbation. For certain classes of problems, it is still possible to prove the 
great POINCARE theorem, but for another group of problems, there appears a 
possibility of finding a certain limited number of well behaved first integrals, 
in addition to the energy integral. These are the exceptional cases announced 
by PoINcarE in the statement of his theorem. These exceptional cases practically 
correspond to problems where a limited number of sets of variables can be 
separated away. 

A typical example corresponds to a problem with conservation of the total 
momentum or moment of momentum. 

There are many other interesting remarks to be found in Porncar&’s remark- 
able book, and they apply not only to celestial mechanics but to general mechanical 
problems. In his third volume (Chap. 28—31), PorncarEé defines periodic 
solutions of the second kind, which correspond to another example of degeneracy: 
some trajectories exhibit a period t, but other trajectories, very close to the first 
ones, may not be periodic at all or may have a much longer period 7 7, a multiple 
of the original tr. An example of these conditions was found in the problem of 
degeneracy discussed in equations (17), (18). 

We can say, in conclusion of these two sections, that the very fundamental 
properties discovered by Poincare in his great theorem are provoked by degen- 
eracy conditions when the Hamiltonian gives rise to coupling between the different 
sets of variables. 
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9. Degeneracy conditions and the possibility of finding a Hamilton-Jacobi 
transformation function 


We shall now generalize the preceding remarks and show that the Hamilton- 
Jacobi transformation function S itself may, in general, be very difficult or 
even impossible to obtain. We shall prove that the Hamilton-Jacobi method 
works wonderfully for separated variables but usually becomes impracticable 
» when the variables are not separable. 


We start from the following formulas®, which can be taken as a definition 


| of the transformation fanction S (S was called F in my book on tensors) of 
~ equation (20): 


os os 
(29) Oak =p,; hence dS= Od 19, = > Pp dp. 
We thus can write 
Tk qk qk 
(30) SSP Ad pe) Ong, ct En ae 
1k 0 9k 0 9k 0 


The kinetic energy F,;, was called J in my book. The integral is taken 
along a natural trajectory running from an initial point ¢,) at t=0 to a final 
point g, at ¢. This integral is completely defined when we specify the total 
energy E of the system under consideration; once we know E, and we give the 
position q, (¢) at time ¢, then the potential energy V(q,) is known, and the difference 
E—V gives F,;,. This, in turn, determines the velocity along the trajectory. 
The total time interval required for a run from qd,» to g, need not be specified 
separately. For further discussion, see my book® on tensors. 

We wish to consider a stable system, with no particle running to infinity; 
such a system will be multiperiodic with no secular term. We thus assume a 
multipertodic solution, where all quantities (£,;,, for instance) can be expanded 
in multiple Fourier series of type (10): 


(31) Po ee ae eae 


m 


The integration (30) yields® 


mv 


The question now is what happens when the denominator (m - v) 1s zero (complete 
degeneracy) or even simply very small (as small as desired, approximate degener- 


acy) ? 


5 BRILLOUIN, L.: Les tenseurs en mécanique et en élasticité. Paris: Masson 1938; 
see Chap. VIII, paragraphs 12—16, especially p.171. These definitions were system- 
atically used by L. BottzMANN in his lectures on mechanics. 

6 Let us specify clearly the difference between equations (21) and (32). Equation 
(21) contains the “‘unperturbed”’ frequencies », corresponding to the Hy Hamiltonian. 
In the other equation (32) the formula works with the final » frequencies corresponding 
to the full problem ruled by the complete Hamiltonian H, including the so-called 
perturbation. Either one of these formulas can be used and may be more convenient 
for discussion. 
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It is at once obvious that such conditions are critical and may mean diver- 
gence of the multiple Fourier series, hence impossiblity of building up the func- 
tion S. We must here consider two distinct problems: separated or non-separable 
variables. 

When variables can be separated, there is absolutely no coupling between a 
qg, variable and another g, one. There exist no cross-terms either in the kinetic — 
or in the potential energy. As a result, there must be no cross-terms in the 
Fourier expansion (31) and no troublesome terms in equation (32): 


(33). 3S = Sy (Oy) Sa (Ge) ae ops) Sn (Ge) + pa Ny Ji, (m integers). 


The transformation function is a sum of separate functions of each variable, 
but a function S,(g,) is only defined modulo /, for the simple reason that we 
may use, from the initial g,9 to the final g, position, an arbitrary trajectory 
making m, loops (each of frequency »,). Each S, can be separately expanded 
in a single Fourier series 


(34) Ge ys us C ere t Mg Kt 2 m1 ME Wok “he Ny is 
k / 


My mK 


where we recognize at first sight the absence of any term in 
(m-v\ t= (>) m»,) t. 


The problem may be non-linear in each separate variable, but cross-terms 
are not allowed. As noted in Section 3, the total energy is a function of the 
J, action variables (equation (9)), which are themselves analytic. 


If the variables cannot be separated, the situation is completely different. 
POINCARE’S theorem proves that in the most general case the j, quantities 
cannot be analytic, and our formula (32) shows that the S function diverges. The 
Hamilton-Jacobi procedure collapses. This is what happens if we make no special 
assumption about the Hamiltonian. 


Many authors, including M. Born, overlooked the non-analyticity of the /’s 
in the general case, and assumed that one could still write a formula of type 
(9), giving the energy as a function of the J’s. It would however be completely 
meaningless to express a well behaved expression, the energy, as a function of 
non-analytic quantities, the J’s. The whole method simply does not work in 
the most general problems, when no restriction is used to limit the values of the 
C,, coefficients. 


10. Sketch of a discussion of the possibilities of convergence 
for non-separated variables 


We found that the Hamilton-Jacobi method works perfectly for separated 
variables, but does not work, in general, for non-separated variables: where is 
the exact limit between these two extremes? We shall find problems of ‘‘approxi- 


mate separation’ where a solution is still possible, provided certain conditions 
are fulfilled. 
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Let us consider again a perturbation problem, with a Hamiltonian (14) 
| (44) H=H,+AH, 


i where Hy is the unperturbed system with separation of variables and H, per- 
turbation with coupling terms between variables. 


| We may discuss this problem according to the perturbation method of Sections 
5 and 6. We look for an expansion of our function S (see equation (20)) 


(20) StS AeSHet eS eh 


where Sy corresponds to separation and is represented by equation (33), while 
S;, Sy... contain the whole series of terms as in equation (21): 


nk Zales 2ri(m-w) 
4) P ean 
where v) represents the unpeturbred frequencies®. We discussed such formulas 
at the end of Section 6, and we followed Born’s distinction between (A) non- 
degenerate problems, where no (m - 79) =O relation can be found, (B) degenerate 
problems, with some (m - v9) =0 relations 


This distinction, however, looks rather artificial. 


In a non-degenerate situation (A), we still have an infinite number of in- 
finitely small (m - v9) denominators and the corresponding numerators A,, must 
be extremely small if divergence is to be avoided. In a degenerate case (B), 
as soon as one denominator (m - 9) is zero, there is an infinite number of other 
(pm-¥9) zero denominators (f an integer), and all the corresponding numerators 
must be zero. The distinction is really superficial. 

The actual problem is concerned with the convergence of the series of coeffi- 
cients —4™ , which should be discussed. The multiple Fourier series for S,, 


M+ Vo) 


S, ... are all of the same type (21) where m stands for the set of indices m,, m,..., 


my. These series converge only when the series of the coefficients ey are 

0 
absolutely convergent. If this be the case, the S function can be constructed, and 
the Hamilton-Jacobi method works correctly. 


Am 


mM: 
cannot be proven, the multiple Fourier series for S diverges, the a 
Jacobi method collapses and we are in the “general situation” described by 
PoincaRE. In order to achieve convergence, it is obviously necessary that the 
A,, coefficient become zero as (m - v9)” with an exponent «21; it is also required 
that the A,,’s themselves decrease fast enough for very large values of the 
M1, My... My indices. When such conditions are fulfilled, the Hamilton-Jacobi 
method can be used, and it extends far beyond the case of separated variables. 
These examples may be called problems of approximate separation. 


If the absolute convergence of the multiple series of the coefficients 7 
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11. Discussion of a simple example with two variables; 
degeneracy means instability or resonance 

Degeneracy occurs when one frequency is zero, or when a condition of com- 
mensurability exists between the original frequencies. If the zero frequency 
corresponds to one of the variables of the motion, then it usually means a situation 
of instability. Such is the case for a rigid pendulum oscillating with very large 
amplitudes and reaching the top with zero velocity, whence it may fall back 
either to the right or to the left. The condition of commensurability will be shown 
to correspond to resonance between two or more degrees of freedom. 
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Fig. 1 


Let us sketch the problem for a system with two variables, hence two un- 
perturbed frequencies 9, and »9,. Figure 1 visualizes the general situation: in 
a plane, we use the m7, as abscissae and m9. as ordinates (7, m, integers). 
Degeneracy holds when one of the points of the rectangular lattice happens to 
fall upon the second diagonal DOD". 


(35) (M+ VY) =M,%1 + My%qQ=O OF My %qy= — Myo. 


On our diagram, this is almost the case for m,=4, m,—=—5. The numerators 
Aimm, Must be zero along the diagonal line DOD’. When the variables are 
separated, the A,,.,,, are different from zero only along the axis 1 and the axis 2 
and remain zero everywhere else in the plane. ; 

Approximate separation holds when the Anm, become zero along DOD’ and 
decrease suitably at large distance. ; 

Let us now consider the case of divergence and try to understand the physical 
meaning of the difficulty encountered. Condition (35) means resonance: ihe m 
harmonic of yg; coincides with the m, harmonic of 15. This Say ie 


cause trouble, since it is well known that resonance leads to infinite terms when 
no damping is available. 
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In many problems with no damping terms (neither positive nor negative), 
it happens that coupling between two oscillations at resonance may destroy 
resonance and give rise to two distinct frequencies. This is, in a different 
language, the well known result that any perturbation (coupling) applied upon 
a degenerate system (having internal resonance) will lift the degeneracy (destroy 
resonance). ; 


When two periods m7), and —m,¥79, happen to be almost equal in the 
unperturbed problem, there is no complication so long as the variables remain 
completely separated and uncoupled. But when we introduce intercoupling 
through the AH, perturbation term (equation (14)), we must be prepared for 
the possibility of a large finite energy transfer from v9, to v9. through internal 
resonance, even when the AH, coupling is extremely small. This is another 
way of explaining the difficulties discussed in Section 5. It would be interesting 
to discuss carefully what happens in a variety of cases of internal resonance, 
for different types of non-linear problems with non-linear coupling, and to dis- 
cover progressively the type of interactions able to take place. This general 
program is intended for further investigation. 


Let us now refer again to Note 6. We find difficulties when there appears 
some resonance between the unperturbed frequencies 91,799, ---, on, but we do 
not know whether any such resonance can take place between the final frequencies 
Y,,Vg,---,Y¥y, defined in Section 9. There are indications that exact resonance 
between the final (actual) frequencies » may not possibly occur, and this again 
is a point requiring careful discussion. 


12. Some general conclusions. Determinism versus statistical mechanics 


We discussed problems of classical mechanics, assuming that the model under 
consideration was stable, with no part of it escaping to infinity. The motion 
was considered as represented by a multiperiodic solution. In such a problem, 
the Hamilton-Jacobi method can be used only when the expansion of the trans- 
formation function S is convergent. This condition is automatically satisfied for 
separated variables, and can be satisfied for ‘‘almost separated” variables. In 
the general case the S expansion diverges, and PoINCARE’s great theorem defi- 
nitely proves that the whole method is worthless. 

The problem can be examined from two different angles, giving two distinct 
and even opposing points of view: the pure mathematical aspect first, and the 
viewpoint of a physicist, second. 

The preceding considerations were based on a purely mathematical discussion. 
It was assumed that a trajectory could be defined by “‘given”’ initial conditions, 
the values of the coordinates and momenta at time ‘=0 being exactly known, 
with no possible error. It was also assumed that the laws of motion were abso- 
lutely correct, with no deviation whatsoever. These two assumptions led to the 
consideration of a single trajectory considered as an infinitely thin mathematical 
curve. The periods of the multiperiodic motion were supposed to be exactly 
known, and questions were raised about their commensurability. The answers 
to the problem were distinguished as ‘‘general’” or ‘“‘degenerate’”’ whether the 
periods were incommensurable or commensurable. 
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To all these definitions, the pAysicist raises some very serious objections. | 
His position is the following: 

A. Initial conditions are not ‘given’; they must be measured and observed. 
This means that the coordinates x; and momenta #; are known within errors 
Ax, and Af;. 

B. The laws of motion are no better known. They result from experiments, 
and they summarize a large number of previous observations. The simple laws 
selected for this representation must be considered only as a first approximation. 
The famous laws of NEWTON’s mechanics are no exception, and the proof of 
this statement has been repeatedly given: these laws had to be amended and 
corrected in many instances, for Relativity, for Quantum Theory, etc. ... 

From these premises it is easy to draw some obvious conclusions, which were 
already stated by E. BoreL, M. Born and the present author (see the references 
in footnotes 7 and 8) and which were sketched briefly at the end of Section 1 
of this paper. 

a) The notion of an irrational number is unphysical. It is impossible to prove 
that a physical quantity is rational or irrational. 

b) The notion of commensurability is unphysical, when the quantities to be 
compared are known only to within some experimental errors. 

c) It is impossible to study the properties of one single (mathematical) tra- 
jectory. The physicist knows only bundles of trajectories, corresponding to slightly — 
different initial conditions. E. BorEL, for instance, computed that a displacement 
of 1cm, on a mass of 1gram, located somewhere in a not too distant star (say, 
Sirius) would change the gravitational field on the earth by a fraction 101. 
The present author went further® and proved that any information obtained 
from an experiment must be paid for by a corresponding increase of entropy 
in the measuring device: infinite accuracy would cost an infinite amount of 
entropy increase and require an infinite energy! This is absolutely unthinkable. 

d) Let us simplify the problem, and assume that the laws of mechanics are 
rigorous, while experimental errors appear only in the determination of initial 
conditions. In the bundle of trajectories defined by these conditions, some may 
be “non-degenerate” while others may ‘‘degenerate’. The bundle may soon 
explode, be divided in a variety of smaller bundles forging ahead in different 
directions. This is the case for a model corresponding to the kinetic theory of 
gases. BOREL computes that errors of 107° on initial conditions will enable one 
to predict molecular collisions for a split second and no more. It is not only 
“very difficult”, but actually impossible to predict exactly the future behavior 
of such a model. The present considerations lead directly to BoLTzMANN’s 
statistical mechanics and the so-called “ergodic” theorem. 

The problem is of such importance that it requires a few more words of 


explanation. Recalling the results of Section 10, we must distinguish between 
three different possibilities: 


’ Borer, E.: Introduction géométrique & quelques théories physiques, p. 94. Paris: 
Gauthier-Villars 1914. — Born, M.: Phil. Quart. 3, 139 (1953). — Kgl. Danske Viden- 
kab. Selskab., Mat-fys. 30, No. 2 (1955). — J. de Physique 20, 43 (1959). 

8 BRILLOUIN, L.: Information and Control 1, 2 (1957). — Nature 183, 501 (1959). — 
Annals of Physics 3, 243 (1958) 


Shortcomings of the Hamilton-Jacobi Method 93 


1. Separated variables. We must use the variables w, and J, introduced in 
Section 3. The j, retain constant values through the bundle and are defined 
within errors A J, that depend only upon the errors made about initial conditions. 
These Aj, remain constant in course of time. The problem was discussed in 
great detail elsewhere®. The angle variables w, exhibit a very different behavior. 
The error Aw, increases steadily in course of time, and after a while it may 
become larger than 2%, which means that we know no longer what this angle 
may be. 

In such problems, the J, are well behaved functions, but the angles w, are 
not. Of course, predictions about the w, values may be good for a certain period 
of time, right at the start; or even (if initial conditions are very exactly defined) 
the predictions may be correct for a rather long interval of time (as is the case 
in Astronomical problems), but in the long run the uncertainty over the phase 
angles must always become dominant. 


2. Almost separated variables. Similar conclusions hold but for the fact that the 
J, @, variables now result from a contact transformation based on the function 
S instead of a point transformation. 


3. General case (POINCARE). No S function can be constructed, and the total 
energy is the only well behaved function in a system of coordinates where the 
center of mass is at rest. 

For separated or almost separated variables, we obtain determinism at the 
beginning, and after a certain time interval the angles w, are no longer observable, 
but the momenta /, are still well behaved integrals of the motion. 

In the general case, the total energy is the only quantity about which predictions 
can be made. (See the reference in footnote 9, Chap. vi, p. 214.) 

e) We may choose another simplified alternative and assume initial conditions 
to be rigorously “‘given’”’ but fix our attention on the fact that the laws of 
mechanics are only a first approximation. These classical laws correspond to 
an “‘unperturbed”’ Hamilton function Hy, which may be completed by some 
unknown small perturbation 1 H,. What is the role played by this perturbation ? 
This is practically the question answered by POINCARE in the famous theorem 
quoted in Section 2 and 7. In the general case, the total energy comes out as 
the only well behaved analytic solution of the initial conditions. It appears that 
the Poincaré theorem contains the justification of BOLTZMANN’s statistical 
mechanics, which should apply when (and only when) the total energy of the 
system actually remains the only well behaved first integral. 


f) All these approximations must be carefully weighed and discussed. For 
some problems (kinetic theory of matter) they make predictions almost impossible. 
For other types of problems (astronomy) the situation is not so hopeless. The 
fundamental question to be discussed is the stability of the solution for small 
changes in the initial conditions or in the Hamiltonian. This question has been 
generally overlooked. 

g) Classical mechanics assumes that it is reasonable to speak of two tra- 
jectories that may be infinitely close together. If then there is a finite difference 


» BRILLOUIN, L.: Vie, Matiére et Observation, Chap. V. Paris: Albin Michel 
1959. — Information and Control, vol. 2, p. 45. 
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on a certain quantity between these trajectories, it means a real discontinuity 
and a non-analytic function. ei ; 

Wave-mechanics, applied to a stable problem, gives proper y functions, defined © 
by quantum numbers. These functions (or the corresponding semi-classical tra-_ 
jectories) can never be infinitely close; they remain at finite intervals, and all 
quantities normally exhibit finite differences. é : 

The transition from wave mechanics to classical mechanics is obtained by 
the B. W. K. method, but it should be emphasized that this method has been 
discussed only for problems with separated variables. The transition from wave 
mechanics to classical mechanics for non-separated variables is an open question. 
There are definite indications? that in some problems (scattering, collisions) 
the transition is not a smooth one and that geometrical optics, for instance, 
does not result from wave optics even for infinitely small wave lengths. There © 
is a whole field of unsolved problems in this connection. 


The research presented here was done under Contract Nonr 266(56) and was first 
described in a Quarterly Report dated July 31, 1959. 


10 BRILLOUIN, L.: J. of Applied Physics 20, 1110 (1949). 
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